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ABSTRACT 


MLA VCR,  a  computer  program,  was  developed  by  the  Naval  Ocean  Systems 
Center  (NOSC)  for  calculating  the  signal  levels  of  electromagnetic  waves  propagating  in 
a  multilayer  tropospheric  waveguide  environment  over  seawater.  The  program  is  an  ex¬ 
tension  of  the  XWVG  which  is  a  trilincar  ducting  program.  Modifications  of  the  XWVG 
were  carried  out  to  handle  multilayer  tropospheric  ducts. 

A  number  of  modifications  and  improvements  on  the  program  made  over  the  past 
several  years  were  not  documented.  A  detailed  documentation  of  M  LAYER  was  also 
not  available.  The  objective  of  this  study  is  to  develop  a  technical  documentation  for 
M  LAYER  using  the  program  as  baseline.  The  study  aims  to  put  together  the  theoretical 
formulations  (specific  to  M LAYER)  into  a  complete  self-contained  document.  This  is 
to  facilitate  potential  users  with  better  appreciation  of  the  capabilities,  limitations,  ap¬ 
proximations  and  assumptions  used  in  the  mathematical  modelling  techniques.  As  far 
as  possible,  the  same  terminologies  and  functional  variables  used  by  Baumgartner  (*n  the 
XWVG  development)  and  by  Pappert  (in  the  MLAYER  development)  are  adopted  to 
enable  one  to  relate  this  document  to  the  program.  Step-by-step  derivation  of  certain 
equations  was  carried  out  and  checked  for  compatibility  with  the  algorithm  in  the  pro¬ 
gram.  An  in-depth  scrutiny  of  each  program  element  was  also  conducted  and  a  de¬ 
scription  for  each  is  provided.  As  a  result  of  a  detailed  analysis  of  the  respective 
algorithm  in  the  program,  the  documentation  for  the  evaluation  of  the  modal  function 
was  eventually  prepared.  Additional  materials  were  gathered  from  technical  reports  and 
papers  to  supplement  the  development  of  this  document. 

The  MLAYER  supporting  programs  (Microsoft  program  maintenance  utility 
"makefiles")  were  modified  to  enable  the  program  to  run  on  Microsoft  FORTRAN  ver¬ 
sion  5.0.  MLAYER  was  tested  and  ran  successfully  on  Microsoft  FORTRAN  version 
5.0.  and  C  compilers  version  5.0. 
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I.  INTRODUCTION 


Propagation  of  radio  waves  through  the  troposphere  is  strongly  affected  by  the 
complex  nature  of  this  propagation  medium.  The  effects  on  radio  wave  propagation  are 
due  to  changes  in  the  refractive  index  with  height,  absorption  by  the  atmosphere  and  its 
constituent  gases,  attenuation  and  scattering  due  to  climatic  conditions. 

Although  weather  conditions  could  result  in  attenuation  of  radio  wave  transmission 
through  the  absorption  of  energv  by  gases,  water  vapour  or  scattering  by  rain,  certain 
meteorological  conditions  can  lead  to  an  enhancement  of  signal  strength.  This  peculiar 
phenomenon  known  as  ducting  is  due  to  the  inversions  of  the  refractive  index  profile. 
Propagation  in  a  tropospheric  duct  can  result  in  an  unusual  increase  in  range. 

To  predict  the  anomalous  propagation  effects,  knowledge  of  the  characteristics  of 
the  refractive  index  and  its  effects  on  propagated  fields  is  required.  In  the  mathematical 
formulation  of  the  model,  the  concept  of  modified  index  of  refraction  is  used  to  describe 
wave  propagation  over  a  flat  earth. 

In  a  uniform  atmosphere  with  constant  index  of  refraction,  electromagnetic  waves 
travel  in  straight  lines  over  the  curved  surface  of  the  earth.  This  situation  is  unaltered  if 
the  earth  is  considered  as  flat  and  the  rays  curved,  provided  the  relative  curvature  be¬ 
tween  the  earth  and  the  rays  is  preserved.  To  account  for  the  upward  bending  of  rays 
over  the  flat  earth,  the  constant  index  of  refraction  is  transformed  to  a  modified  index 
which  increases  with  height  and  is  given  by 

m(z)  «  n  +  -f  ,  (1) 

where  n  is  the  index  of  refraction,  z  is  the  height  above  the  surface  of  the  earth  and  a  is 
the  radius  of  the  earth.  Equation  1  is  still  valid  in  a  stratified  atmosphere  where  the  index 
of  refraction  is  a  function  of  height.  In  this  case,  n  is  replaced  by  n(z)  in  Equation  1. 

The  modified  index  concept  allows  one  to  transform  a  spherically  stratified  refrac¬ 
tive  structure  into  a  planar  layer  by  preserving  the  relative  curvature  between  the  normal 
of  the  radio- wavefront  and  the  earth’s  surface.  Therefore  wave  propagation  over  a 
spherical  earth  can  be  reduced  to  that  over  a  flat  earth  if  the  index  of  refraction  is  re^ 
placed  by  the  modified  index  of  refraction.  For  z/a  <3  1,  the  angle  of  inclination  of  a  ray 


at  any  point  in  space  and  the  distance  (optical  length)  along  the  ray  path  will  remain 
invariant  in  the  transformation. 

The  flat-earth  approximation  is  adopted  in  M LAYER  [Ref.  1  and  2].  This  ap¬ 
proximation  considerably  simplifies  the  mathematical  solution  of  the  problem.  The  ac¬ 
curacy  of  the  flat-earth  approximation  was  investigated  by  Pekeris  [Ref.  3].  It  was  found 
that  for  a  standard  atmosphere,  the  fractional  error  in  the  height  gain  function  is  ap¬ 
proximately 


2ffv*6  -s/a 

5  Xa 3/2  ' 


(2) 


Equation  2  shows  that  the  fiat-earth  approximation  will  break  down  at  high  frequency 
and  at  great  height  above  the  earth  s  surface.  Tor  example,  a  10  GHz  signal  at  a  height 
of  2000  m  will  give  a  fractional  error  of  approximately  1 14%.  The  work  of  Pekeris  also 
indicates  that  the  flat-earth  approximation  is  valid  to  within  two  percent  for  ranees  up 
to  about  one-half  of  the  radius  of  the  earth  [Ref.  3]  . 

As  typical  values  of  m(z)  differ  from  unity  by  about  300  parts  per  million,  it  is 
sometimes  more  convenient  to  introduce  the  modified  rcfractivity  to  describe  the  phe¬ 
nomena  involved  in  the  formation  of  tropospheric  ducts.  The  modified  rcfractivity  is 
defined  as 

M.(z)  =  [m( r)  -  1]  x  10*  .  (3) 


In  this  program,  the  modified  rcfractivity  profile  of  the  troposphere  is  modelled  as  a  sc¬ 
ries  of  continous  piecewise  linear  segments  (see  Figure  1). 

The  linear  approximation  simplifies  mathematical  manipulations  to  a  large  extent. 
In  general,  any  complicated  profile  can  be  approximated  by  several  linear  segments  and 
the  accuracy  of  the  modelling  can  be  improved  by  increasing  the  number  of  linear  seg¬ 
ments  and  reducing  the  thickness  of  each  layer. 

If  Xf,(z)  is  the  value  of  the  modified  refractivity  in  the  /'*  layer,  then  for  a  linear 
modified  rcfractivity  profile, 


M,(z) 


M  {zj)  + 


«/.U, 

dz 


(2  “  2,) 


zt  £  z  £  r,+,  . 


(4) 


From  Equation  3,  the  square  of  the  modified  index  of  refraction  is 


3 


mf  (;)  2s  1  +  2  x  10  6 Mt  (z)  -  jtjt  ,  (5) 

where  r/(  is  introduced  to  account  for  attenuation  of  the  signal  due  to  the  absorption  by 
atmospheric  gases.  The  derivative  of  in}  (z)  is  given  by 


2  x  10"6 


=  <*/  • 


(6) 


Equations  4.  5  and  6  will  be  used  in  the  mathematical  model  to  represent  the  character¬ 
istics  of  the  propagation  medium. 

\1  LAYER  also  incorporates  in  the  solution  of  the  modal  function,  the  effects  of 
surface  roughness,  atmospheric  absorption  and  the  variation  of  the  complex  index  of 
refraction  of  seawater  as  a  function  of  temperature,  salinity  and  frequency. 

The  program  uses  the  Shellman-Morfitt  complex  root-finding  routine  to  locate  the 
modes  that  propagate  in  the  tropospheric  ducts.  This  routine  attempts  to  find  all  wav¬ 
eguide  modes  with  attenuation  rates  below  a  user -specified  value. 

Overflow  problems  commonly  encountered  in  the  evaluation  of  the  modal  function 
are  avoided  by  performing  numerical  computations  with  extended  complex  arithmetic. 
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II.  MATHEMATICAL  FORMULATION  OF  THE  MODEL 

A.  HORIZONTALLY  POLARIZED  WAVE  PROPAGATION 

In  formulating  the  problem,  the  flat-earth  approximation  with  modified  index  of 
refraction  is  adopted.  The  cylindrical  coordinate  system,  (r,  < f> ,  z) ,  is  used.  The  vertical 
height,  z,  is  measured  from  the  surface  of  the  fiat-earth. 

A  horizontally  polarized  source  may  be  approximated  as  a  radiating  magnetic  di¬ 
pole  oriented  in  the  z  direction  and  located  at  r=0  and  z  -  zT  .  The  magnetic  Hertz 
vector  associated  with  the  dipole  is  given  by 

n  (r,  z)  =  n  (r,  2)  z  ,  (7) 

where  n  (r,  z)  satisfies  the  following  equation: 

V2  n  (r,  z )  +  kW(z)  n (r,  z)  =  —An  {<5(.v)  <5(»  S(z  -  zr) }  ,  (8) 

where  k  is  the  wave  number  in  free  space,  <5(.)  is  the  Dirac  delta  function  and  zT  is  the 
location  of  the  source.  The  right  hand  side  of  Equation  8  corresponds  to  a  point  source 
of  strength  —  4*  located  at  r=0  and  z  —  zT  .  The  time  dependence  e"°'  is  assumed 
throughout  this  thesis. 

The  electric  and  magnetic  field  can  be  expressed  in  terms  of  the  Hertz  vector  as 

E  —  —  y'cuA/pV  x  fl  (9) 

and 

H  =  V  x  V  x  U  ,  (10) 

where  n0  is  the  free  space  permeability.  Solution  to  Equation  8  can  be  represented  as  a 
contour  integral: 

n  (r,  z)  -  (rp)y(/)(:)  ,  (11) 
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where  Hf  (/•/»)  is  the  zeroth-order  Hankcl  function  of  the  second  kind  and  jlp.z)  is  the 
height  gain  function.  The  contour  of  integration,  C,  as  shown  in  f  igure  2  can  be  de¬ 
formed  into  C,  and  C2  ,  also  shown  in  Figure  2. 
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1  he  height  gain  function  in  the  i"'  tropospheric  layer,  denoted  as  f(p, satisfies  the 
following  differentia!  equation: 


2  /  x  P 
mi  fc)  — ~y 

K 


}/(/>. 2)  =  “  ( z  -  zT) 


(12) 


In  the  ground,  the  height  gain  function,  fg(p,  z),  is  given  by  the  solution  of  the  fol¬ 
lowing  equation: 


(13) 


where  nt  is  the  complex  index  of  refraction  of  the  ground  and  is  related  to  the  relative 
permitivity,  r.t.  and  conductivity,  os,  of  the  ground  as 


2 

na 


0„ 

c0(O 


(14) 


The  height  gain  function,  f.(p ,  z),  must  represent  an  outgoing  exponentially  decaying 
wave  as  z  -*  —  oo,  i.c., 

Jj  (/>,*)  =  V'  ,  (15) 


where 

y  -  kjn*  -  p2 


(16) 


and 


P  =  J  •  (17) 

For  the  wave  to  attenuate  into  the  ground,  the  branch  of  the  square  root  in  Equation 
16  should  be  chosen  such  that 


/m(y)  <  0  .  (IS) 

If  the  ground  is  not  dissipative,  the  branch  point  and  the  integration  contour  satisfying 
Equation  18  is  shown  in  Figure  3. 
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Note  that  the  A,  's  and  5/s  arc  different  from  those  in  Reference  2  where  the  height  gain 
function  is  expressed  as  5,  (A,  +  A,  k3) .  Equation  25  corresponds  to  that  implemented  in 
the  program. 
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Figure  3.  Branch  of  %!nj  —  / ? 2 

1  lie  constants  of  proportionality  of  liquations  26  and  27  arc  of  no  consequence  as 
they  all'ect  only  B,(p)  and  A,(p)  of  liquation  25  but  not  J,(p,  z).  1  lie  height  gain  functions, 
J'\p,z),  used  in  the  computations  of  the  mode  sum  and  path  loss  arc  normalized  as  in 
liquation  77.  1  lie  solution  of  Stokes’  equation  is  discussed  i  Chapter  IV. 
l  or  the  highest  layer,  i=  1,  the  height  gain  function  is  given  by 

J',(p,z)  —  B,  (p)  h2  (<//)  ,  (28) 

where  /i, (q,)  is  the  modified  llankcl  function  of  older  one-third  and  is  proprotional  to 
the  Airy  function  as 

M‘ 7/)  ~  Ai(-q,e~J2^)  .  (29) 
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As  in  liquations  2(>  and  2~.  the  constant  of  proportionality  of  liquation  29  it;  also  omit¬ 
ted  hi  the  program. 

Liquation  25  represents  the  solution  ol'fiji.  r)  for  every  layer  except  the  one  in  which 
the  source  is  located.  In  the  layer  where  the  source  is  located,  the  height  gain  function 
is  given  by 

/( ./'•-)  =  (/>)  [.1;  (p)  k)  (t/;)  +  /<;  (<?;)]  +  fis\p.(]i)  ,  (30) 

where  /'■•(/».< /)  is  the  particular  solution  of  the  inhomogenous  liquation  12.  Therefore 
the  complete  solutions  of/(/».  r)  are  given  by 

fi  (/».  -M  =  By  (/))  [.?:  ( / » 1 A j  (</,)  +  k:  (</.)]  +  q,)  oiU  ;  /  =  1  to  I  -  1  (31 ) 

and 

J]  (/».  :)  =  B.  (p)h:  (c/,)  +fy'\p.  <j;)  <i;A/  ;  i  =  1  .  (32) 

where 
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B;  [•  (;  1  >1+7  ('V#.  <lu+\ >  <\w 

=  -0+1  ^l^/.-i-!  ./-H  )  +  ^:l</;-j-l,i-H^] 


+/  (/7- '/;+!./+ l\+i.' 


M 


and 


^  [•■!/  A-',  (</W+i)  +  A-'2  (^„+1):  [-^T  ]2/3  “/ 


-j-  U'{s\p-  </,,•+ 1>1  <Vv 


—  ^(+1  I>0-H  ^  i^/i+l.i+l)  4"  ^  2W:+1./+1^ 


+  “  [/'*  (/7-  ‘/,-rl.i+O]  ‘0+1,  U 


I  rom  liquation  37. 

*  l  <  (li,i+ 1 )  fy-'h  +  k 2 ( 1  I I'l 

—  ^l(</;+j.!+l)A+l  '0+1  “  ^'2(^+1.  (+!^/+! 


— /'  (/*•  </<+!. (+|)  ‘0+1..\/  J  '/i,  (+l)  ‘O'/ 

From  liquation  38. 

A',^./+|)(-^)2/-  MM/  +  ^(<//.t4|)(^)2%4 


(•V 


^iWf+l,  i+l)^  ft(+)  )  ’  “/• 
—  'M'/.+  i,  i+i)^  a;+1  ^  a 


/+]  "/+ 1  '0+1 


/+!  ^0+1 


_  I/'  M  <7;+!.  j+l)]  jj+l.M 


-4-i /V<7,/+,)]<M  . 


(38) 


(39) 


(40) 


II 


For  i  =  1-1. 


and 


[- (O’/-!. +/  (/’• 

=  /?,  h2  [<-hi)  +J(  \l>.  C ('/./)t)/U 


(41) 


+  77  [/'"(/>. 


U 


(42) 


((>■  9/  i)¥i  m 

v  -/  /  i~ 

From  liquation  41, 

A)((7/_|  /)  /-)/_!•!/_!  + 

=f<\f,-l]ijh)i.M~f(  \i>’  cll-),l^l-\M  . 

From  Equation  42. 

k'\Uli-\ y  ^/- 1  ^/-i  -4/-]  +  1./)^  ^  a/-i  #/-i 


(43) 


=  [/(  IP' <7/,/)](^/,.u  ”  t/(  . 


In  the  above  equations.  <7  *  is  defined  as 

<7u-  =  9;  (’  =  +a)  •  (45) 

Equations  36,  39,  40,  43  and  44  can  be  written  in  matrix  notation  as: 

A  |  =  'A  ,  (46) 


12 


13 


14 


(55) 


^  ^72 .2  ^ 

j*i  •  (56) 

^3.2  =  k’2{l]u)^  *TT-  ^  '  «1  .  (57) 

A.u  =  -A',^)^  y/3a2  ,  (5S) 

A3  j  -  A  2^/2.;^  "JT  ^  a2  •  (59) 

A2/-2.2/-3  =  *l(tf/-l./)  •  (60) 

&21-2.21—2  =  A'2(<7/_i./)  i  (9 1 ) 

^2/-2.2/-l  =  ~  A2(<7/,/)  «  (62) 

^2/— 1,2/— 3  =  A'i(?/_ti/)^  R/_]  ^  ‘  «/_]  ,  (63) 

^2/— 1,2/— 2  =  A'2(<7/_1,/)^  k^_,  ^  1  al- 1  (6-1) 
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and 


A’/— 1,2/— 1  — 


-*■*»»(  i),n 


(65) 


The  elements  of  the  column  vector,  0  ,  are 


/' i  =  9i,i)6ia/  ~  77  [7'5’l/h  9i,i)]‘s 


I _ r  r(s\ 


1A/ 


(66) 


C'2  ~  /(  \/6  92, 2)^2 A/  ~ /'  \P<  9i,2)6l  A/  • 


6), 


(67) 


^3  =  77(/(5\p,  92,2)]‘?2A/-  jr[/'(°(P-  9u']6 


1  A/ 


(6$) 


•A/-!  ~  /(  (/’•  9/,/)6/t/— /(  V'*  9/-1, /)•’/- 1, A/ 


(69) 


and 


0/  = 


C2 


l/V  9/./)]6/a/  ~  77  l/'V  9/-,. /)]<>/-, .a/ 


(70) 


To  determine  the  modes  that  propagate  in  the  tropospheric  ducts,  the  roots  of 


II  A  ||  -0  (71) 

have  to  be  located.  This  is  discussed  in  Chapter  III.  The  symbol.  ||  A  ||.  denotes  the  de¬ 
terminant  of  A  . 

Trom  Equation  11,  the  solution  of  Equation  S  in  the  layer  is 

n  (r.z)  =  4-  f  pdpH0W{rp)f  [,>,:)  •  (72) 

*■  Jc 

From  the  theory  of  residues.  Equation  72  can  be  expressed  as 
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n(r.-)  = 


-i^v. 


T~  M',;)  +  *T  pdpllo™  {rp)f(p.z) 


+Jf. 


(7?) 


+  -7  I  pdpll0{2)  {rp)f[p,z)  , 


where  b,  (m)  are  the  residues  of  pII0G'{rp)f{p,  z)  evaluated  at  the  poles  of/  (p, :)  in  the 
fourth  quadrant.  The  contour  C,  is  the  quarter-circle  of  infinite  radius  and  C2  is  the 
contour  enclosing  the  branch  cut  of  the  radical  of  Equation  16.  Since  Hf'(rp)  decays 
exponentially  in  the  lower  half  complex  p-plane. 


(74| 


For  small  height,  z.  the  integral  over  contour  C;  represents  the  contribution  by  the  sur¬ 
face  wave  to  the  total  field.  Beyond  the  horizon,  this  field  is  small  compared  to  the  total 
field  and  can  be  neglected.  Therefore  Equation  73  becomes 


Tl(r.  z)  =  -  njYjbj  {»>) 


(75) 


where  b,  (m)  is  given  by 

b-t  {»i)  =  II0  a \rpm)  gm  (z)  gm  (rr)  (76) 

and  g,„  (z)  is  the  normalized  height  gain  function  given  by  [Ref.  2  ] 


where  q(z)  —  q(z)  for  z,  <  z  <  r,.,. 

According  to  the  program  source  code,  the  integral,  II  =  f*fm:[q{z)]dz.  is  evalu¬ 
ated  by  considering  that  /,,[</(r)]  is  a  solution  of  Stokes's  equation, 


d\rm 

ch/2 


+  llfm 


0 


?S) 
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This  leads  to 


d 

clq- 


‘ l ! 


({/,,■  t 


r  ‘  +  (  — r 


dq 


_  f  -  ,  -> ,  /•  ‘JUi  ■">1 

-.A/i  +  -‘/.A/r  +  -  .  ■  , 

</l/  1(1/  (/.y- 


=  f  *  +  2  — 

Jm  '  -  ; 


AA, 


d\L 


lid  [_  rA/4 


r  +  <?/m 


(79) 


=  A,A  . 


Note  that  Equations  7S  and  79  imply  that  the  derivative  of f„  is  continous  at  zT  .  That 
is.  fjq)  does  not  include  the  contribution  from  the  last  term  of  Equations  31  and  32. 
The  integral  of  A,  :[</(c)]  in  the  layer  is  given  by 


where 


(SO) 


(81) 


From  Equation  (SO),  it  can  be  seen  that  the  contribution  of  A/,  to  II  at  z  =  z,  (i  ±  I) 
is  given  by 


M  =  ~~r  yhfm 2 1  *//)  + 
]i  L 


#„,(«,)  1:1 r-1-' 

d<h  U, 


li.t  <1 i- 1 ,/ 

— ~  +  — - 

(I  i  <1  i-\ 


rm2(<iu) 


(82) 


df. 

Bv  makine  use  of  the  fact  that  is  continous  at  the  laver  boundaries, 

dz 


which  implies 


<k- 1'  “77 —  Ul.  1‘7/-|.;0  = 


<W. 


and 


J  — Um  (‘/.-I.P]}2  =  l  —  H  4-r  V 


lUh-  ] 


I  \  j  I  dtli 


Uw.  I | 


Substituting  Equation  So  into  Equation  S2.  A/,  becomes 


M  = 


<hj  ‘h-u 

*4"  ; 

‘!  i  ‘I  ;-\ 


+ 


7T 


Eor  i=  1. 


—  /  Vl.l./m  W|,A  {  ?  CAm  (</],!  )]} 

Vi  V  1  (.  <■  ‘1\  I 


The  intceral.  II ,  is 


u  =  y>; . 

i*l 


The  second  term  of  the  denominator  of  Equation  77  can  be  written  as 
fm 7  [y|(l>]  ./  </;•  \  _  ,/A  “[t/f0-*]  /  d;  </</,  \ 

f*n\  \  dp  J p—pm  *-  pw  y  hii  i  c/p  Jp—pm 

.  X,  Tv((h]  dy  /  _p_  \  2/3  2  p„, 
~  ^  2  pm  c/c/,  V  a>  /  A; 

_  ,  /», '  [y(0)]  /  a  \  2/3 

-  ^  A-2  Iv  </>/, 


( S-4 ) 


(85) 


(86) 


(S7) 


(SS) 


(S9) 


l  or  cases  of  interest.  |  />„  r  |  g>  1  .  The  Ilankcl  function  in  Equation  76  can  be  replaced 
by  its  asymptotic  approximation: 
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Therefore  liquation  75  becomes 


To  allow  for  spherical  spreading,  the  r-12  term  in  liquation  91  can  be  replaced  by 
[u  sin  ( ~  )]  !  1  .  a  being  the  earth's  radius,  liquation  91  becomes 


nir.z)  = 


1_ 


a  sin  (  —  ) 


i!  ’  -1 


'vT'  -i/e 

S’  /  Cm  Sr 


(:)  gm  (z,)c-Jr”r 


(92) 


In  the  cylindrical  coordinate  system,  (r.  >/>.  z  ).  the  transverse  electric  field  component, 
is  given  by 


1C 

S.' 


-  <o/<0 


r  n 

£  y 


DilTerentiating  Equation  92  with  respect  to  r  gives 


(93) 


m 


Neglecting  higher  order  terms  of  —  ,  liquation  94  becomes 


in. 

<  r 


o  smf  —  )  _j 


■2  -;-.  J  V'  1  2 
f  ’  /  />#.,  & 


//j  f  >*/•  V-  7  ' 


The  Tree-space  solution  to  Equation  8  is  given  by 


nf5  — 


—;k  .  r 
C  N 


c~i¥r 


T  -l 

r~  +  z~ 


;  r  p  z 


Thus. 


fl/s  jk  ~ikr 


-//•:- 


cr 


i  • 


and  similarly,  the  higher  order  term  of  —  in  Equation  97  is  neglected,  giving 

<  Tl/.S  jk  t_ ji<r 

r  c 


( r 


Normalizing  the  field  to  the  free  space  broadside  field,  the  mode  sum  becomes 


E, 


E-r 


rs 


k'a  sin(  — ) 


!  '2 


V 


Pm  12  Em 


(7  C 

c*m  ^  7  f  L 


The  coherent  mode  sum  in  decibels  is  eiven  bv 


ECMS  =  10  loc 


E 


=  10  l02 


r.s 


2rr 


f _ _ 

l  A2«  sin(  ) 


Pm  2  Sm  (.=  )  Em  (-  /)  *  W 
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(95) 


(96) 


(9‘ 


(98) 


(99) 


(100) 


The  power  sum  or  incoherent  mode  sum  (El MS)  signal  can  also  be  calculated  if  the 
phase  of  each  term  is  ignored: 


ELMS 


10  log 


ir.r 


u: 


a  sin 


Si, 


1,2  £»,(')  &»|(‘ 7-)  c 


J 


(MU) 
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The  path  lossses  are 


Coherent  path  los-  =  32.45  +  2o  log  r  +  20  loci'—  ECMS 


and 

Incoherent  path  loss  =  32.45  +  20  log  r  +  20  log  f—  TIMS 

where  r  is  the  range  in  kilometers  and  f  is  the  frequency  in  Mil/. 

B.  VERTICALLY  POLARIZED  WAVE  PROPAGATION 

A  vertically  polarized  wave  may  be  treated  as  due  to  a  radiating  clectr 
entcd  in  the  z  direction  and  located  a*  r=  o  and  :  —  . 

The  electric  field  and  magnetic  field  can  be  approximated  by  [  Ref.  lj 

//  =  jf. Qi.aV  x  II 

and 


/  •;  =  v  x  v  x  n  , 

where  IT  is  the  electric  Hertz  vector. 

Calculation  of  the  electric  Hertz  vector  is  identical  to  that  of  the  nu 
vector  except  for  the  following  modifications: 

1.  The  height  gain  function  satisfies  the  following  differential  equation: 


l  -v  L . 

2.  Continuity  condition  at  i  =  1  becomes 

0)  =  0; 


1  - 

-  -tr 


This  icads  to 


2  .'3 


y\k  d£/l.l 


~j'i 


/cmm 

■) 


^i(‘/i.i) 


t  m2) 


(inf) 


ic  dipole  ori- 


(104) 


(K>5) 


ignetic  Hertz 


(  I0(q 


(10~> 


(loS) 


*  *s 


^1.2  —  ( 


M^V'/i.i 


-J, - —  AyU/,,) 

ii.  T 


(K»yj 


and 


'-'C  —  jv  -  J\  P  (/'•  ‘/u si  ~  -  [f\  '  *  (/A  <7i ,i )]l’i  u  • 
//“  '  - 

The  transverse  magnetic  field  component.  //  .  is  given  by 

n  '~n 

<.  r 

Trom  Equations  94  to  9S,  it  can  be  seen  that 


(110) 


(111) 


/u 

rs 

Ah/  sim  —  ) 


V 


AW-'?)  *’ 


(112) 


It  follows  from  liquation  1 12  that  ECMS.  I. IMS.  coherent  path  loss  and  incoherent  path 
loss  are  given  by  liquations  loo.  KM.  1  f >2  and  103  respectively. 

C.  DETERMINATION  OF  THE  HEIGHT  GAIN  FUNCTION 

T  he  heieht  gain  function. /(/>. :),  is  given  by 


Ji  {/'•  r)  —  Bj  [.-IjAq  (</.)  4-  As  (</;)]  :  /  —  1  to/  —  I 


(lie) 


or 


f  (p,  z)  =  Bj  lh(q-)  :/=/  .  (114) 

The  coefficients.  A,  and  B„  are  functions  of  p  as  in  Equation  25. 

Determination  of  the  values  of  A,  and  B,  is  dependent  on  the  evaluation  of  the  in¬ 
tegral  (or  summation)  of  Equation  SS.  It  was  found  that  if  the  attenuation  rate  of  a 
particular  mode  is  small  (<  0.1  dB  km),  the  integration  Equation  S8  is  carried  out  from 
the  top  layer  to  the  first  layer  with  the  height  gain  function  at  the  .  p  layer  set  to  a 
certain  finite  value.  This  will  prevent  the  integral  from  blowing  up. 
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On  the  other  hand,  if  the  rate  of  attenuation  is  large  (>  0.1  dB  km),  the  integration 
is  performed  from  the  first  layer  to  the  top  layer  with  the  height  gain  function  at  the  first 
las  er  set  to  a  certain  finite  value.  This  will  ensure  that  the  value  of  integral  will  not  be¬ 
come  too  small. 

Therefore,  for  attenuation  rates  greater  than  0.1  dB  km.  A,  's  and  B,  s  arc  deter¬ 
mined  successively  starting  from  A,  and  Bt  up  to  A,  and  B, .  At  i  =  1  and  z  =  0. 

fjp.  2  =  0)  =  /,(/>.  z  =  0)  (115) 


and 


-47  [/,(/>.  2)],=0  =  -^-[/;(p.2)]r=0 


(116) 


Thus. 


Ag  —  B\\_Aik] (c/i.j)  ^':('7i,i)] 


(117) 


and 

jv.u  =  B,  WM  +  *>/).,)]( 

Dividing  Equation  US  by  liquation  117  gives 

(y  "jj"  '*i^,2(</i,|)  —  i .i ) 

si  j  =  ,  ”  -  - 

//^i(‘7i.i)  ~ 

(f  i  k  ;(</], i)  —  ] ' 

—  -  > 

—  q  i  ^’iO/i.i) 


(US) 


(U9) 


where 


In  the  program.  f(p,  z,)  is  set  to  1.  This  gives 


■  T  ^ 2 U/ 1 , i ) 


(120) 


(121) 
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For  i  =  2  to  I  - 1 .  A,  s  and  B,  's  arc  determined  by  matching  conditions  at  every  layer 
boundary.  At  layer  boundary  7.  =  r  , 


and 


77  [/(P.r.)]  =  ~  [y;_,(/>.r:)] 


From  Equation  12 


A  [-I;  ^\U]i,-)  +  /\'2f 

=  A—  1  [  -  ^  i— •  ^  1  ( <7;— ! ./  *  +  M‘7;-  Ij  0  ' 


Thus, 


A  « 


-  F'mF/;.;  .)  -r  A;(<7;  ■  j 


The  coefficient,  .1,.  is  determined  from  Equation  123  as  follows: 

#;[■'!/*' 1  (*/«)  +  )]</'/ 

=  /< ' , ( </, _],,■)  + 

Substituting  Equation  125  into  Equation  126  gives 

A  j  f  <7;-_ ,  )  +  k2{<Ji-\j) 

Ai  j  +  /.2((/;; }  1  (4u ) + a'2(*j  >y< 

=  L'*i-l,*’l(<7/-u)  +  k'lUll-] ,/Wi-l  • 


Thus. 


(123) 


(124) 


(125) 


(126) 


(127) 


■h 


"tj 

d-h 


(12S) 


where 


tl‘U  =  q'ik'iUhjWt- 1  ki(‘h-u)  + 

-  7'i-i  k'A<h-u )  +  )1 


(120) 


and 


dA-,  =  A,(<7u  )[-•!,_,  a',(<7i-u  )  +  *'2(‘7i-u  )] 

-  ?'/  a'i(<7;,<  )['0-i  ai(<7/-i,/  )  +  Mtt-u  )1  . 


(130) 


l  or  i  =  1  and  at  the  layer  boundary,  z  =  z,  . 

fi  0}*  */)  =  fi~\  O’-  r/ ) 


(131) 


This  gives 


fl/- ;[-■!/- 1  AlO'/ J-1,7  1  A;0//-1.;  )] 

/;2f'//V  ) 


(13: 


Tor  attenuation  rates  less  than  or  equal  to  0.1  dB  km.  .1,  s  and  /?, 's  arc  determined 
successively  starting  from  A,  and  fi,  down  to  A,  and  fi,  . 

The  height  gain  function  at  the  top  layer  is  given  by 

fl(p.Z,)  =  /(;(<//)  -  (133) 


where  the  coefficient,  fi,.  is  set  to  one.  Matching  the  boundary  condition  at  z  =  z,  gives 


_ M'Z/.O _ 

'0-1  Al0//— 1,/  )  +  ^ 


(134) 


and 
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/-i 
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(135) 
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where 


—  <■/  i  h  ?(Qij  )  ) 

—  <//_]  h-A<faj)k  2^ ?/— 1 ,/ ) 


d-h-\  —  <7'/-i  hiUh.i  )k'\{<li-\,i ) 

(137) 

—  <7*/  h'litfu  )^i(‘7/-i,/ ) 

Similarly,  the  coefficients.  .■/,  and  /?„  are  obtained  by  matching  boundary  conditions  at 
c  =  z,  for  i=  1-2  to  1. 

Thus. 


I*t  "  */4 


_£/+1  1  4-  1  ./4*  1  )  +  k,Uj,+  ]J+i  ) 

•V'i(‘/,\.+1  )  +  *i<‘7/./+l  ) 


where 


n-h  ~  Li  (+1  1  1C-1/+1  ^  l((7;+l,<+)  )  "b  A"  2(^4  1 ./+]  )1 

—  '/  (  A  l(l1i.i+ 1  )[•  1/+I  A'l(<7;-H,/+1  )  +  ^:(<7(+l,;+l  1] 


<7/1;  —  <7'/  A  i(<7i(+i  )[^i+i  ^i(<7(+i,(+i  )  +  ^2(c7(+i,/+i  )] 
—  tfi+l  ^l(</j,i+l  K'^+l  A  t(9i+l,i+l  )  "b  ^(tfi+U+l  )] 


D.  TEST  FOR  EVANESCENCE 

The  functions.  /<,(</)  and  k2(q).  in  Equation  26  and  27  can  also  be  expressed  in  terms 
of  the  modified  Ilankcl  functions  of  order  one  third  as  follows  (sec  Appendix  A): 


*,(</)  =  -j2U2)'!*Ai{-qej:!:i-) 


_  I  2.  n  «>  /  2  3/:\ 

-\T?  j  ">»  (-«  j 


=  *1  ('/) 


A-2(4/)  =  2(1 2)'i6eM(,Ai(-q) 


=  h2{q)  ~  cj4~^  h^q)  , 


where 


*,W  -  (f./5;i),f3H,)3,"(t"3;!)  (144) 

and 

h2(q)  =  j2U2)'!t  Ai{~qc^) 

-(  4v'  -)  ’//  .  =  (  fv-)  . 

For  large  |  <7 1,  Ar,(</)  and  A,(</)  can  be  approximated  by  the  asymptotic  expansion  of  the 
Ilankel  function  [Ref.  4  and  Appendix  A].  That  is, 


,,->,1/6  -1/2  -1/4  (  2  3/2  arr 

(12)  71  q  exp  ./(  y  q - jy 


J-  <  arg(<7)  <  ~ 


,,->.1/6  -1/2  -1/4  .(2  3U  5tt  \ 

(12)  ■  n  1  q  exp  -  <r  -  7T  j  : 


0  <  arg  (<7)  <  2 n  . 

From  Equations  146  and  147,  it  can  be  deduced  that  for  large  |  q  |  and  in  the  region 


2S 


(148) 


4- 

-=r-  <  are(</)  <  -7- , 

>  w  3 

k}U])  will  be  exponentially  large  and  A_,(</)  exponentially  small.  When  this  occurs, 
/i,(t/)  and /;,(</)  become  linearly  dependent  numerically,  i.e..  from  liquation  143. 

/hiq)  ^  7)  .  (149) 

This  implies  that  A,  (q)  and  k2  (</)  are  also  linearly  dependent  numerically.  Under  this 
condition,  the  matrix  of  the  modal  function  will  become  singular  (in  a  numerical  sense) 
for  all  values  ofq.  Pappert  and  Goodhart  defined  the  fields  as  evanescent  when  liquation 
149  holds  [Ref.  5  and  6). 

When  evanescence  occurs,  the  height  gain  function  is  represented  by  an  exponen¬ 
tially  decaying  field,  i.e.. 


f(P'=)  =  *:(</;)  • 


( 1 50) 


It  is  necessary  to  incorporate  test  conditions  in  the  program  to  check  if  the  rf  field 
has  become  evanescent.  This  is  achieved  by  tracking  the  real  part  of  the  exponent  of 
in  Equation  146. 

Since  evanescence  may  only  occur  in  the  region  specified  by  Equation  148,  it  can 
be  deduced  by  simple  transformation  of  Equation  14S  that 


< 


(151) 


This  shows  that  one  needs  only  to  test  for  evanscence  when  the  real  part  of  the  exponent 
of  A,  (q)  is  positive.  That  is. 


Re 


>  0 


(152) 


As  the  coefficients  of  the  height  gain  function  are  dependent  on  how  the  integral  of 
Equation  SS  is  evaluated,  it  follows  that  different  test  conditions  exist  for  upward  and 
downward  integration.  Eor  upward  integration,  conditions  for  non-cvanscence  in  the  i'h 
layer  arc 


IN 


or 


Rr  £  y ./'  (‘7/.1+I  )3'2]  -  *{  y./'  {*li/  2]  <  7 

for  /?<’£  yy  (</,-..>+, )3,‘J  ^  0 

and  ^<-|^yJW;.,)3/2J  ^  0 


for  /?«:[" -i-y^)3''2 
and  Rrf 


<  0 

Uo 


(154) 


or 


>  0 


( 1 55) 


or 


r(+1  <  45  m  .  (156) 

In  principle,  the  right  side  of  liquation  153  should  be  as  large  as  possible.  However,  a 

large  quantity  will  lead  to  severe  cancellation  error  when  extracting  k2  (<7)  (which  is  a 

small  number)  from  a  linear  combination  of  two  large  numbers.  h}  (<])  and  h:  (</)  . 

*> 

Equation  154  is  obtained  from  Equation  153  by  setting  Rc[  y./(‘7..)3  ■’]  to  zero  when 
its  value  is  negative.  This  is  true  because  when  the  real  part  of  the  exponent  of  A,( <7)  is 
negative,  the  condition  of  evanscence  will  not  occur. 

J)]. 

With  regards  the  condition  in  Equation  155.  it  can  be  shown  that  when  y—  is  small 
(which  is  usually  the  case)  and  Equation  155  is  true,  then 

Re  (a,-)  >  0  .  (157) 

This  implies 

»i2(r/+1)  >  ni~  (r;)  ,  (15S) 


which  leads  to 
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Ri'L(h.i+ 1]  >  AY[</;,] 


(159) 


or 

/^[t»/^,+i)3;2]  <  7^’  [i'-/,V/;')3:]  (l6°) 

or 

j^  —  /  (?iV+,  )';2J  -  Re  £ -jj  (q:^a  <  0  .  (161) 

By  comparing  liquation  161  with  Equation  153.  it  is  clear  that  Equation  161  is  a  condi¬ 
tion  for  non-evanescence. 

Finally,  the  condition  in  Equation  156  was  based  on  the  fact  that  at  lower  altitude, 
the  magnitude  ofq  is  small  enough  such  that  hx(q)  and  h2  (q)  arc  well  behaved  in  the  re¬ 
gion  stated  in  Equation  14S. 

By  the  same  argument,  conditions  for  non-evanescence  in  the  /"■  layer  for  downward 
integration  are  as  follows: 

Rc  i.'— ,)  :  ~  <  7 

for  (<7/-, ./-i)'  2]  ^  0  (162) 

and  ^  0 


or 


for  /?*7[‘3*^(?f-u)3'2J  <  0 
and  ^  0 


(163) 


or 


(164) 
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or 


~j  <  45/??  .  (165) 

E.  EVALUATION  OF  THE  MODAL  FUNCTION 

When  evaluating  the  relative  field  strength  of  the  electromagnetic  waves  using 
Equations  100  and  101.  the  eigenvalues,  p  =  pn  ,  Tor  which 

II  A  ||  =  0  (166) 

are  required.  The  determinant,  ||  A  ||,  is  known  as  the  modal  function.  Each  eigenvalue 
represents  a  different  electromagnetic  mode  that  propagates  in  the  multilayer  tropos¬ 
pheric  waveguide  environment. 

The  roots  of  Equation  166  are  found  by  the  Shellman-Morlitt  complex  root-finding 
routine  to  be  discussed  in  Chapter  III.  Instead  of  searching  for  solutions  of  Equation 
166  in  the  p-spaee,  it  is  more  convenient  to  conduct  the  search  in  the  </u -plane.  The 
variable,  p.  is  related  to  qux  by 


In  searching  for  the  roots  of  the  modal  function,  the  determinant  must  be  evaluated 
many  times.  Therefore  an  efficient  method  for  accomplishing  this  is  required.  In  partic¬ 
ular.  the  Laplace's  expansion  is  used  to  obtain  the  determinant.  The  Laplace's  expansion 
is  based  on  the  following  theorem  [Ref.  7). 

Theorem:  Laplace's  Theorem  states  that  if  one  selects  any  r  rows  of  a  de¬ 
terminant.  |  .1  |  ,  forms  all  possible  r-rowed  minors  from  these  r  rows,  multi¬ 
plies  each  of  these  minors  by  its  co-factor  and  then  adds  the  results,  one 
obtains  |  A  |. 

The  method  also  takes  advantage  of  the  presence  of  zero  elements  in  the  determi¬ 
nant  to  reduce  numerical  computations. 

The  understanding  of  this  technique  can  be  best  achieved  with  an  example.  Con¬ 
sider  a  trilinear  refractivity  profile  with  the  following  modal  function: 
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Let  the  last  two  rows  be  the  r  rows  from  which  minors  are  to  be  formed.  Then  applying 
Laplace's  Theorem,  the  modal  function  is  given  by 
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,  (169) 


If  a  four-layer  piecewise  linear  refractivity  profile  is  used,  the  modal  function  will  be 
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Notice  that  Equation  170  is  built  from  Equation  16S  by  adding  eight  non-zero  entries, 
namely.  c,Sf..  a,<,  %.  a.,-.  a-<.  a-,  and  a-  .  Therefore,  tire  results  from  Equation  108 
can  be  used  to  compute  Equation  170,  i.e.. 


II  A||4  =  ||A!|,  . 
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In  Equation  171.  the  determinant  of  the  5x5  matrix  is  also  evaluated  using  the 
Laplace's  expansion. 

Without  loss  of  generality,  the  evaluation  of  modal  functions  of  multilayer  piece¬ 
wise  linear  refractivity  profile  is  performed  by  systematically  increasing  the  number  of 
rows  and  columns  of  A  and  applying  the  Laplace's  expansion. 

The  derivative  of  the  modal  function.  ||  A  ||.  with  respect  to  is  also  required  by  the 
Shcllman-Morlitt  root  finding  routine.  This  is  accomplished  by  expanding  ||  A  ||  using  the 
Laplace's  expansion  and  performing  differentiation  using  the  chain-rule. 


F.  EFFECTS  OF  SURFACE  ROUGHNESS 


The  materials  in  this  section  arc  adapted  from  Chapter  III  ofXWVG  [Ref.  Ij.  It  is 
included  in  this  document  for  completeness. 

Surface  roughness  could  have  significant  effect  on  the  signal  levels  at  large  ranges, 
especially  when  the  frequencies  of  the  rf  waves  are  above  several  gigahertz. 

The  effects  of  surface  roughness  is  modelled  by  Baumgartner  in  XWYG  (Ref.  1]  as 
a  lossy  infinitesimally  thin  bottom  layer  of  constant  modified  rcfractivity  shown  in  Fig¬ 
ure  4. 

In  the  bottom  layer.  0  <  z  <  zr  the  height  gain  function,^  (p,  z),  can  be  written 
as 

/u  (>.:)  =  A0  +  Rc~Jn:  ]  .  (172| 

where 

H  =  (0)  -  ft2  ,  (173) 

/>2  =  A-  (174) 

A* 

and  R  is  the  plane  wave  reflection  coefficient. 

By  matching  boundary  conditions  at  each  layer  boundary  as  was  done  in  Section 
A.  it  can  be  seen  that  all  previous  results  still  hold  if  the  following  substitution  is  made: 

v  — +  v  =  n  — - f  1  /  *> t 

i  is  .1  n  t  *  '  -  / 


In  general,  the  bump  heights  of  the  rough  surface  are  assumed  to  be  Gaussian  dis¬ 
tributed.  In  this  case,  the  refection  coefficient,  R.  can  be  represented  as  the  product  of 
the  smooth  surface  plane  wave  reflection  coefficient  and  a  surface  roughness  factor  (Ref. 
8],  For  horizontal  polarization,  the  reflection  coefficient.  R,  is  given  by 


Rc  (<?,  ,)  >  0 
Re(<7u)  <0 


(176) 


where  6  is  the  root  mean  square  surface  bump  heights  and  RH  is  the  Fresnel  reflection 
coefficient  for  a  horizontally  polarized  wave  given  by 
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where  R  is  the  Fresnel  relleetion  coefficient  for  a  vertically  polarized  wave  given  by 
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For  horizontal  polarization.  can  be  expressed  as  f  ollows  [Ref.  1]: 
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(is: 
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and 
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For  small  values  of  ( <7,  ,  | ,  «(<7U)  and  b(qlA)  can  be  expanded  in  series  as 
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The  derivative  of  ;,w  with  respect  to  qu  can  be  obtained  by  dilrerentiatinc  liquation 
ISO.  i.e  . 


C''M!  k  , ,  .  ra  ,  ,  rk 

- -  =  — -  %i  i)  — - -  fl(<7i  ,)“ - 

((h.\  />-(</,,)  L  < '/u 


c  *  h 

Similarlv.  — - and  — —  are  obtained  from  liquation  1S1  and  1S2  respectively; 

(</u  r  (/ ii 


3S 


and 


<J1  _  T];:  •,  ai  \  1/3  .  3  7i 

T  t7:-'  I  2  A  ‘  V  A  j  4T  A  '  V  A 


«7i.i 


tA7  +  — A‘ii- 
J  8T- 


T^rA-V  +  — -  Tk2S2 

(>  I  1  /-  T'J 


16T 


(t)". 


+  t 

+4- 


-VA'W-^- 

24T~ 


,  6  5  ,2  -2 

A  d  —  ~  A  () 

12ST 


15T 


twr — 4tv+— tv 
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256T5 


(191) 


'1  \3  4 

<7l,1 


i  \  n  .•  i  5 


T"  '  ‘(m  '  +  -• 


The  corresponding  results  for  vertical  polarization  are: 


>.s  r 


^  f  / : .  ] ) 
£(‘7i.i) 


(192) 


and 


,  sr 


A 


C<> 


1.1  /’*(</,,,) 


/’(^/i.i)  - ‘'Wi.i)  r/> 


<'</!, 1 


(193) 


where 


s («,.,)  -  (t)1"  [t+(-^-)w <>»•»> 


and 


*<*.,)  -  tanh(4)[T  +  (-^- (195 


5) 


Tor  small  values  of  |  </, ,  )  .  f/(</; ,)  and  />(</,,)  becomes 
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",  X  ,  J-\  -.1/3  1/2  j  ,  2r2,  J\  x2/3  1  ,6 -A.  *1  '  3 

)  <h.)  |  A*M~ <7i.i  --T  *  «M~  )  <7i.i 

2  ,  lOclO,  -.10/3  5  1 

+  Uk  *  (T]  ^  " . / 


,  »r«->)  -1/2  f.  .  1  »l  ,2/3 
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1  /  _?4_  .4/3  2  ]  ,  Jl  ,2  3 

8T2  (  /•  )  ?U  +  „T1  (  *  )  ?U 


16T3  * 


5  ,  al  vS/3  4  .  7  ,  al  ,10/3  5 

NT*  '  k  ’  *■'  +  ^Z~F  l  — )  9u 


M«u>  =  1  +  “ 4“  T1'2  </,  i  '''  {/;2<t2 


2  -2  ,  a 1 


P(t)  +lTr,nt^ 


-  I  —  4-  — ! —  /,vr 


3A  l  +T^A‘5‘  ("f}  f/l- 


+  |_~  6T  kS  +7^r"J(“)  *>• 


2  ,.10  =10  ,  1  ,  "  ;fi  >  ,  2  c:  ,  *1  .3  4 

+  TT  *  15  ^ - r  A  0 - r  k  <>  (  —  (/,  , 

13  24T"  12ST  J  A  • 


i  10  ;  10  1  ,/'c6  .  7  ,2-2  ,  y!  .11/3  5  , 
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The  derivative  of  «(</,  ,)  and  A>(c/,.t )  are  obtained  from  liquation  196  and  147  respec- 
tivclv: 


canlHy) 
+  *Vi-?1.,'»s«h!(4) 
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and 


,?Li2l  T1/2  7  — 1/2  f  1  ,2  ,2  al  .1/3 
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— >7l a-y|^y  A  J  (  “  ' 
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G.  COMPLEX  INDEX  OF  REFRACTION  Or  SEA  WATER 

The  materials  in  this  section  are  adapted  from  Chapter  IV  of  XWYG  (Ref.  1], 

The  complex  index  of  refraction  is  one  of  the  physical  parameters  required  for 
evaluating  the  modal  function  and  height  gain  functions.  The  complex  index  of  refrac¬ 
tion  of  sea  water,  n..  is  a  function  of  frequency,  temperature  and  salinity.  It  is  related  to 
the  effective  relative  dielectric  function.  e,„-.  and  the  effective  conductivity. 
a,„ .  through 


rlS  ce/r  j  ( •>£ 


where  c9  is  the  permittivity  of  free  space  and  a>  is  the  frequency  at  which  n.  is  evaluated. 
The  effective  relative  dielectric  function  and  the  effective  conductivity  are  related  to  the 
complex  dielectric  function,  e,  and  the  ionic  conductivity,  o.  through 


e  =  £  -  ji  , 


tpff  —  £ 


ocfr  =  a  +  u)  £0  £' 


42 


In  liquation  205,  a  is  assumed  to  be  real.  This  is  a  good  approximation  for  frequencies 
much  less  than  10-  GHz  [Ref.  lj. 

The  theory  of  an  ideal  polar  dielectric  in  an  alternating  electromagnetic  field  was 
first  formulated  by  Debye  [Ref.  9].  In  the  theory,  the  complex  dielectric  function  is  ex¬ 
pressed  as  a  function  of  a  single  relaxation  time.  t.  as  follows: 

£  =  £'  -jc"  ,  (206) 


* 

1  +  (0~  t" 

and 


I  +  t  >:  z; 


where  is  the  part  of  the  dielectric  function  due  to  the  atomic  and  dielectric  polariza¬ 
tion.  cs  is  the  static  dielectric  function  and  u)  is  the  frequency  of  the  electromagnetic 
wave.  Equations  206  to  20S  represent  the  fall  in  the  value  of  the  dielectric  function  from 
c,  to  This  fall  is  accompanied  by  a  single  broad  absorption  band  in  the  vicinity  of 
the  characteristic  wavelength.  given  by 


(209) 


The  characterisation  of  the  variation  with  frequency  of  the  complex  dielectric  function 
of  water  in  terms  of  a  single  relaxation  time  is  valid  for  frequencies  up  to  300  GHz 
[Ref.  1).  Tor  higher  frequencies,  more  than  one  relaxation  time  may  be  required. 

In  the  model  of  Klein  and  Swift  [Ref.  lOj.  is  assumed  to  be  4.9  for  both  pure 
water  and  seawater.  With  this  assumption,  the  static  dielectric  function  for  pure  water 
is  given  by 

r.S  P  =  88.045  -  0.4147/,  -  6.295  x  10~4  /,2  +  1.075  x  10“5  /c  3  (210) 


and  the  static  dielectric  function  of  seawater  is  given  by 

£5,5  =  £5.<z  a  (■S*  l()  >  (211) 
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where 


cSa  =  S7.134-  1.949  x  l')"1  tc  —  1.276  x  10'*//  +  2.461  x  1<)~  /.  3  (212) 


and 


a  (S.  ic )  =  1.000  +  1.613  x  10~5  .S  tc  -  3.656  x  10~3  S 


+  3.210  x  lO"5  S2  -  4.232  x  10""  5 3 


(213) 


In  the  above  equations,  i,  is  the  temperature  in  degrees  Celcius  and  S  is  the  salinity  in 
grams  of'  salt  per  kilogram  of  seawater.  The  relaxation  time  of  seawater  is  given  by 


~s  =  ~p  h  (5,  K  )  • 


!14) 


where  r.  is  the  relaxation  of  pure  water  gi\  cn  by 

zr  =  1.76S  x  l<r"-6.0S6x  1 0" 13  rc  +  1.104  x  10'1J  ic  2  -  S.l  1 1  x  10_17/e3  (215) 


and 

b{S,  i.)  =  1.000  +  2.282  x  10"5  .S  tc  -  7.63S  x  10"4  S 
-  7.760  x  lif"  S2  +  1.105  x  10-8  53  . 

The  ionic  conductivity  of  seawater  in  siemens  per  meter  is  given  by 

a(tc.  S)  =  ct(25.  5)  exp(  —  Ay)  . 

where 

A  =  25  -  /,  . 


(216) 


(217) 


(21S) 


and 


<?  =  2.033  x  10-2  +  1.266  x  10"4  A  +  2.464  x  10~6  A2 
-  S[  1.849  x  10"s  -  2.551  x  10"7  A  +  2.551  x  10"8  A:] 


a(25,  S)  =  S  [0. IS2521  -  I.46l92x  10"3  .S  +  2.09324 
x  10" 5  S2  -  1.28205  x  10"'  S3  ]  . 


(219) 


(220) 
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H.  ATMOSPHERIC  ABSORPTION 


The  materials  in  this  section  arc  adapted  from  References  11.12  and  13. 

The  gas  molecules  in  the  air  absorb  energy  from  electromagnetic  waves.  The  ab¬ 
sorption  is  due  to  the  quantum  mechanical  resonances  of  the  gas  molecules.  These  res¬ 
onances  at  specific  frequencies  then  broaden  into  bands  as  a  result  of  collisons  of  the 
gas  particles.  In  the  lower  atmosphere,  the  density  of  gas  particles  is  higher.  This  means 
that  collisons  are  more  frequent  and  the  absorption  bands  arc  broader. 

In  the  troposphere  (height  <  10  km),  only  the  lines  bands  of  oxygen  and  water 
vapour  are  important.  The  principal  resonance  of  oxygen  is  at  60  GHz  and  that  of  water 
vapour  is  at  22.2  GIlz.  The  extent  of  absorption  by  oxygen  molecules  is  dependent  on 
the  pressure,  p.  and  temperature,  t.  In  the  case  of  water  vapour,  absorption  is  also  de¬ 
pendent  on  the  water  vapour  density.  //.  Tor  specific  values  of  p.  t,  and  />.  the  absorption 
coclUcicnt.  gci  .  is  computed  as  [Ref.  1 1] 


£m.>,  =  O.IS20/.\>  (Wjkm  ,  '221) 

where  f  is  the  frequency  and  .V,  is  the  absorption  spectrum.  The  absorption  spectrum. 
.V, .  is  related  to  the  line  spectra  of  absorption.  SI',  the  continuum  spectra.  .Yr .  and  the 
liquid  water  extinction.  ,Y„  .  by  the  following  expressions  [Ref.  1 1J: 

A'.  =  VsF  +  \  +  .Y.  ppm  .  (222) 

i 

The  summation  of  SI  is  done  over  44  spectral  lines  of  oxygen  and  29  spectral  lines  of 
water  \apour,  [Ref.  II].  i.c.. 

Ysr  =  (fl^3  exp  to  (Od  -  0)]  r 

i  i=  i 

(223) 

+  2^h i  (/')  c  O'  5  exp  [b2  (i)  (1  -  0)]  F  , 

i=  i 

where  «,(/),  a2  (/).  />,(/)  and  b2  (i)  arc  spectroscopic  coefficients  derived  from  experiments 
and  spectroscopic  parameter  compilations.  Values  of  a, (7).  a2  (/).  (?,(/)  and  b2  (/)  arc  listed 
in  Table  1. 


The  dry  air  pressure,  p.  at  a  height  of  z  above  the  earth  s  surface  is  given  by 
[Ref.  12] 

f  ,  0.006.V  1 5.256  ,  n 

r  -  1111  - —  J  kra  ■ 

The  relative  inverse  temperature,  0 .  is  related  to  the  temperature.  tk .  in  °K  by 

too 


(224) 


e  = 


(225) 


The  water  vapour  partial  pressure.  <?„.  in  kPa  is  related  to  the  saturation  water  vapour 
pressure.  es ,  and  relative  humidity,  rh.  by  the  following  equations: 


eM.  =  es  {rh)  x  10  “  , 


es  =  0.6105  exp 


(226) 


kPa  (227) 


and 


tc  —  tk  —  273.15  °C  . 


(228) 


In  Equation  223,  F  is  the  shape  factor  of  the  broadened  resonance  lines  given  by 
[Ref.  11] 


(f_\ 

.e  -  (vn  ~f)c 

l 

+' 

1 

\  r°  / 

_  (' o -fr  +  s2 

(V0  +./)2  +  S2  J 

:  3S  and 

a 

•> 

<T 

[  vo  ) 

.  (‘o -J) 2  +  S2 

j 

(vo  +JY  +  s'  - 

(229) 


(230) 


for  oxygen  line  39  to  44  and  the  water  vapour  lines. 

In  Equation  229  and  230,  v0  is  the  molecular  line  centre  frequency,  g  is  the  width 
and  c  is  the  overlap  interference. 
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For  oxvsen  in  air. 


g  =  a,  (/>  +  1.3  c)  ()'  'v  . 
c  =  ch  p  6t,s 


and  lor  water  vapour  in  air, 


g  =  b3  (4,8  e  +  />)  0 


,0.6 


The  line  data  a}.  a,  and  /\,  are  listed  in  Table  1. 
The  continuum  spectrum.  .Vf .  is  given  by 

.Vc  —  1.9/rpO'-'  x  lo ~r  +  6.25 fpO2  x  10" 


(231) 

(232) 

(233) 


l  r  +  gu 


-  +  2.1  pO  -  x  10" 


/>/v?i(234| 


and 


g0  =  0.01 2(/»  +  1.3  e)0°'9  .  (235) 

The  liquid  water  extinction.  .Y„  ,  is  evaluated  from  the  following  expressions: 

A;,.  =  - - -  ppm  for  /  <  300  GHz  (236) 

(e'  +  2)*  +  (£")* 

and 

Au.  0.55  u/-0’1  $~6  ppm  for  /  >  300  GHz  ,  (237) 

where  w  is  the  liquid  water  concentration  in  grams  per  cubic  meter.  The  real  and  im¬ 
aginary  part  of  the  dielectric  function  of  warer.  z’  and  c",  are  given  by  Equations  20?  and 
20S  respectively.  It  is  assumed  that  z.  and  t  given  in  Equations  207  and  20S  arc  given 
by  Equations  210  and  215  respectively. 

The  absorptivity  of  tropospheric  gases,  tj,  is  related  to  the  atmospheric  absorption, 
gcb ,  in  dB  km.  by  the  following  expression  (Ref.  13): 

Sab  =  -y"  10*  i]  Iog10  e  .  (23S) 
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Table  1.  SPECTROSCOPIC  COEITICIENTS  FOR  OXYGEN  («,-«<)  AND 
WATER  (/;,  -  />,)  [FROM  RET.  1 1] 


Center 


frequency,  ►„ 

Strength,  n, 

1  enipernture 

Width, 

Interference,  o. 

Temperature 

Quantum  number 

(Cllz) 

(kllz/kfa) 

exponent,  a , 

(Gllx/kl’a) 

(l/krn) 

exponent,  o5 

identification,  K  * 

0" 

0.307  E-3  (ppm/kTa) 

0.12  (0.056)  E- 

1 

i* 

to  37* 

30.47381 

0.94  E-6 

0  969  El  1 

0.86  E-2 

0  520  E-2 

0.179  E+l 

37- 

50.98742 

0.244  E-5 

0,869  E+l 

0  87  E-2 

0.550  E-2 

0.169  E+l 

35' 

51  50311 

0.604  E-5 

0.774  E+  1 

0.89  E-2 

0.560  E-2 

0.177  E+  1 

33' 

52.02124 

0.141  E— 4 

0.684  E+l 

0.92  E-2 

0.550  E-2 

0.181  E+l 

3 1 " 

52.54227 

0  308  E-4 

0.600  E  +  1 

0.94  E-2 

0.569  E-2 

0  179  E  4  1 

29- 

53.06683 

0.637  E-4 

0.522  E+l 

0.97  E-2 

0.528  E-2 

0.189  E  +  l 

27- 

53.59570 

0.124  E-3 

0.448  E+l 

0.100  E-l 

0.544  E-2 

0.183  E  +  l 

25“ 

54.12997 

0.2265  E-3 

0  381  E  l  1 

0.102  E-l 

0.480  E-2 

0  199  E+  1 

23" 

54.67115 

0.3893  E-3 

0  319  EH 

0.105  E-l 

0.484  E-2 

0  190  E  l  1 

2 1 " 

55.22137 

0.6274  E-3 

0.262  E+  1 

0.1079  E-l 

0  417  E-2 

0.207  El  1 

19" 

55.78382 

0.9471  E-3 

0.212  E+l 

o.i  no  E-i 

0.375  E-2 

0.207  E  l  l 

17- 

56.26477 

0.5453  E-3 

0.100  E-t 

0.1646  E-l 

0.774  E-2 

0.890 

Dl* 

1* 

56.36339 

0.1335  E-2 

0.166  E+l 

0.1144  E-l 

0.297  E-2 

0.229  El  1 

15- 

56.96818 

0.1752  E-2 

0.126  E+l 

0  1181  E-l 

0  212  E-2 

0  253  E+l 

13 

57.61249 

0.2125  E-2 

0.910 

0.1221  E-l 

0.940  E-3 

0.376  E+  1 

i  r 

58.32389 

0.2369  E-2 

0.621 

0.1266  E-l 

-0  550  E-3 

-0  III  E  l  2 

D2 

r 

58.44658 

0.1447  E-2 

0.827  E-l 

0.1449  E-i 

0.597  E-2 

0.790 

3* 

59.16422 

0.2387  E-2 

0.386 

0.1319  E-l 

-0.244  E-2 

0.700  E-l 

V 

59.59098 

0.2097  E-2 

0.207 

0.1360  E-l 

0.344  E-2 

0.490 

5* 

60.30004 

0.2109  E-2 

0.207 

0.1382  E-l 

-0.435  E-2* 

0.680 

D3 

5' 

60.43478 

0.2444  E-2 

0.386 

0.1297  E-l 

0.132  E-2 

-0.120  E+l 

7* 

61.15057 

0.2486  E-2 

0.621 

0.1248  E-l 

-0.360  E-3 

0.584  E+l 

9* 

61.80016 

0.2281  E-2 

0.910 

0.1207  E-l 

-0.159  E-2 

0.286  E+l 

11* 

67  41122 

'  0.1919  E-2 

0. 126  E  +  l 

0.1171  E-l 

-0.266  E-2 

0.226  E+l 

D4 

13* 

62.48626 

0.1507  U-2 

0.827  E-l 

0.1468  E-l 

-0.503  E-2r 

0.850 

3‘ 

62.99797 

0  1492  E-2 

0 .166  E  +  l 

0.1139  E-l 

-0.334  E-2 

0  218  EH 

15* 

63.56852 

0.1079  E-2 

0  212  E+l 

0. 1 108  E-  1 

-0.417  E-2 

0  196  E+l 

17* 

64.12778 

0.7281  E-3 

0.262  E+l 

0.1078  E-l 

-0.448  E-2 

0.200  E+l 

19* 

64.67886 

0.4601  E-3 

0.319  E+l 

0.105  E-l 

-0.515  E-2 

0  184  E+l 

21* 

65.22412 

0  2727  E-3 

0.381  E+l 

0  102  E-l 

-0.507  E-2 

0.192  E  +  l 

23* 

65.76474 

0  152  E-3 

0.448  E+l 

0.100  D-l 

-0.567  E-2 

0  178  E  +  l 

25* 

66.30195 

0.794  E-4 

0.522  E  +  l 

0.97  E-2 

-0  549  E-2 

0  184  E+l 

27* 

66.83663 

0.391  E-4 

0.600  E+l 

0.94  E-2 

-0.588  E-2 

0.174  E+l 

29* 

67.36933 

0.181  E-4 

0.684  E4  1 

0.92  E-2 

-0.560  E-2 

0  177  E+  1 

31* 

67.90051 

0.795  E-5 

0.774  E  l  l 

0.89  E-2 

-0.580  E-2 

0.173  E  l  1 

33* 

68.43054 

0.328  E-5 

0.869  E+l 

0.87  E-2 

-0.570  E-2 

0.165  El  l 

35* 

68.95972 

0.128  E-5 

0.969  E+l 

0.86  E-2 

-0.530  E-2 

0.174  E+l 

37* 

118.75034 

0.9341  E-3 

0.000 

0.1592  E-l 

-0.441  E-3 

0.890 

r 

49 


Table  I.  (continued) 


Ccnler 

frequency,  »-0 
(Gllr) 

Strength,  a, 
(kllr/kl'a) 

1  enipetaluie 
exponent,  <?, 

Width,  a, 
(Gllz/kl'a) 

Quantum  number 
identification  (O, ) 

Lower 

Upper 

368-198)50 

0.679  E-4 

0  200  E-l 

0.156  E-l 

1.  1 

2.  3 

424.76)120 

0  6)8  U-3 

0  122  E-l 

0  147  E-l 

2.  1 

2,  3 

487.249371 

0.235  E-3 

0.122  E-l 

0.147  E-l 

2,  1 

3.  3 

715.393150 

0.996  E-4 

0  891  E-l 

0  144  E-l 

3,  3 

4.  5 

773.8)9732 

0.571  11-3 

0  798  E-l 

0.140  E-l 

4,  3 

4.  5 

834.145790 

0.180  E-3 

0.798  E-l 

0.140  E-l 

4,  3 

5.  5 

Identification  (11,0) 

*'o 

*, 

>>, 

<>, 

Lower 

Upper 

•  22.235080 

0.105 

0  214  U  1  1 

0  281  E-l 

5.  2,  3 

6.  1.  6 

68.052 

0  180  fl-2 

0  8  75  I:  I  1 

0  280  E-l 

3,2.  1(1)' 

4.  1.  3 

•  18)310091 

0  238  E  l  1 

0  653 

0  282  E-l 

2.  7,  0 

3.  1.  3 

321.225644 

0.460  E-l 

0.616  E-l  1 

0  220  E-l 

9.  3,  6 

10.  2.  9 

•  325152919 

0  155  Ell 

0.152  E  1  1 

0  290  E-l 

4,  2,  2 

5.  1.  5 

•  380.197)72 

0.123  E  +  2 

0  102  E  l  1 

0  285  E-l 

3,  2.  1 

4.  1.  4 

386.778 

0.400  E-2 

0.733  E  1  1 

0  160  E-l 

II.  2.  10 

10.  3,  7 

4)7.34667 

0  630  E-l 

0  502  E  l  1 

0  150  E-l 

6,  6.  0 

7.  5.  3 

439.150812 

0.921 

0.356  E  i  1 

0  175  E-l 

5.  5.  0 

6.  4.  3 

44)  018295 

0  191 

0.502  1:4  1 

0.148  E-l 

6,  6,  1 

7.  5.  2 

■  448  001075 

0  107  E  l  2 

0.137  El  1 

0.246  E-l 

3.  3,  0 

4.  2.  3 

470.888917 

0.328 

0.357  E  l  1 

0.181  E-l 

5.  5.  1 

6,  4,  2  • 

474.689127 

0.124  E4  1 

0  2)4  EH 

0  210  E-l 

4,  4.  0 

5.  3.  3 

488  49113) 

0  256 

0.281  El  1 

0.222  E-l 

7.  1.  7 

6.  2.  4 

504  219 

0.380  E-l 

0.669  E  l  1 

0  127  E-l 

7.  7.  0 

8.  6.  3 

505  126 

0.120  E-l 

0.669  E  1  1 

0  130  E-l 

7.  7.  1 

8.  6.  2 

•  556.936002 

0  526  E4  3 

0  114 

0  317  E-l 

1.  0.  1 

1.  1.  0 

•  620.700807 

0  521  Efl 

0  2)4  El  1 

0  216  E-l 

4,  4.  1 

5.  3.  2 

658.340 

0.460 

0  776  E  l  1 

0  328  E-l 

1.0.  1(1) 

1.  1.  0 

•  752.0)3227 

0  259  E  l  3 

0.336 

0  302  E-l 

2.0.  2 

2.  1.  1 

836.836 

0  120  E-l 

0  811  E  l  1 

0  170  E-l 

11.2,  9 

10.  5.  6 

859  810 

0  150  E-l 

0,799  El  1 

0  270  E-l 

2.0.  2(1) 

2.  1.  1 

899.380 

0  910  E-l 

0.781  E  l  1 

0.300  E-l 

1.  1.  1  (1) 

2.  0.  2 

903.280 

0  640  E-l 

0  835  I:  l  1 

0  280  E-l 

2.2.  1(1) 

3.  1.  2 

907.773 

0  1 79 

0  504  II  i  1 

0  201  E-l 

8.  3.  5 

9.  2.  8 

•  916.169 

0.890  E  l  1 

0  137  II  l  1 

0  249  E-l 

3.  3.  1 

4.  2.  2 

970  320 

0  940  E4  1 

0  184  El  1 

0  246  E-l 

4.  3.  1 

5.  2.  4 

•  987.910 

0  145  I;  l  3 

0  180 

0  299  E-l 

1.  1.  1 

2.  0.  2 

•  1097.368 

0.840  El  3 

0  656 

0.335  E-i 

3.  0,  3 

3.  1.  2 

Read  0.307  ti  — 3  as  0  307  x  I0'\ 

'Nonresonant  O,  spectrum  (equation  (21)). 

*131  denotes  doublet 

'  Ro.itnkrani's  |I975|  first-order  solution  for  tire  60  C5 1 1  x  band  shape  hos  a  shortcoming:  the  listed  values  of  n,(K  «•  3", 
5')  have  to  be  teduced  5%  (which  is  of  negligible  consequence)  to  assure  that  «(0,)  &  0.000  for  /  >  160  Gilt  even  when 
m  0  (equation  (21). 

'Here  (I)  denotes  first  vibrationally  excited  stale. 
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III.  SHELLMAN  -  MORFITT  ROOT  FINDING  ROUTINE 


A.  GENERAL  OUTLINE 

'I  he  Shellman-Morfitt  root  finding  routine  was  developed  by  C.I1.  Shellman  and 
D.G.  Morfitt  to  locate  ELF  VI. I  LF  mode  constants  of  wave  propagating  in  an  earth- 
ionosphere  waveguide  environment.  The  theory  of  the  algorithm  described  in  this 
chapter  is  based  on  References  1  and  14. 

1  he  complex  root  finding  routine  is  used  to  find  the  zeros  of  the  modal  function. 
The  routine  will  locate  all  simple  zeros  of  an  analytic  function  in  a  prescibed  rectangular 
region  of  the  complex  z-plane.  The  principle  of  the  root  finding  method  is  based  on  the 
following  theorem  |Rcf.  15): 

Theorem  (Argument  Principle): 

Let  I)  be  a  simply  connected  domain  and  f(z)  be  anal}  tic  in  D  except  at  a 
finite  number  of  poles.  Let  F  be  a  closed  contour  in  I)  not  passing  through 
any  of  the  zeros  or  poles  cr  fvz ).  Then,  the  accumulated  phase  change.  A  o. 
of  [(/)  around  T  traversed  in  a  clockwise  direction  is  given  by 

A  6  =  2  -  ( . Y.  -  AT)  .  (240) 

where  .V.  is  the  number  of  zeros  and  .V.  is  the  number  of  poles  (counting 
multiplicities)  enclosed  by  T. 

from  the  above  theorem,  if  fiz)  has  only  zeros  in  the  region  enclosed  by  T  .  then 
every  constant  phase  contour  associated  with  a  zero  crosses  T. 

An  analytic  function  f(z)  can  be  expressed  as 


/(•')  =  \  I]'  +  CM/)]2  expO'O)  ,  (2411 

where  0  is  the  phase  of  fiz)  given  by 


0  =  arg  [./(.-)] 


(242) 
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If  l(z)  does  not  have  any  pole  or  branch  point,  tlien  the  constant  phase  curve.  0  =  0C. 
radiating  from  a  zero  of  f(z)  must  cross  the  closed  contour.  T,  enclosing  the  zero  at  least 
once.  In  addition,  no  other  zero  of  f(z)  may  be  on  this  phase  curve.  This  is  illustrated  in 
I  ’igure  5. 

A  line  of  constant  phase  which  crosses  the  contour  T  may  be  followed  until  it  leads 
to  a  zero  or  until  it  crosses  the  contour  again.  A  zero  of  f(z)  can  be  determined  from  the 
intersection  of  the  curves 

I  in  (/ )  =  0  ;  (0.  =  0  or  r)  (243) 

and 

Rc{f)  =  0  :  (d.  =  or  \  .  (2441 


ss 


Closed  contour 


C  onstant  Phase  Line  of  1(7)  (From  Kef.  I] 


B.  ROOT  SEARCH  fROCEDl  RE 


A  search  rectangle  is  set  up  in  the  complex  /-plane  so  that  the  simple  zeros  of  ft/) 
will  he  located  by  the  routine.  The  corners  of  this  rectangle  arc  denoted  by: 

•  rL  -  value  of  the  real  part  ofz  at  the  left  edge  of  the  rectangle. 

•  iR  -  value  of  the  real  part  ofz  at  the  right  edge  of  the  rectangle. 

•  ir  -  value  of  the  imaginary  part  ofz  at  the  top  edge  of  the  rectangle. 

•  ts  -  value  of  the  imaginary  part  ofz  at  the  bottom  of  the  rectangle. 

The  search  rectangle  is  divided  into  small  mesh  squares  of  length  J.  Restrictions  on 
the  size  of  i)  will  be  addressed  later. 

A  new  search  rectangle  is  generated  from  the  specified  search  rectangle  by  express¬ 
ing  the  edge  of  the  new  rectangle  in  terms  of  mesh  units.  One  mesh  unit  is  of  length  0. 
Thus 


and 


Jl  = 


Jr  — 


(245) 

(246) 

(247) 


h  - 


Im 


(24S) 


where  Int(x)  denotes  the  integer  part  of  x.  In  order  that  the  original  search  rectangle  falls 
completely  within  the  new  search  rectangle,  the  new  rectangle  is  made  one  mesh  unit 
larger  on  all  sides.  That  is. 


h 


j-'c  - 

I  Jl  ~ 


<L  ^  0 

1[  <  0 


(249) 


J 


k 


f  -'r  +  2 
(  Jr  +  1 


' R  ^  0 
'll  <  (» 


(25H) 
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JT  +  2 


ir  >  (» 

lr  <  0 


\h  ~  1 

-  2 


'B  2>  0 
//,  <  0 


The  search  rectangle  with  mesh  grid  is  shown  in  Figure  6.  The  pius  and  minus  signs 
at  the  corners  of  the  mesh  squares  are  the  signs  of  Im(l).  The  solid  curves  represent  the 
phase  line.  Im(fi  =  U.  while  the  dashed  curves  represent  the  phase  line.  Rc(0  =  0.  The 
zeros  of  liz)  are  located  at  the  crossing  of  the  lm(f)  =  0  and  Re(0  =  0  curves.  A  local 
mesh  coordinate  system  is  shown  in  Figure  7.  The  lower  left-hand  corner  of  the  mesh 
square  is  taken  as  the  origin  of  that  mesh  square's  local  coordinate  system.  The  value 
of  f( z )  evaluated  in  the  local  mesh  coordinates  of  the  k!h  mesh  square  is  denoted  by 
/(i)  (A  .  c>T  .  Thus,  the  values  of  f(z)  at  the  corners  of  mesh  square  1  arc 


f  '(0.  1)  =  f{JLS  +  jJrt) 


/  MU.  0)  =  /[./,_  a  +  j  (,/r  -  l)o]  . 


/0,(1.U)  =  /[</.  +  l)d  +  j(JT-  1)0] 


/  (1.  1)  =  /[</,  +  W+PrS']  ■  (256) 

A  basic  assumption  made  on  f(z)  is  that  lm(0  and  Re(f)  are  linear  functions  of  the 
local  mesh  coordinates,  y>  and  e>.  .  This  means  that  in  the  mesh  square  k.  f(z)  can  be 
expressed  as  [Ref.  lj 

fa)(/k.  )  =  uk  +  bk  ok  +  rkyk  +  <Ik  wk  /..  . 


(257) 


riguie  6.  Illustration  of  search  rectangle  with  niesli  grid  (I'rom  Ref.  I] 

where 


Ok  =  /''■’(  0,0)  , 

(258) 

'k  =  /(V, 1)  -  /(,)(«,  0)  , 

(259) 

ck  =  /WH,0)  -/<A)( 0,0) 

(260) 

dk  =  /(*>.0)  +/(A)(I.I)  ~/(A)(0,  I)  -/W(l,0)  .  (261) 

This  assumption  puts  an  upper  bound  on  the  size  of  the  mesh  unit,  5.  I  he  mesh 
unit  should  be  sufficiently  small  to  cnsuic  that  the  linear  variation  of  l(z)  holds.  How¬ 
ever,  ()  should  not  be  excessively  small  to  reduce  computer  run  time. 
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Figure  7.  Local  mesli  coordinate  system  [From  Ref.  1J 


I  lie  zeros  of  f[z)  are  searched  by  examining  the  imaginary  part  of  f(z)  counter¬ 
clockwise  around  the  rectangular  contour.  For  the  example  shown  in  Figure  6,  the  edges 
of  the  search  rectangle  arc  examined  for  crossing  of  the  phase  curve,  lm(f)  =  0,  starting 
with  mesh  square  1,  edge  1.  The  values  hn[f  w  (0,  1)]  and  //»[/  (l)  (0,  ())]  are  first  com¬ 
puted.  If  they  arc  of  the  same  sign,  then  there  will  be  no  lm(l)  =  0  curve  passing  through 
edge  1. 

Mesh  square  2  is  investigated  next.  I  he  values  oT  /m[/  <J)(0,  1)J  and  /;»[/  m  (0,  0)] 
arc  computed.  1  he  value  of  7/u[/ ,J)  (0,  I)]  is  identical  to  //»[/  (l)  (0,  OJ]  .  It  the  signs 
of  Im\_J' (J)  (0,  1)J  and  //?»[/'  <2>  (U,  0)j  arc  opposite  as  shown  in  Figure  6,  then  the 
lm(f)  =  0  curve  enters  mesh  square  2  along  edge  I.  The  crossing  point  is  given  by  [Ref. 

II 

*2  -  0  (262) 
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and 


hnjf(1){  0,0)] 

(0,  1)]  -  7m[/  (2)  (0,0)] 


(263) 


The  next  step  is  to  determine  if  a  zero  is  within  mesh  square  2.  The  values  of 
/m[/  0>(I,O)]  and  //«[/  c>  (1,  1)]  arc  first  computed.  Two  tests  must  be  made: 


1.  A  test  is  made  to  determine  if  there  is  one  or  two  lm(Q  =  0  curves  entering  and 

leaving  the  mesh  square.  If  Im\_f  <*!  (0.  0)]  and  hn\_f  w  (1,  1)]  are  of  the  same  sign, 
c.g.,  minus,  and  //#*[/*  ( 0.  1)]  and  lm\_f  <*>  ( 1 , 0)]  both  have  opposite  signs,  c.g., 

plus,  then  there  are  two  lines  of  lm(  0  —  0  entering  and  leaving  the  square.  This 
condition  is  depicted  in  mesh  square  14  of  Figure  6.  Otherwise,  there  is  only  one 
line  of  Jm(f  )  =  0  entering  and  leaving  the  mesh  square.  This  is  illustrated  in  mesh 
square  2  of  Figure  6. 

2.  A  test  has  to  be  performed  to  determine  if  there  is  at  least  one  Re(0  =  0  entering 

and  leaving  the  mesh  square.  If  (0. 0)]  ,  (0.  1)]  ,  Rc[J'w  (1. 0)]  and 

M/“>(1,0]  arc  l^lc  samc  !^cn  l^ere  ’s  n0  *'nc  ^c(0  =  0  passing 

through  the  mesh  square.  Otherwise,  there  is  at  least  one  Rc(l)  =  0  curve  entering 
and  leaving  the  mesh  square. 


Note  that  the  above  tests  will  fail  if  there  arc  two  or  more  lm(f)  =  0  or  Re(f)~0 
crossing  the  mesh  square  at  the  same  edge.  As  shown  in  Figure  8,  two  1111(0  =  0  crossing 
tiic  samc  edge  of  a  mesh  square  results  in  a  zero  being  missed.  In  this  situation,  the  mesh 
size  has  to  be  reduced  in  order  to  resolve  the  two  zeros. 
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figuie  8.  Crossing  of  two  Iin(l)  =  0  at  the  same  edge  of  the  inesli  square 


I  lie  program  (see  subroutine  fzerox  in  Chapter  V)  also  checks  lor  the  presence  of 
more  than  one  zero  on  the  same  phase  line.  Such  a  situation  may  occur  as  in  figure  9. 
1  he  program  treats  such  an  occurrence  as  an  error  ami  will  proceed  to  reduce  the  mesh 
size.  1  he  root-search  will  start  all  over  again. 

One  can  infer  that  the  presence  of  more  than  one  zero  on  the  same  phase  line  is  due 
to  the  inability  to  resolve  two  phase  lines  crossing  a  mesh  square  at  the  same  edge.  1  hese 
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zeros  actually  lie  on  different  phase  lines  and  are  'legitimate'  7.eros.  It  is  the  author  s 
opinion  that  these  zeros  arc  acceptable  and  a  mesh-size  reduction  accompanied  by  a  ic- 
run  of  the  root  search  routine  unnecessary. 


Figure  9.  Two  zeros  found  on  (he  same  phase  line 


The  line  lm(f)  =  0  is  traced  from  mesh  to  mesh  until  it  exits  the  contour  rectangle. 
For  the  example  in  Figure  (3,  the  curve  All  is  followed  through  mesh  squares  2.  3,  7,  II, 
10,  14  and  13.  In  mesh  square  11,  both  lm(()  =  0  and  Rc(l)=0  pass  through  the  mesh 
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square.  The  point  of  intersection  can  be  determined  from  Liquation  25".  The  math¬ 
ematical  details  are  presented  in  section  C. 

The  contour  AB  is  followed  out  of  mesh  square  11  into  mesh  squares  10  and  14. 
In  mesh  square  14,  there  are  two  lm(fi=0  curves  and  one  Rc(f)  =  0  curve.  A  crossing 
point  is  chosen  to  be  a  zero  of  fiz)  if  it  lies  on  the  lm(0  =  0  curve  currently  being  fol¬ 
lowed.  Therefore,  when  the  contour  AB  is  being  traced,  the  zero  in  mesh  square  14  will 
be  ignored.  However,  this  zero  will  be  located  when  contour  CD  is  being  traced. 

When  contour  AB  exits  the  search  rectangle  at  point  B.  the  value  of  the  leading 
edge  (counterclockwise)  of  the  mesh  square  is  noted.  In  this  example,  the  value  of  the 
leading  edge  is/'15’  (0,  1).  Identifying  the  exit  points  would  prevent  subsequent  re-enter¬ 
ing  of  the  search  rectangle  via  these  points. 

Alter  the  contour  AB  exits  at  B.  the  search  for  zeros  of  f(z)  is  continued  by  looking 
for  more  crossings  of  lm(f)  =  0  with  the  search  rectangle.  Since  mesh  square  2  was  the 
hist  mesh  square  that  was  checked  for  such  crossings,  the  next  mesh  square  to  be  ex¬ 
amined  will  be  mesh  square  3. 

The  next  crossing  of  lm(0  =  0  with  the  search  rectangle  is  at  mesh  square  12.  This 
contour  CD  is  again  followed  through  to  mesh  21  where  the  exit  point  D  is  noted.  The 
zero  along  CD  in  mesh  square  14  is  considered  as  an  acceptable  zero  of  f(z)  since  it  lies 
on  the  contour  Cl)  which  is  the  current  contour  being  traced.  The  rectangular  contour 
search  will  continue  as  described  above  until  the  search  gets  to  mesh  square  1.  edge  4 
where  it  stops. 

C.  DETERMINATION  OF  THE  COORDINATES  OF  ROOTS 

The  root  of  liz)  is  the  point  of  intersection  of  lm(f)  =  0  and  Rc(n  =  0.  This  can  be 
determined  from  Equation  257. 

By  equating  Im(f)  of  Equation  257  to  zero,  the  following  equations  are  obtained: 


//??(«/-)  +  Jni(bk)(Ok 
hn(ck)  +  Im(dk)u)k 


and 


4-  [ml  ('{.)/ L. 
hn(bk)  +  lm{dk)/k 


(264) 


(265) 


Similarly,  by  equating  Rc(f)  to  zero. 
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y.k  - 


Relak)  +  Re{h.  )(■>,, 

Rc[Cj,\  T  /u'(  </;)<•);. 


Reia,)  -)-  Re(ek )■/_,. 

0>k  =  '  Wk)  +  Wk)yk  '  (267) 

Note  that  Equations  264  and  265  are  equivalent;  so  arc  Equations  266  and  267. 
By  equating  the  right  hand  side  of  Equation  264  and  that  of  Equation  266.  the  following 
quadratic  equation  in  w,  results: 


where 


f’l* "  +  Q\V)lc  +  B  ,  —  0  , 


1']  =  Re  (fy)  -  ///?(/’/  )  fie  (<4) 


£4  =  Re  (ak)  Itn(i/k)  +  Re  (hk)  hn{<\) 

-  lm{a,)  Re  Ulk)  -  Re  {ck) 


JE,  =  Re  (ok)  Im(ck)  —  Im(ak)  Re  (ek)  . 

Similarly,  from  Equation  265  and  267.  the  quadratic  equation  in  yk  is  obtained: 


where 


y.k  +  Ql  '/k  +  2  =  0  . 


P2  =  Re(ck)  Im(iik)  —  Itn(ck)  Re(c/k)  , 

(J2  =  /?<’(«*)  7w(<4)  4-  /?<•(<-..)  /m(7y) 

-  Im(ok)  Re{i //;)  -  //;/(q.)  /M/y) 


H'2  =  Re(ak)}f?i{bk)  —  Im(ak)  Relbk)  . 
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According  to  the  program  source  code,  if 

!  I  <  I  A  I  .  (2"M 

then  the  root-finder  routine  uses  liquations  204  and  26S  to  solve  for  yt  and  o,  .  If  the 
condition  in  liquation  276  is  not  satisfied,  then  y .  and  w,  will  be  determined  from 
liquations  265  and  272. 

However,  it  should  be  noted  that  if  liquation  276  is  true  and  if  |  !  approaches 

zero,  one  of  the  solutions  of  liquation  26S  will  become  infinite.  By  the  same  token,  if 
liquation  276  is  not  true  and  if  |  /\  |  approaches  zero,  one  of  the  solutions  of  liquation 
272  will  become  infinite.  The  program  (sec  subroutine  quad  in  Chapter  V)  prevents  such 
occurrences  by  incorporating  a  test  such  that  the  program  returns  only  one  solution 
when  i  |  or  |  l\  \  approaches  zero.  i.e..  liquation  26S  or  272  reduces  to  a  linear 
equation. 

On  the  other  hand, .unless  both  | /’,  |  and  |  /',  |  become  very  small,  the  above 
problem  could  be  avoided  if  y ,  and  o.  are  determined  from  liquations  265  and  272  when 
Uquation  2~6  is  true;  otherwise  liquations  264  and  26$  are  used.  This  modification  is 
recommended. 

A  solution  of  liquations  265  and  272  or  264  and  268  will  be  valid  if  it  lies  within  the 
current  mesh  square,  i.e.. 


0  <  Xs  <  1 


(277) 


and 


0  <  w,  <  1  . 


(278) 


where  r.  =  {y,.  o>.)  is  the  solution  . 

In  addition,  the  solution  must  also  lie  on  the  line  lm(0  =  0  currently  being  traced. 
The  test  for  this  condition  is  derived  below. 

From  Equation  257,  lm(f)  =  0  and  Rc(f)  =  0  are  equations  of  an  equilateral  hyper¬ 
bola  with  vertical  and  horizontal  asymptotes  given  by 

Im  (  /’, ) 
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and 


Im  (<■  ) 
c  Im  (<<  ! 

Lot  {/c  o>L)  be  tlic  point  at  which  Im(0  =  U  enter  the  mesh  square.  If' 


l/.E  -  X(.)  (/,  -  /(•)  >  <> 


and 


(0>£  ~  t)C)  (Wj  -  CJC)  >  0  . 

then  the  solution  of f(z )  lies  on  the  current  line. 

Figure  H>  illustrates  the  ease  where  e  is  a  proper  solution  and  Figure 
case  where  r,  is  not  a  proper  solution  of  !'(/)  =  o. 


(2  SO) 


(2S1) 


(2$2) 


shows  the 
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D.  SOU  '1  ION  on  HE  QUADRATIC  EQUATION 

Procedures  lor  the  solution  of  the  quadratic  equation  are  required  for  solving 
Equations  2<>S  and  2~2. 

Consider  the  quadratic  equation  of  the  for  in 

ax'  +  2!\\  +  (•  =  <>.  ( 2  S  3 ) 

The  solution  is 

2h  ±  ^  4 h "  —  4we 


1  -  1  - 

\ 


I  :  ( 


=  —  [- 


(2S4) 


where 


--  .  (2S5) 

/»• 


I  fence. 


and 


\ 


I  4-  c 


Ifc  <7  1.  then 


(1  4  £)'  : 


(2S6| 


(2S7) 


(2SS) 


Thus. 


Factoring 


Let 


and 


In  general 


out 


^  'j  .  Equation  2S9  becomes 

2  b  [  2 

+-o(+-o 


or 

L  A2  J 


(2S9) 


(290) 


h  =  I 


(291) 


(292) 


F  = 


(+-')(+- 


(293) 


(294) 


os 


.V- 


(295) 


ana 


- ~  ~  {  /[•  4-  6  4-  /t  4* 

2  b 


f 


Once  a*,  is  found,  x2  can  be  found  from  Liquations  286  and  2S7.  By  adding  Equations 
286  and  287,  x2  is  given  by 

.v2  =  -  2—  -  x,  .  (296) 

E.  NEWTON  RAPHSON  ITERATION 

The  zeros  of  f( /)  found  by  the  Shellmun  and  Morfitt  routine  are  otily  approximate 
solutions.  I  he  Newton  Raphson  iteration  is  used  to  improxc  the  accuracy  of  the  zero 
locations.  The  new  root.  is  related  to  the  previous  root.  r  .  bv 

rM-:  =  :i  +  •  (-l)'T) 

where 


Ar, 


/U-;l 


(298) 


Iteration  of  Equation  297  is  continued  until  A:  is  smaller  than  a  pre-assigned  tolerance. 
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IV.  SOLUTION  OF  STOKES'  DIFFERENTIAL  EQUATION 
A.  STOKES'  EQUATION 

When  the  modified  refractivity  profile  is  modelled  by  a  scries  of  piecewise  linear 
segments  in  the  mathematical  formulation  of  wave  propagation  in  a  multilayer  tropos¬ 
pheric  waveguide  environment,  the  height  gain  function  in  each  layer  can  be  obtained 
from  solution  of  Stokes'  Equation. 


(-/,)  =  0 


:  «  =  1. 


(-991 


Solutions  of  Stokes'  Equation  are  commonly  given  in  terms  of  Airy  functions  or  in  terms 
of  modified  Hankel  functions  of  order  one-third.  Several  possible  solutions  exist  and  the 
choice  of  solution  is  dependent  on  numerical  and  physical  considerations.  Since  //»(</.) 
is  proportional  to  the  attenuation  rate  of  the  electromagnetic  signal,  physical  consider¬ 
ations  dictate  that  lm{q)  >  0  .  In  this  program,  the  height  gain  function  is  expressed 
as  a  linear  combination  of  £,(</  )  and  A,  (</)  which  arc  linearly  independent  for 
ImU] )  >  0  .  The  choice  of  /;,(</),  A ,((/.)  or  equivalently  {.-//(  —  </  c  5 ).  .hi  —  </  )i  is 
to  assure  that  at  least  one  of  the  two  functions  do  not  become  very  small  numerically  in 
the  upper  half  of  the  complex  (/  -plane.  Eor  the  top  layer,  i.c..  i=  1.  the  boundary  condi¬ 
tion  at  ;  -+  oo  requires  that  the  height  gain  function  represents  an  outgoing  wave. 
Hence  in  the  topmost  layer,  /:-(</)  or  equivalently  Ai(  —  <7  e~'2'  !)  is  chosen  to  satisfy  the 
radiation  condition. 

B.  EVALUATION  OF  AIRY  FUNCTIONS 
I.  General  Outline 

The  algorithm  for  evaluation  of  complex  Airy  functions  is  adopted  from  Z. 
Schultcn  and  D.G.Yl.  Anderson  [Ref.  16],  Evaluation  of  the  complex  function  over  the 
entire  /-plane  was  achieved  with  a  Taylor  series  and  generalised  Gaussian  quadrature 
method.  A  Taylor  series  is  used  only  for  a  small  region  near  the  origin  within  an  ellipse 
with  foci  (  -1.5S.  3.95)  and  (U.2S,  -2.11).  major  axis  of  9,0  and  minor  axis  of  0.4. 
Outside  the  ellipse,  the  Gaussian  quadrature  approximation  is  used.  The  elliptical 
boundary  was  chosen  to  satisfy  restrictions  imposed  by  computational  methods. 
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Analysis  of  the  Airy  functions  need  only  be  done  in  the  upper  hall'  complex 
plane  as  those  in  the  lower  half  plane  can  he  computed  by  the  following  conjugate 
relations: 

,!/(:)  =  Ai’(z’)  (300) 


and 


(/AH:)  _  dA/'lz') 

dz  dz 


(301) 


Therefore,  an  algorithm  can  be  developed  to  compute  Ah/)  for  o  <  arg  r  <  2, 7, '3  and 
by  conjugate.  Ai(/t  for  — 2r/3  <  arg  r  <  0  can  be  computed. 

For  2-/3  <  arg  r  <  4r/3,  the  following  connection  formula  of  Airy  functions 
can  be  used: 

/!/(.-)  =  ej~!3  Ai(zc~i:r:-  |  +  2  Ai  {zc  ;2~  3 )  .  (302) 


2.  Evaluation  bv  Taylor  Series  Expansion 

Tor  small  values  of  1/  |,  the  Airy  functions  are  evaluated  by  the  scries  expansion 
[Ref.  1  and  4j. 

Ailz)  =  (',/(.-)  -  C;g(r)  (303) 

and 

,l/'(r)  =  O/'U)  -  Gg'(r)  .  (304) 


where 


/(r)  = 


?  !■: 


(3  AT 


(305) 
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The  number  of  terms  required  for  convergence  of  Uk  scries  is  determined  by  a  tolerance 
factor.  £.  such  that 


/!-] 


<= 0 


(311) 


where 


Ai{z) 


n 


*'=') 


(312) 


Similiur  test  for  comergence  is  implemented  for  Ai'{z)  .  In  the  program,  r,  is  taken  to  he 
the  machine  epsilon,  i.e.. 


3.  Evaluation  by  Generalised  Gaussian  Quadrature  Method 

Integral  expressions  for  Ai(z)  is  given  by  (Ref.  1('| 


....  I  -in  '  plxuix 

At  (r)  =  - —  r  t’  •  w  — - 

■>  *  -r  .v 

**  \  %  ii 


for  |  e  (  >  0  .  |  arg  z  j  <  —~ 


where 


/>(  a)  = 


„-iT  -,-if 


a  .v 


The  function.  />{.v).  is  non-negative  and  exponentially  decreasing.  The  integral  in 
Equation  314  can  be  evaluated  by  the  generalised  Gaussian  quadratic  approximation. 


l*-»3  rA, 

/T.V)  </.V  V  U; 

:  +  .v  ~  Zj  :  + 


where  w  is  the  quadrature  weight  and  .r,  is  the  abscissa.  Hence 


AH:) 


—  r  e  e  ,  > 
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and 


The  quadrature  weights,  ir,  .  and  abscissae,  .v  ,  had  been  calculated  by  Z.  Schultcn  ct  al. 
[Ref.  16)  and  are  given  in  Tables  1!  to  IV.  I  hese  tables  had  been  used  for  computation 
of  single  precision  values  up  to  seven  significant  digits. 
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Table  -•  NUMBER  OF  TERMS  REQUIRED  IN  QUADRATURE  FORMULA 
[AFTER  REE.  16j  _ 


2 

nrg  z 

Ai(z) 

2  rr 

2  terms 

1 

\z\  >  11 

<  3 

2tt 

5  < 

M  <  ll 

77 

4  terms 

2.5  < 

Ul  <  5 

4  terms 

Table  3.  4  TERM  GENERALISED  GAUSSIAN  INTEGRATION  FOR  AIRY 

FUNCTIONS  [FROM  REF.  16 J 


Table  4.  2-TERM  GENERALISED  GAUSSIAN  INTEGRATION  TOR  AIRY 
FUNCTIONS  [FROM  REF.  I6[ 


1.0592469382112378 
3.6800601 8661 53044(-02) 


3.1927194042263958(  — 02) 
9.6807280595773604(— 01) 


C.  EXTENDED  COMPLEX  ARITHMETIC 

T  he  extended  complex  arithmetic  was  introduced  by  Baumgartner  in  XWYG  to 
handle  complex  numbers  of  large  magnitude  to  avoid  overflow  problem. 

It  was  discussed  in  Section  D  of  Chapter  11  that  for  large  values  of  |  q  |  .  the  mag¬ 
nitude  of  A, (<7)  and  k2{q)  may  become  exponentially  large  or  exponentially  small. 

Numerical  evaluation  of  the  modal  function  of  Equation  166  can  easily  yield  com¬ 
plex  numbers  with  magnitudes  as  large  as  lo-|0,)0  or  as  small  as  10~inpn .  Numbers  of  these 
magnitudes  arc  outside  the  numerical  limits  of  most  computers. 

To  overcome  this  problem,  each  complex  function,  k}(q)  (or  k2  [q)),  is  represented 
by 

A,  =  kxe*  ,  (320) 

where  A,  is  the  complex  amplitude  and  <f>  is  the  real  exponent.  In  the  program.  A,  is 
"normalized”  such  that 

c~‘  <  max  [  |  Pc  (A,)  |  ,  |  hn{kx)  |  ]  <  e+l 


(321) 


and  4>  has  integer  values.  By  evaluating  the  complex  amplitude  and  real  exponent  sepa¬ 
rately.  overflow  or  underflow  problems  are  avoided.  Functions  expressed  in  the  form 
as  in  Equation  320  are  known  as  extended  complex  numbers.  Numerical  manipulation 
such  as  addition  and  multiplication  involving  extended  complex  numbers  are  performed 
in  the  following  manner.  If  r,  =  r,  and  zl  -Zj^  are  two  extended  complex  num¬ 
ber,  then  their  product  is 


2  =  2, 


A  <5 

=  z  e 


(322) 


where 


A  A  A 

z  -  Zj  z2 


and 


4*  —  4*i  +  4*2  • 

The  sum  of  two  extended  complex  numbers  is  given  by 

2  -  Z,  +  Z2 


(323) 


(324) 


(325) 


where 

z-z,  +  (326) 

and 

4*  -  4*i  ;  4*i  4>i  (327) 

or 

2-22  +  2,^-^  ;  4>,  <  <t>,  (328) 

and 

4*  —  4>i  ’  4>\  <  4*i  ■  (329) 
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V.  DESCRIPTION  OF  MLAYER  PROGRAM  ELEMENTS 

A  short  description  of  the  functions  of  cael:  program  element  in  MLAYER  is  pro¬ 
vided  in  this  Chapter.  The  program  had  been  tested  and  ran  successfully  on  the  Micro¬ 
soft  FORTRAN  version  5.0  compilers  and  Microsoft  C  compilers  version  5.0.  The 
maximum  number  of  layers  and  modes  are  35  and  150  respectively.  These  limits  are  the 
upper  bounds  for  operation  on  personal  computers  with  640  kilobytes  ormain  memory. 

A.  MAIN 

i.  Description 

Main  is  the  controlling  program  element  for  MLAYER.  It  calls  appropriate 
subroutines  to  calculate  modal  eigenvalues,  mode  sum,  path  loss  and  radio  horizon  dis¬ 
tance.  Main  also  calls  subroutine  wvgstdin  to  read  in  data  from  input  file,  iilcin.  Three 
output  files  are  opened,  namely,  filein.out,  filein.cig  and  filein.plt. 

I  dein.out  provides  a  documented  output  of  input  parameters  and  a  list  of  mo¬ 
dal  eigenvalues  found  in  each  search  rectangle.  Filein.out  also  contains  a  list  of  eigen¬ 
values  in  increasing  order  of  real  parts  with  the  corresponding  values  in  the  0-plane  and 
respective  attenuation  of  each  mode  in  decibels  per  kilometer. 

Filein.cig  is  an  exact  replica  of  the  input  file  except  that  the  parameter,  mvfile, 
is  set  to  one  and  it  contains  a  list  of  eigenvalues  found.  Once  filein.cig  is  created,  it  can 
be  used  as  the  input  file  to  re-run  MLAYER  for  various  transmitter; receiver  range  and 
height  configurations. 

Filein.plt  contains  eight  columns  of  data,  namely  range,  transmitter  height,  re¬ 
ceiver  height,  coherent  mode  sum.  incoherent  mode  sum,  coherent  path  loss,  incoherent 
path  loss  and  radio  horizon  distance.  The  data  are  arranged  in  this  manner  to  facilitate 
graph  plotting. 

Samples  of  filein,  filein.out,  filein.eig  and  filein.plt  are  given  in  Appendix  B. 

In  addition,  main  also  computes  the  following: 

•  The  derivative  of  the  modified  rcfractivity  with  respect  to  z  for  each  layer. 

•  The  derivative  of  the  atmospheric  absorption  with  respect  to  z  for  each  layer. 

•  The  absorptivity  and  its  derivative  with  respect  to  z  for  each  layer. 

•  The  derivative  of  q  with  respect  to  z  for  each  layer. 

•  The  modal  attenuation  rate  for  each  mode. 
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•  The  coefficients  required  for  computation  of  surface  roughness. 

•  The  eigenvalues  in  the  0-plane.  The  eigenvalues  in  the  0-plane  arc  related  to  those 
in  the  t/,  ,-plane  by 

0  =  sin  [v  1  -  p2  ]  .  (330) 

The  conversion  to  0-plane  is  done  solely  for  comparison  of  results  with  other  models 
where  root-search  was  performed  in  the  0-plane. 


2.  Inputs 


Input  data  are  read  in  by  subroutine  wvgstdin  in  the  following  order: 


•  filein 

•  nifile 

•  fqmzin 

•  delfq 

•  nfreq 

•  mpol 

•  aloss 

•  seatmp 

•  seaslt 

•  iflgab 

•  airtrnp 

•  rh 

•  wgpm3 

•  rmsbht 

•  ztinit 

•  delzt 

•  nzt 


:  Pathname  of  input  file. 

:  If  mfile  =  0.  read  input  data  and  calculate  eigenvalues. 

I f  mfilc  =  I,  read  input  data  and  eigenvalues. 

:  Initial  frc'p'.ncy  in  megahertz. 

:  Frequency  increment  in  megahertz. 

:  Number  of  frequencies  for  which  modes  will  be  found. 

:  If  mpol  =  0,  wave  is  horizontally  polarized. 

If  mpol=  1.  wave  is  vertically  polarized. 

:  Maximum  rate  of  attenuation  in  decibels  per  kilometer  of  modes 
that  will  be  found. 

:  Seawater  temperature  in  degrees  Celsius. 

:  Seawater  salinity  in  grams  of  salt  per  kilogram  of  seawater. 

:  If  ifigab  =  0,  calculate  the  height  invariant  atmospheric  absorption. 

If  iflgab  =  1,  read  in  array  containing  values  of  atmospheric 
absorption  at  layer  boundaries. 

:  Air  temperature  in  degrees  Celsius  at  reference  height,  zref. 

Airtrnp  is  used  for  calculating  the  height  invariant  atmospheric 
absorption  when  iflgab  =0. 

:  Relative  humidity  in  percent  at  reference  height.  Used  for  calculating 
the  height  invariant  atmospheric  absorption  when  iflgab  =  0. 

:  Water  concentration  in  grams  per  cubic  meter  at  reference  height. 
Used  for  calculating  the  height  invariant  atmospheric  absorption 
when  iflcab  =  0. 

:  Rms  bump  height  in  meters  of  the  sea  surface. 

:  Initial  height  of  transmitter  in  meters. 

:  Height  increment  of  transmitter  in  meters. 

:  Number  of  transmitter  heights  at  which  field  strength  is  to  be 


79 


computed. 

•  zrinit  :  Initial  height  of  tcceivcr  in  meters. 

•  delzr  :  Height  increment  of  receiver  in  meters. 

•  nzr  :  Number  of  receiver  heights  at  which  field  strength  is  to  be  computed. 

•  xinit  :  Initial  range  separation  of  the  transmitter  and  receiver  in  kilometers. 

•  dclx  :  Increment  in  range  separation  in  kilometers. 

•  nx  :  Number  of  range  separation  values  at  which  field  strength  is  to  be 

computed. 

•  zref  :  Reference  height  at  which  air  temperature,  relative  humidity  and 

water  concentration  are  specified. 

•  nzlayr  :  Number  of  linear  segments  used  to  model  the  modified  refractivity 

profile. 

•  zi(j)  :  Array  containing  values  of  heights  in  meters  at  layer  boundaries. 

•  zim(j)  :  Array  containing  values  of  modified  refractivity  at  layer  boundaries. 

•  zigab(j)  :  Array  containing  values  of  atmospheric  absorption  in  decibels 

per  kilometer  at  layer  boundaries. 

3.  Outputs 

The  outputs  to  filein.out  (logical  file  16)  are: 

•  nrmodc  :  Number  of  modes  found. 

•  qcigcn(j)  :  Array  containing  values  of  modal  eigenvalues  found  in  the  complex 

-plane. 

•  theta  :  Complex  eigenangle  referenced  to  ground  level  for  each  mode. 

•  atnu  :  Rate  of  attenuation  in  decibels  per  kilometer  for  each  mode. 

The  outputs  to  lilcin.eig  (logical  file  17)  are  : 

•  nrmode  :  Number  of  modes  found. 

•  qeigen(j)  :  Array  containing  values  of  modal  eigenvalues  found  in  the  complex 

qx  ,-plane. 

The  outputs  to  filein.plt  (logical  file  15) 

•  ccms  :  Coherent  mode  sum  field  strength  relative  to  free  space  in  decibels. 

•  cims  :  Incoherent  mode  sum  field  strength  relative  to  free  space  in  decibels. 

•  ccpl  :  Coherent  mode  sum  path  loss  in  decibels. 

•  eipl  :  Incoherent  mode  sum  path  loss  in  decibels. 
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4.  Calling  subroutine 

•  None 

5.  Subroutines  called 

•  wvgstdin 

•  chkntod 

•  seah2o  ;  commented  out 

•  fndinod 

•  casin 

•  ao2h2o  ;  disabled 

•  modsum 

•  dhoriz 

•  function  ibstrip 

The  call  to  subroutine  scah2o  was  commented  out  in  revision  3.0  of  the  pro¬ 
gram.  To  use  the  subroutine,  remove  the  comment  statement.  Subroutine  ao2h2o  was 
disabled  by  setting  iflgab  to  one.  To  enable  ao2h2o.  set  illgab  to  zero. 

6.  Common  block  areas 

•  com  1 

•  com  2 

•  com  3 

•  pap  1 

•  datum 

B.  SUBROUTINE  ABCOEF 

1.  Description 

* 

Subroutine  abeoef  calculates  the  coefficients,  A,  and  of  the  height  gain 
function  given  by  Equations  25  and  2S.  Computations  of  the  coefficients  arc  done  in 
extended  complex  arithmetic.  The  normalizing  integral  of  Equation  77  is  also  evaluated. 
The  coefficients  arc  then  normalized  by  the  following  relation: 


SI 


By  noting  that 

it;  d; 

dl'm  (lil\,\  '  (/l'm 

(332) 

‘t;  /  k  \  2/3  2/<„, 

«/«/,.!  ‘  V  *1  )  k2 

it  can  be  seen  that  the  denominator  of  Equation  331  is  identical  to  that  of  liquation  77. 

In  calculating  the  coefficients,  the  decision  to  integrate  up  or  down  is  deter¬ 
mined  by  the  rate  of  attenuation  given  by 

art  =  /?<■[/)„,]  x  ?n  log  exit)*  dB'jkm  .  (333) 

If  the  attenuation  rate  is  greater  than  0.1  dB  km.  upward  integration  is  performed;  oth¬ 
erwise  downward  integration  is  carried  out. 

At  the  same  time,  the  subroutine  also  checks  for  evanescence  using  liquations 
153  through  156  for  upward  integration  and  liquations  162  through  165  for  downward 
integration. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  abcoef  (zero.  m). 

3.  Inputs 

The  input  variables  are 

•  zero  :  An  eigenvalue  in  the  <7,  ,-planc  . 

•  m  :  Mode  number. 

4.  Outputs 

The  output  variables  are 

•  acoefa(i.m) :  Two  dimensional  array  element  of  complex  amplitude  of  A  in  the 

layer  and  m'"  mode. 

•  acocfc(i.m) :  1  wo  dimensional  array  element  of  real  exponent  of  A  in  the 

i":  layer  and  ni  mode. 
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•  bcoefati.mi:  Two  dimensional  array  element  of  complex  amplitude  of  B  in  the 

layer  and  nr  mode. 

•  bcoefe(i.m) :  Two  dimensional  array  element  of  real  exponent  of  B  in  the 

/'  layer  and  ui  mode. 

5.  Calling  program  element 

•  modsum 

6.  Subroutines  called 

•  XCdnj 

•  xcadd 

7.  Common  block  areas 

•  com  1 

•  com  2 

•  pap  1 

•  pap  2 

C.  Sl'BROLTINC  ADDX 
1.  Description 

This  subroutine  adds  two  extended  real  numbers  according  to  the  following 

equation: 


r  =  r,  +  .  (334) 

where 

r  =  za  exp  (r,i  .  (335) 


ri  —  zio  exp  (::,<,) 


(336) 


and 


~i  =  z2a  C.\p(r2c)  .  (337) 

Calling  statement 

The  subroutine  is  called  by  the  statement 

Cali  addx  (/a.  /e.  zla.  /le.  z2a.  z2e). 
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3.  Inputs 

The  input  variables  are 

•  zla  :  Real  amplitude  of  the  extended  real  number,  r,  . 

•  zlc  :  Real  exponent  of  the  extended  real  number,  r  . 

•  z2a  :  Real  amplitude  of  the  extended  real  number,  r,  . 

•  z2e  :  Real  exponent  of  the  extended  real  number,  r,  . 

4.  Outputs 

The  output  variables  arc 

•  za  :  Real  amplitude  of' the  extended  real  number,  z. 

•  ze  :  Real  exponent  of  the  extended  real  number,  z. 

5.  Calling  program  element 

•  modsum 

6.  Subroutine  called 

•  normre 

7.  Common  block  area 

•  none 

D.  SUBROUTINE  A02H20 

1.  Description 

This  subroutine  calculates  the  atmospheric  absorption  coefficient  in  decibels 
per  kilometer  due  to  the  quantum  mechanical  resonances  of  oxygen  molecules  and  water 
vapour.  Computed  values  are  good  for  frequencies  in  the  range  1  Ghz  to  1UU0  GHz. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  ao2h2o  (atmabs,  fliz,  tc,  rh,  zm,  wgpm3). 

3.  Inputs 

The  input  variables  are 

•  (hz  :  Frequency  in  Hertz. 

•  tc  :  Temperature  in  degrees  Cclcius. 

•  rh  :  Relative  humidity  in  percentage. 

•  zm  :  Height  above  sea  level  in  meters. 
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•  \vgpm3 


:  Water  concentration  in  grains  per  cubic  meter  at  reference  height. 


4.  Output 

I  he  output  variable  is 

•  atmabs  :  Atmosphetric  absorption  coefficient  in  decibels  per  kilometer. 

5.  Calling  program  element 

•  main 

6.  Subroutine  called 

•  none 

7.  Common  block  area 

•  none 

E.  SUBROUTINE  CIIKMOD 

1.  Description 

This  subroutine  checks  for  modal  eigenvalues  with  zero  value  and  discard  them. 
A  zero-value  eigenvalue  may  sometimes  lead  to  a  dividc-by-zero  error  in  the  mode  sum 
calculation. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  chkmod  (qeigen.  nrmode). 

3.  Inputs 

The  input  variables  are 

•  qeigen(j)  :  Array  containing  location  of  modal  eigenvalues  in  the  complex 

t/,., -plane. 

•  nrmode  :  Number  of  modes. 

4.  Outputs 

The  output  variables  are 

•  qcigcn(j)  :  Array  containing  location  of  non-zero  eigenvalues  in  complex 

<7,  ,-planc. 

•  nrmode  :  Number  of  modes. 


85 


5.  Calling  program  element 

•  main 

(>.  Subroutine  called 

•  none 

7.  Common  block  area 

•  none 

F.  SUBROUTINE  DETR 

1.  Description 

This  subroutine  evaluates  the  modal  function  determinant  of  Equation  47  in 
extended  complex  arithmetic.  The  algorithm  is  discussed  in  Section  E  of  Chapter  II. 

2.  (  ailing  statement 

I  he  subroutine  is  called  by  the  statement 
Call  detr  (deta.  dote,  n/layr.  d.  do.  dpi,  dple,  dm!,  dmle.  dp2.  dp2e.  dm2.  dm2c). 

3.  Inputs 

1  he  input  variables  arc 

•  n/layr  :  Number  of  linear  segments  used  to  model  the  modified  rcfractivity 

profile. 

•  d  :  Array  containing  complex  amplitudes  of  elements  along  the  main 

diagonal  of  the  modal  function  determinant. 

•  d '  :  Array  containing  real  exponents  of  elements  along  the  main 

diagonal  of  the  modal  function  determinant. 

•  dpi  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

one  above  the  main  diagonal  of  the  modal  function  determinant. 

•  dple  :  Array  containing  real  exponents  of  elements  along  the  diagonal  one 

above  the  main  diagonal  of  the  modal  function  determinant. 

•  dml  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

one  below  the  main  diagonal  of  the  modal  function  determinant. 

•  dmle  :  Array  containing  real  exponents  of  elements  along  the  diagonal  one 

below  the  main  diagonal  of  the  modal  function  determinant. 

•  dp2  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

two  above  the  main  diagonal  of  the  modal  function  determinant. 

•  dp2e  :  Array  containing  real  exponents  of  elements  along  the  diagonal  two 

above  the  main  diagonal  of  the  modal  function  determinant. 

•  dm2  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

two  below  the  main  diagonal  of  the  modal  function  determinant. 

•  dm2e  :  Array  containing  real  exponents  of  elements  along  the  diagonal  two 


below  the  main  diagonal  of  the  modal  function  determinant. 


4.  Outputs 

The  output  variables  are 

•  deta  :  Complex  amplitude  of  the  modal  f  unction. 

•  dete  :  Real  exponent  of  the  modal  function. 

5.  Calling  program  element 

•  fctvlx 

6.  Subroutine  called 

•  xcadd 

7.  Common  block  area 

•  none 

G.  SUBROUTINE  DHOR1Z 

1.  Description 

The  subroutine  calculates  the  radio  horizon  distance  between  the  transmitter 
and  receiver.  The  radio  horizon  distance  in  kilometers  is  given  bv 


d  —  \  ~  •  y  “T  U  :R  M  \  z  V  ]  1  ( 1 
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where  a  is  the  radius  of  the  earth,  and  and  zT  are  the  receiver  and  transmitter  heights 
respectively. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  dhoriz  (dhz,  zrevr.  zxmtr). 

3.  Inputs 

The  input  variables  are 

•  zrevr  :  Receiver  height  in  meters. 

•  zxmtr  :  Transmitter  height  in  meters. 

4.  Output 

I  lie  output  variable  is 

•  dhz  :  Radio  horizon  distance  in  kilometers. 


5.  Calling  program  element 

•  main 

6.  Subroutine  called 

•  none 

7.  Commom  block  area 

•  none 

H.  SUBROUTINE  DXDETR 

1.  Description 

This  subroutine  evaluates  the  modal  function  determinant  and  its  derivative  in 
extended  complex  arithmetic.  Values  of  the  modal  function  determinant  and  its  deriva¬ 
tive  are  required  by  the  Newton-Ruphson  root-finding  routine.  The  methods  for  evalu¬ 
ating  the  modal  function  determinant  and  its  derivative  arc  discussed  in  Section  E  of 
Chapter  II. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  dxdetr  (deta.  dctc.  dxdeta.  d.xdctc.  nzlayr.  d.  de.  dpi.  dple,  dm!,  dmlc.  dp2.  dp2e. 
dm2.  dm2e.  dxd.  dxde.  dxdpl.  d.xdplc.  d.xdml.  dxdmle.  dxdp2.  dxdp2e.  dxdm2.  dxdm2c). 

3.  Inputs 

The  input  variables  are 

•  n/layr  :  Number  of  linear  segments  used  to  model  the  modified  refractivity 

profile. 

•  d  :  Array  containing  complex  amplitudes  of  elements  along  the  main 

diagonal  of  the  modal  function  deteinmuuii. 

•  de  :  Array  containing  real  exponents  of  elements  along  the  main 

diagonal  of  the  modal  function  determinant. 

•  dpi  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

one  above  the  main  diagonal  of  the  modal  function  determinant. 

•  dple  :  Array  containing  real  exponents  of  elements  along  the  diagonal  one 

above  the  main  diagonal  of  the  modal  function  determinant. 

•  dml  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

one  below  the  main  diagonal  of  the  modal  function  determinant. 

•  dmlc  :  Array  containing  real  exponents  of  elements  along  the  diagonal  one 

below  the  main  diagonal  of  the  modal  function  determinant. 

•  dp2  :  Array  containing  complex  amplitudes  of  elements  along  the  diagonal 

two  above  the  main  diagonal  of  the  modal  function  determinant. 


dp2e  :  Array  containing  real  exponent''  of  elements  alone  the  dfieonu!  two 

above  the  mam  diagonal  of  the  modal  Junction  determinant. 

dm2  :  Arrax  containing  complex  amplitudes  of  elements  alone  the  diagnna 

two  below  the  mam  diagonal  o!  the  modal  function  determinant. 

dm2e  :  Array  eontainine  real  exponents  cf  elements  alone  1 1 : e  diagonal  two 
below  the  main  diagonal  of  the  modal  function  determinant. 

dxd  :  Arrax  eontainine  complex  amplitudes  of  the  derivatives  of  the 

elements  alone  the  main  diagonal  of  the  modal  function  determinant 

dxde  :  Arrax  containing  real  exponents  of  the  derivatives  oi  the  element'' 

along  the  main  diagonal  of  the  modal  function  determinant. 

dxdpl  :  Array  containing  complex  amplitudes  of  the  derivativ  es  of  the 

elements  along  the  diagonal  one  above  die  main  diagonal  of  the 
modal  function  determinant. 

dxdpl e  :  Array  containing  real  exponents  of  the  derivatives  of  the 

elements  along  the  diagonal  one  above  the  main  diagonal  of  the 
modal  function  determinant. 

dxdml  :  Array  containing  complex  amplitudes  of  the  derivatives  of  the 

elements  along  the  diagonal  one  below  the  main  diagonal  of  the 
modal  function  determinant. 

dxdmlc  :  Array  containing  real  exponents  of  the  derivatives  of  the  elements 
along  the  diagonal  one  below  the  main  diagonal  of  the  model 
function  determinant. 

d.xdp2  :  Array  containing  complex  amplitudes  of  the  derivatives  of'  the 

elements  along  the  diagonal  two  above  the  main  diagonal  of  the 
modal  function  determinant. 

d.xdp2c  :  Array  containing  real  exponents  of  the  derivatives  of  the  elements 
along  the  diagonal  two  abov  e  the  main  diagonal  of  the  modal 
function  determinant. 

d\dm2  :  Array  containing  complex  amplitudes  of  the  derivatives  or,.h? 

elements  along  the  diagonal  two  below  the  main  diagonal  of  the 
modal  function  determinant. 

d\dm2c  :  Array  containing  real  exponents  of  the  derivatives  of  the  elements 
along  the  diagonal  two  below  the  main  diagonal  of  the  modal 
function  determinant. 


4.  Outputs 

The  output  variables  are 

deta  :  Complex  amplitude  of  the  modal  function, 

dote  :  Real  exponent  of  the  modal  function. 

dxdcta  :  Complex  amplitude  of  live  derivative  of  the  modal  function, 
dxdcte  :  Real  exponent  of  the  derivative  of  the  modal  function. 


5.  Calling  program  element 

•  fdfdtx 

6.  Subroutine  called 

•  xcadd 

7.  Commom  block  area 

•  none 

I.  SUBROUTINE  FCTVLX 

1.  Description 

This  subroutine  calculates  the  values  of  the  elements  in  the  mat  rix  of  the  modal 
function  of  Equation  47  and  invokes  subroutine  detr  to  evaluate  the  modal  function.  The 
subroutine  shifts  the  imaginary  part  of  qlA  by  qshift  so  that  the  real  axis  will  not  fall  on 
the  mesh  line  of  the  search  rectangle  set  up  by  the  Shellman-Morlitt  routine.  This  allows 
the  zeros  on  the  real  axis  to  be  located.  It  is  known  that  zeros  near  or  on  the  negative 
real  axis  contribute  significantly  to  the  fields  within  a  duct. 

However,  the  shift  implies  that  the  Shellman-Morfitt  routine  will  search  for 
zeros  of  ||A  (z,.,)||  instead  of  ||A(<7U)II  where 

'i.i  -  <7i.i  -  j qshift  .  (339) 

To  compensate  for  the  error,  the  shift  back  to  qul  is  done  in  subroutine  fndmod  for  all 
zeros  found. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  fctvlx  (ql  lin,  deta,  dctc). 

3.  Input 

The  input  variable  is 

•  qllin  :  Location  in  the  complex  <?,  , -plane  where  the  modal  function  is  to 

be  evaluated. 

4.  Outputs 

The  output  variable  are 

•  deta  :  Complex  amplitude  of  the  modal  function. 

•  dctc  :  Real  exponent  of  the  modal  function. 
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5.  Calling  program  element 

•  fzerox 

6.  Subroutines  called 

•  xcdai 

•  xcadd 

•  dctr 

•  surf 

7.  Common  block  areas 

•  cotnl 

•  com2 

•  conij 

•  com4 


J.  SUBROUTINE  FBFDTX 

1.  Description 

This  subroutine  calculates  the  values  of  the  elements  in  the  matrices  of 

Il4(c/u)l!  and  ~ —  ||A(</,,,)!|  and  then  invokes  subroutine  dxdctr  to  evaluate  the  modal 
f-  lh.\ 

function  and  its  derivative.  The  subroutine  also  shifts  the  imaginary  part  of  qu  by  qshift 
as  in  subroutine  fctvlx.  It  is  similarly  compensated  by  a  "shift- back"  in  subroutine 
fndmod. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  fdfdtx  (qllin.  deta.  dete,  dxdeta,  dxdete). 

3.  Input 

The  input  variable  is 

•  ql  lin  :  Location  in  the  complex  <7,.,-plane  where  the  modal  function  and  its 
derivative  is  to  be  evaluated. 


4.  Outputs 

The  output  variables  are 

deta  :  Complex  amplitude  of  the  modal  function. 

dete  :  Real  exponent  of  the  modal  function. 

dxdeta  :  Complex  amplitude  of  the  derivative  of  the  modal  function. 

dxdete  :  Real  exponent  of  the  derivative  of  the  modal  function. 


91 


5.  Calling  program  element 

•  f/crox 

6.  Subroutines  called 

•  \cdai 

•  xcadd 

•  dxdetr 

•  surf 

7.  Common  block  areas 

•  coml 

•  coin2 

•  com? 

•  com4 

K.  SUBROUTINE  FIND  EX 

1.  Description 

The  incorporation  of  the  effects  of  surface  roughness  into  the  mathematical 
model  results  in  discontinuity  of  the  modal  function  along  the  imaginary  <•/,  ,- axis  .  The 
Shcllman-Morfitt  routine  will  Tail  if  the  edge  of  the  search  rectangle  lies  on  the  imaginary 
axis. 

Subroutine  findfx  will  determine  if  any  edge  of  the  search  rectangle  is  on  the 
imaginary  axis.  If  so,  a  small  offset  of  1  x  10 •'*  is  introduced.  If  the  left  edge  is  on  the 
imaginary  axis,  then  the  modal  function  is  evaluated  at  Re(qu)  =  +  offset.  If  the  right 
edge  is  on  the  imaginary  axis,  the  modal  function  is  evaluated  at  Re(q , ,)  =  —  offset  . 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  findfx  (jr,  ji.  f.  fe.  tleft.  tright). 

3.  Inputs 

The  input  variables  arc 

•  jr  :  Real  part  of  the  mesh  square  coordinates  in  mesh  units. 

•  ji  :  Imaginary  part  of  the  mesh  squares  coordinates  in  mesh  units. 

•  tleft  :  Real  part  of  qu  at  the  left  edge  of  the  search  rectangle. 

•  tright  :  Real  part  of  </, ,  at  the  right  edge  of  the  search  rectangle. 
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4.  Outputs 

The  output  variables  are 

•  f  :  Complex  amplitude  of  the  modal  function. 

•  fe  :  Real  exponent  of  the  modal  Junction. 

5.  Calling  program  element 

•  fzerox 

6.  Subroutine  called 

•  fctvlx 

7.  Common  block  area 

•  Imccom 


L.  SI  DROUT1NE  FND.MOD 
1.  Description 

This  subroutine  sets  up  the  rectangular  region  in  the  complex  t/;  ,-plane  for  the 
Shellman-Moriiu  root  finding  routine  to  locate  the  zeros  of  the  modal  function. 

The  initial  search  rectangle  is  set  up  with 


•  the  left  edge  of  rectangle,  tlclt  =  <>.<». 

•  the  right  edge  of  rectangle,  tright  =  tleft  +  tstep.  where  tstep  is  the  length  of  the 
rectangle. 


the  top  of  the  rectangle,  ttop  given  by  |Ref.  5] 


ttop  = 


2  x  in  ?  A  loss 
A2o  log  c 


(3-40; 


where  Aloss  is  the  maximum  attenuation  rate  in  decibels  per  kilometer  of  those 
modes  to  be  found. 

•  the  bottom  of  the  rectangle,  tbot  =  —  tol.  where  tol  is  the  tolerance  to  which  zeros 
are  to  be  located. 


In  order  that  the  zeros  on  the  real  axis  can  be  located,  a  small  offset  of  qshift 
is  introduced  in  subroutines  fdfdtx  and  fctvlx  so  as  to  avoid  having  a  mesh  line  on  the 
real  axis.  This  offset  is  corrected  in  the  subroutine  fndmod  by  the  statement 


zeros(A)  =  zeros(A)  —  j  qshifi 


(341) 


After  the  initial  rectangle  has  been  searched  for  zeros  of  the  modal  function,  a 
new  search  rectangle  is  set  up  to  the  left  of  the  initial  rectangle.  If  no  zero  is  found  in 
three  consecutive  search  rectangles,  the  search  for  roots  to  the  left  of  the  initial  rectangle 
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is  stopped.  Fndmod  next  forms  search  rectangles  to  the  right  of  the  initial  rectangles. 
If  no  zero  is  found  in  three  consecutive  search  rectangles,  the  search  for  root'-  to  the 
right  of  the  initial  rectangle  is  stopped. 

Since  subroutine  fzerox  extends  the  search  rectangle  by  one  mesh  unit  on  all 
sides,  fndmod  also  checks  for  and  eliminates  all  zeros  outside  the  current  search 
rectangle. 

After  all  the  zeros  have  been  found,  fndmod  sorts  them  in  order  of  increasing 
real  parts  and  stores  them  in  the  array  called  zeros. 

There  are  two  versions  of  fndmod.  namely,  fndmod20.for  and  fndmodSO.for. 
FndmodSO  sets  up  a  smaller  mesh  and  takes  longer  time  to  find  the  modes.  I'ndmodSO 
should  be  used  only  if  fndmod20  fails  to  find  all  modes. 

2.  Calling  subroutine 

The  subroutine  is  called  by  the  statement 
Call  fndmod  (qeigen.  nrmode.  dmdz,  detadx,  zim.  nzlayr,  aloss.  waveno). 

3.  Inputs 

The  input  variables  are 


•  dmdz 

•  detadz 

•  zim 

•  nzlayr 

•  aloss 

•  waveno 


Arrav  containine  values  of 


d\t 

ik 

</>l 


Arrav  containine  values  of——  . 

</_- 

Array  containing  values  of  M(z)  at  layer  boundaries. 

Number  of  linear  segments  used  to  model  the  modified  refractivity 
profile. 


:  Maximum  attenuation  rate  (in  decibels  per  kilometer)  of  modes  to 
be  found. 


:  Wave  number. 


4.  Outputs 

The  output  variables  are 

•  qeigen  :  Array  containing  the  zeros  of  the  modal  function. 


•  nrmode  :  Actual  number  of  modes  found. 


5.  Calling  program  element 

•  main 

6.  Subroutine  called 

•  fzerox 

7.  Common  block  areas 

•  conv4 

•  coinx 

M.  SUBROUTINE  FZEROX 

1.  Description 

Subroutine  fzerox  is  tlic  main  subroutine  used  for  finding  the  zeros  of  a  com¬ 
plex  function  using  the  Shellman-Morlitt  routine  as  described  in  Chapter  111.  Fzerox 
will  convert  the  coordinates  of  the  search  rectangle  edges  into  mesh  units  and  proceed 
with  the  search  for  sign  changes  of  lm(f)  =  0  along  each  edge.  '1  he  maximum  number  of 
crossings  o(Tm(f)  =  (>  allowed  with  any  edge  ol  the  search  rectangle  is  loti. 

Fzerox  allows  a  maximum  number,  maxnsq.  of  mesh  squares  to  lie  on  any  one 
phase  curve.  Jm(f)  =  0.  If  this  number  is  exceeded,  the  program  will  reduce  the  mesh  size 
by  one-half  and  start  over.  If  the  problem  persists  after  the  mesh  size  has  been  reduced 
by  maxnt  times,  an  error  message  will  be  written  to  the  output  file  and  the  program 
stops.  In  the  program,  nia.xnq  is  set  at  four  times  the  number  of  mesh  squares  along  the 
longest  side  of  the  search  rectangle. 

Fzerox  also  checks  for  modes  found  on  the  same  phase  line.  If  so.  the  program 
will  reduce  the  mesh  size  and  start  over.  If  the  problem  remains  after  mesh  size  has  been 
reduced  by  times,  an  error  message  will  be  written  to  output  file  and  the  program 
will  stop. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  Fzerox  (licit,  tright.  tbot.  ttop,  tmsho,  tol.  mprint.  zeros). 

3.  Inputs 

The  input  variables  are 

•  tlcft  :  Value  of  the  real  part  of  z  at  th^  left  edge  of  search  rectangle. 

•  tright  :  Value  of  the  real  part  of  z  at  the  right  edge  of  search  rectangle. 

•  tbot  :  Value  of  the  imaginary  part  of  z  at  the  bottom  edge  of  the  search 

rectangle. 


•  tol  :  Tolerance  to  which  zero'  arc  to  be  found.  Zeros  located  closer  than 

tol  cannot  be  distinguished. 

•  niprint  :  A  flag  for  debugging  output.  If  mprint  =  0,  no  debugging  will  be 

given.  M print  =  1  will  activate  debugging  printout. 

4.  Outputs 

The  output  variables  are 

•  zeros  :  Array  containing  the  complex  zeros  of  f(7-)  in  the  specified 

search  rectangle  in  the  complex  z-plane. 

•  nrz  :  The  actual  number  of  complex  zeros  of  f(z)  found. 

5.  Calling  program  element 

•  fndmod 

6.  Subroutines  called 

•  find  lx 

•  quad 

•  nomshx 

7.  Common  block  areas 

•  newmsh 

•  tmccom 

N.  SUBROUTINE  HTGAIN 

1.  Description 

This  subroutine  evaluates  the  normalized  height  gain  function  given  by 
Equation  77. 

2.  Calling  Statement 

The  subroutine  is  called  by  the  statement 

Call  htgain  (htga.  htge.  zero,  z,  m). 

3.  Inputs 

The  input  variables  are 

•  zero  :  Modal  eigenvalue  in  qt  ,-space  for  mode  m. 

•  z  :  1  lcight  (in  meters)  above  ground  at  which  height  gain  is  to  be 

evaluated. 
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•  111 


:  Mode  index. 


4.  Outputs 

The  output  variables  are 

•  htga  :  Complex  height  gain  factor  in  extended  arithmetic. 

•  luge  :  Real  height  gain  exponent  in  extended  arithmetic. 

5.  Calling  program  element 

•  modsum 

6.  Subroutines  called 

•  xedai 

•  xcadd 

7.  Common  block  areas 

•  com: 

•  papl 

•  pap2 

O.  SUBROUTINE  MODSUM 

1.  Description 

This  subroutine  calculates  the  field  strength  relative  to  free  space  given  by 
Equations  U>0  and  KM. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  modsum  (corns,  eims.  xm,  rngfac,  zr.  zt.  qcigen.  ncount). 

3.  Inputs 

The  input  variables  are 

•  xm  :  Range  in  meters. 

•  rngfac  :  Range  factor  calculated  in  the  main  program. 

•  zr  :  Receiver  height  in  meters. 

•  zt  :  Transmitter  height  in  meters. 

•  qeigen  :  Array  containing  eigenvalues  in  q,  ,-spacc  . 

•  ncount  :  Counter  set  in  the  main  program  to  avoid  unnecessary  calulations  of 

acoefa.  acoefe.  bcoefa  and  bcoefe.  If  ncount  =  1.  subroutine 
modsum  will  call  subroutine  abcoef.  otherwise  subroutine 
modsum  will  compute  modesum  only. 
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4.  Outputs 

The  output  variables  are 

•  eems  :  Coherent  field  strength  relative  to  free  space. 

•  cims  :  Incoherent  field  strength  relative  to  free  space. 

5.  Calling  program  element 

•  main 

6.  Subroutines  called 

•  abeoef 

•  htgain 

•  xcadd 

•  addx 

7.  Common  block  areas 

•  coml 

•  com2 

•  pap2 

P.  SUBROUTINE  NOMSUX 

1.  Description 

Subroutine  nomshx  takes  the  approximate  locations  of  the  complex  zeros  of 
f<[z)  provided  by  fzerox  and  improves  the  accuracies  of  the  zero  locations  using  New- 
ton-Raphson  iteration. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement  . 

Call  nomshx  (tol.  zeros,  nrz). 

3.  Inputs 

The  input  variables  are 

•  tol  :  The  tolerance  to  which  zeros  are  to  be  located.  Zeros  located  closer 

than  tol  cannot  be  distinguished. 

•  zeros  :  Array  containing  the  approximate  locations  of  the  complex  zeros  of 
ffZ). 

:  The  number  of  complex  zeros  stored  in  array  zeros. 


•  nrz 


4.  Output 

‘I  he  output  variable  is 

•  zeros  :  An  array  containing  the  iterated  values  of  the  complex  zeros  of  fly). 

5.  Calling  program  element 

•  fzerox 

6.  Subroutine  called 

•  fdfdtx 

7.  Common  block  area 

•  Newmsh 

Q.  SUBROUTINE  NORME 

1.  Description 

1'his  subroutine  "normalizes"  the  complex  extended  numbers,  z  =  za  exp(zc)  such  that 
exp  (  —l.O.i  <  max  [  |  Rc(:u)  |  ,  |  lm(:a)  j  ]  <  exp  ( 1.0)  (342) 

and  ze  lias  integer  values. 

2.  Calling  statement 

The  program  is  called  by  the  statement 

Call  norme  (za.  ze). 

3.  Inputs 

The  input  variables  are 

•  za  :  Complex  amplitude  of  the  extended  complex  number,  z. 

•  ze  :  Real  exponent  of  the  extended  complex  number,  z 

4.  Outputs 

The  output  variables  are 

•  za  :  Complex  amplitude  of  the  "normalized"  extended  complex  number, 

z. 

•  ze  :  Real  exponent  of  the  "normalized"  extended  complex  number,  z. 

5.  Calling  program  elements 

•  xcadd 

•  xedait 

•  xedaig 
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6.  Subroutine  called 


•  none 

7.  Common  block  area 

•  none 

R.  SUBROl  TINE  NORMRE 

l.  Description 

This  subroutine  "normalizes"  the  real  extended  number.  7  =  za  exp(zc).  such 

that 


exp  (-l.(i)  <  |  :a  j  <  exp  (+1.0) 


[3431 


and  ze  has  integer  values. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  normre  (/.a.  ze). 

3.  Inputs 

The  input  variables  are 

•  za  :  Real  amplitude  of  the  real  extended  number,  z. 

•  ze  :  Real  exponent  of  the  real  extended  number,  z. 

4.  Outputs 

The  output  variables  are 

•  za  :  Real  amplitude  of  the  "normalized"  real  extended  number,  z. 

•  ze  :  Real  exponent  of  the  "normalized"  real  extended  number,  z. 

5.  Calling  program  element 

•  addx 

6.  Subroutine  called 

•  none 

7.  Common  block  area 

•  none 
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S.  SUBROUTINE  QUAD 

1.  Description 

This  subroutine  finds  the  roots  of  the  quadiatie  equation  of  the  form  given  by 
Equation  2S3.  If  |  e  |  given  in  Equation  2S5  is  less  than  U.3.  the  roots  are  computed  with 
Equations  295  and  296.  Otherwise  liquations  286  and  2S7  arc  used.  In  addition,  if 
|  j  <  If)-3"  .  subroutine  quad  returns  only  1  solution  given  by  Equation  295. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  quad  (a.  b.  c.  sol.  nrsol.  niprint). 

3.  Inputs 

The  input  variables  are 

•  a  :  The  coeilicient  of  .v:  in  the  quadratic  equation. 

•  b  :  The  coeilicient  of  2\  in  the  quadratic  equation. 

•  c  :  The  constant  term  in  the  quadratic  equation. 

4.  Outputs 

The  output  variables  arc 

•  sol  :  An  array  containing  the  real  roots  of  the  quadratic  equation. 

•  nrsol  :  The  number  of  real  roots  found. 

•  niprint  :  A  flag  for  debugging.  If  mprint  is  not  equal  to  ?cro.  values  of  the 

coefficients  a.  b.  and  c  will  be  written  to  file. out. 

5.  Calling  program  element 

•  f/ero.\ 

6.  Subroutine  called 

•  none 

7.  Common  block  area 

•  none 

T.  SUBROUTINE  SEAH20 
I.  Description 

This  subroutine  evaluates  the  real  effective  relative  dielectric  constant  and  the 
real  effective  conductivity  of  seawater  as  a  function  of  temperature,  salinity  and  fre¬ 
quency.  The  program  is  based  on  the  equations  presented  in  section  G  of  Chapter  II. 
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2.  Calling  statement 

The  subroutine  is  called  b\  the  statement 

Call  seah2o  ( sigelh  epself.  t.  s.  freq). 

3.  Inputs 

The  input  variables  are 

•  t  :  Temperature  of'  seawater  in  degree  Celcius. 

•  s  :  Salinity  of  seawater  in  grams  of  salt  per  kilogram  of  seawater. 

•  freq  :  Frequency  in  Hertz. 

4.  Outputs 

The  output  variables  are 

•  sigelf  :  Real  effective  conductivity  in  siemens  per  meter. 

•  epself  :  Real  effective  relative  dielectric  constant. 

5.  Calling  program  element 

•  main 

6.  Subroutine  called 

•  none 

7.  Common  block  area 

•  none 


U.  SUBROUTINE  SURF 

I.  Description 

This  subroutine  evaluates 

model. 


and  —7- —  when  surface  rouehness  is  included  in  the 


2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  surf  (ql  1,  gamma,  dgamdq). 

3.  Input 

The  input  variable  is 

•  ql  1  :  Location  in  the  complex  <jru-planc  where  the  modal  function  and  its 

derivative  arc  to  be  evaluated. 


4.  Outputs 

The  output  variables  arc 

•  gamma  :  This  is  y,  given  by  Equation  ISO  and  192  for  horizontal  and  vertical 
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•  dgamdq 


polarization  respectively. 


:  This 


given  by  liquations  1ST  and  193  lor  horizontal  and 


vertical  polarization  respectively. 


5.  Calling  program  elements 

•  fctvlx 

•  iiildtx 

6.  Subroutines  called 

•  function  ctaiih 

•  function  csesh2 


7.  C  ommon  block  areas 

•  com! 

•  conk 

•  comb 


V.  SUBROUTINE  XAINEG 

1.  Description 

This  subroutine  evaluates  the  Airy  function.  Ai(z).  and  its  derivative,  .-w'tzi.  for 
complex  z  in  extended  complex  arithmetic  for 


<  are  (r )  <  0 


C'4-ll 


1  he  subroutine  will  test  if  z*  is  within  the  ellipse  with  foci  (  — 15S.  3.95 >  and 
(0.28.  —2.1 1 ).  major  axis  of  9.0  and  minor  axis  of  6.4.  If  r  is  within  the  ellipse,  xaineg 
will  invoke  xedait  to  compute  the  Airy  function  and  its  derivative  with  a  Taylor  series. 
Otherwise  xaineq  invokes  xedaig  which  uses  the  Gaussian  quadrature  method. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  xaineg  (z.  aia.  aie.  daia.  daic). 


3.  Input 

The  input  variable  >s 


A  value  in  the  complex  7-plane  at  which  the  Airy  function  and  its 
derivative  are  to  he  computed. 


IM3 


•  z 


4.  Outputs 


•  aia  :  Complex  amplitude  of'  Ait/)  in  extended  complex  arithmetic. 

•  aie  :  Real  exponent  of  Ai(/ )  in  extended  complex  arithmetic. 

•  daia  :  Complex  amplitude  of  Ai  '  (/)  in  extended  complex  arithmetic. 

•  date  :  Real  exponent  of  Ai  '  (/)  in  extended  complex  arithmetic. 

5.  Calling  program  element 

•  xedai 

6.  Subroutines  called 

•  xedait 

•  xcdaie 

7.  Common  block  area 

•  none 

\\.  Sl’BROl  TINE  XAITOS 
1.  Description 

This  subroutine  evaluates  the  Airy  function.  Ah / 1.  and  its  derivative.  Ai  t/ 1  in 
extended  complex  arithmetic  for 

0  <  are  (:)  <  -rr~  .  (345) 


The  subroutine  will  test  if  /  is  within  the  ellipse  with  loci  (-1.58.  3.95)  and 
(0.28.  —2.11),  major  axis  of  9.0  and  minor  axis  of  6.4.  If  /  is  within  ellipse,  xaipos  will 
invoke  xedait  to  compute  the  Ain  function  and  its  derivative  with  a  Tax  lor  series.  Oth¬ 
erwise  xaipos  invokes  xedaig  which  uses  the  Gaussian  quadrature  method. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 


3. 


Call  xaipos  (/.  aia.  aie.  daia.  daic). 

Input 

The  input  variable  is 


•  /  :  A  value  in  the  complex  /.-plane  at  which  the  Airy  function  and  its 

derivative  are  to  be  computed. 


•4.  Outputs 

I  he  output  variables  are 

•  aia  :  Complex  amplitude  of  Aim  in  extended  complex  arithmetic. 

•  aie  :  Real  exponent  of  Ait/i  in  extended  complex  arithmetic. 

•  daia  :  Complex  amplitude  of  Aim  in  extended  complex  arithmetic. 

•  date  :  Real  exponent  of  Ai  (/)  in  extended  complex  arithmetic. 

5.  Calling  program  element 

•  xedai 

6.  Subroutines  called 

•  xedait 

•  xcdaie 

7.  Common  block  area 

•  none 

SUBROUTINE  XCDAI 
1.  Description 

The  subroutine  is  the  driving  program  in  the  evaluation  of  the  Air>  function. 
Aim.  and  its  derivative.  Ai'( / ).  If 


0  <  atg(r)  <  —  .  (34<d 

xedai  invokes  xaipos  to  compute  Ait/)  and  Aim  and  invokes  xuineg  to  compute 
.hire  -  ■'  •')  and  . f/'tre  ~ ')  .  It  then  uses  the  connection  formula  of  Airy  function  in 
Equation  3<>2  and  its  derivative  to  compute  dure  J)  and  d/’fre  d  . 

If 


--ip  <  are(r)  <  0  ,  (347) 

xedai  invokes  xaipos  to  compute  Aiizc  ■1"  ’)  and  d/'(re  ■ )  .  It  also  invokes  xainee  to 
compute  Ait/)  and  Ai  t/).  It  then  uses  the  connection  formula  to  compute  .f/tre  _  •') 

and  d/’(re  *’"  ’)  . 

If 


~~r~  <  arc  (r  l  < 


1 34  Si 


1<>5 


or 


<  are  |e) 


xcdai  calls  xaipos  to  compute  A /(re  ~  ;t  ’’)  and  . I /'(rf  -1,?)  .  It  also  calls  xaincc  to  com¬ 
pute  . //( - c  -'’)  and  Ai'(:e  ■')  .  It  then  uses  the  connection  formula  to  compute  Aip/j 
and  Ai  (/>- 

If 


-V  <  arc  (r)  <  .  (35<0 

xcdai  calls  xaipos  to  compute  Ai(z)  and  Ai’(z)  and  calls  xaincc  to  compute  AH:e  ') 
and  .f.'lce  ‘ J)  .  It  then  uses  the  connection  formula  to  conpute  Ailze  *  •  ?l  and 
J/'U-c  ■’-)  . 

If 

<  arc  (r)  <  ——  .  (351) 

a  2 

xcdai  calls  xaipos  to  compute  Ai{ze  5)  and  Ai’(:c  ■’)  .  It  also  calls  xaincc  to  compute 
Ail./)  and  Ai(z).  It  then  uses  the  connection  formula  to  compute  .hire--*;)  and 
Ai'lzc  -■**  ')  . 

2.  Calling  statement 

1  he  subroutine  is  called  by  the  statement 

Call  xcdai  (z,  aia.  aie,  daiu.  daic.  aipa.  aipe.  daipa.  daipe.  aima. 
aime.  daima.  daime). 

3.  Input 

The  input  variable  is 

•  z  :  A  value  in  the  complex  z-planc  at  which  the  Airy  function  and  it.. 

derivative  arc  to  be  computed. 

4.  Outputs 

The  output  variables  are 

•  aia  :  Complex  amplitude  of  Ai{z)  in  extended  arithmetic. 
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•  aic 


:  Real  exponent  of  Ai{z)  in  extended  complex  arithmetic. 

•  daia  :  Complex  amplitude  of  Ai'(z)  in  extended  complex  arithmetic. 

•  daie  :  Real  exponent  of  Ai'(z)  in  extended  complex  arithmetic. 

•  aipa  :  Complex  amplitude  of  in  extended  complex  arithmetic. 

•  aipc  :  Real  exponent  of  Ai{:e ,2r  i)  in  extended  complex  arithmetic. 

•  daipa  :  Complex  amplitude  of  Ai’(zcJ3n})  in  extended  complex  arithmetic. 

•  daipe  :  Real  exponent  of  Ai'ize^3)  in  extended  complex  arithmetic. 

•  aima  :  Complex  amplitude  of  Ai(ze  in  extended  complex  arithmetic. 

•  aime  :  Real  exponent  of  Ai(ze  in  extended  complex  arithmetic. 

•  daima  :  Complex  amplitude  of  Ai‘{ze  -J3**)  in  extended  complex  arithmetic. 

•  daime  :  Real  exponent  of  Ai\:e  in  extended  complex  arithmetic. 

5.  Calling  program  element 

•  fctvlx 

6.  Subroutines  called 

•  xaipos 

•  xaineg 

•  xcadd 

7.  Common  block  area 

•  none 

Y.  SUBROUTINE  XCDAIG  !. 

1.  Description 

The  subroutine  computes  the  Airy  function,  Ai(z),  and  its  derivative,  Ai'(z),  in 
extended  arithmetic  using  the  Gaussian  quadrature  method. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  xedaig  (z,  aia,  aic,  daia,  daie). 

3.  Inputs 

The  input  variables  is 

•  z  :  A  value  in  the  complex  z-planc  at  which  the  Airy  function  and  its 

derivative  arc  to  be  computed. 
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4.  Outputs 

The  output  variables  are 

•  aia  :  Complex  amplitude  of  Ai(z)  in  extended  arithmetic. 

•  aie  :  Real  exponent  of  Ai(z)  in  extended  complex  arithmetic. 

•  daia  :  Complex  amplitude  of  Ai'(z)  in  extended  complex  arithmetic. 

•  daie  :  Real  exponent  of  Ai'(z)  in  extended  complex  arithmetic. 

5.  Calling  program  elements 

•  xaipos 

•  xaineg 

6.  Subroutine  called 

•  norme 

7.  Common  block  area 

•  norme 

Z.  SUBROUTINE  XCDAIT 

1.  Description 

The  subroutine  computes  the  Airy  function.  Ai(z).  and  its  derivative.  Ai'(z),  in 
extended  arithmetic  using  a  Taylor  series  expansion. 

2.  Calling  statement 

The  subroutine  is  called  by  the  statement 

Call  xedait  (z.  aia.  aie.  daia.  daie). 

3.  Input 

The  input  variable  is 

•  z  :  A  value  in  the  complex  z-planc  at  which  the  Airy  function  and  its 

derivative  are  to  be  computed. 

4.  Outputs 

The  output  variables  are 

•  aia  :  Complex  amplitude  of  Ai(z)  in  extended  arithmetic. 

•  aie  :  Real  exponent  of  Ai(z)  in  extended  complex  arithmetic. 

•  daia  :  Complex  amplitude  of  Ai'(z)  in  extended  complex  arithmetic. 

•  daie  :  Real  exponent  of  Ai  (z)  in  extended  complex  arithmetic. 
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5.  Calling  program  elements 

•  xaipos 

•  xaineg 

6.  Subroutine  called 

•  normc 

7.  Common  block  area 

•  none 

AA.  FUNCTION  PROGRAMS 

1.  Function  casin 

Function  casin  evaluates  the  complex  arcsine  of  a  complex  number.  /.  with  the 

equation 

sin-1  (r)  =  -j  In  [( 1  -:2yn  +  jz]  .  (352) 

The  function  is  called  by  the  statement 

casin  (z). 

The  function  is  used  in  the  main  program  to  convert  the  eigenvalues  from 
</)_, -plane  to  0-plane  . 

2.  Function  csech2 

Function  cscch2  evaluates  the  square  of  the  hyperbolic  secant  of  a  complex 
number,  z.  with  the  equation 

scch2  (r)  =  — - - - -  •  (353) 

c“  +  2  +  e  ~* 

The  function  is  called  by  the  statement 

csech2  (z). 

It  is  used  in  subroutine  surf  for  evaluating  ■  1 a—  ,  and  ~— 

Ct/l.l  (  Cli l.l  C  £/ 1  | 

(see  Section  F,  Chapter  II). 

3.  Function  ctanh 

Function  ctanh  evaluates  the  hyperbolic  tangent  of  a  complex  number,  z.  with 
the  equation 
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tanh  (r)  = 


c*  —  t‘ 

for  |  Rc(:)  ]  >  - Ur 

10  ,  2 

or  |  Iut(z)  |  >  - —=r 

10  v  2 


(354; 


or 


tanh  (r) 


1  -3  _2_  .5  _  _P_  ,7  .9 

3  ‘  15  315  ‘  2835  ‘ 


for  |  Rc(r)  |  <  - 

10,2 


or  |  //;;(r)  |  <  - = 

10  N  2 


(355) 


The  function  is  called  by  the  statement 

ctanh  (z). 

It  is  used  in  subroutine  surf  for  evaluating  a.  b,  a  ,b  and  their  derivatives. 

4.  Function  ibstrip 

Function  ibstrip  finds  the  length  of  a  character  string  by  removing  all  trailing 
spaces  from  the  end  of  the  string  and  returns  the  number  of  characters  in  the  stripped 
string. 

The  function  is  called  by  the  statement 

ibstrip  tinst,  outst). 

where  inst  is  the  input  string  and  outst  is  the  stripped  string. 

The  function  is  used  in  the  main  program  to  check  the  length  of  the  filename. 

AB.  MLAYER  SUPPORTING  PROGRAMS 

M  LAYER  has  several  supporting  programs  that  helps  in  the  running  of  the  pro¬ 
gram.  They  are 

•  mlayer.mak 

•  vwgclean.mak 

•  wvgstrip.mak 

•  wvgclean.c 
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•  wvgstrip.for 

•  wvestdin.for 


1.  Mlayer.niak 

This  is  the  makefile  for  the  creation  of  mla2<)  and  mlaSO  programs.  Mla20  uses 
subroutine  fndmoddO.for  to  set  up  the  rectangular  search  region  whereas  mlaSO  uses 
fndmodSo.for. 

Mlayer.niak  compiles  and  links  all  program  elements  and  creates  mla20.exe  and 
mlaS0.exe  for  the  execution  of  Ml.AYLR  program  with  either  mla20  or  mlaSO. 

To  compile  and  link  MLAYliR.  type 

nmake  mlayer.niak 

Nmake  will  compare  the  modification  dates  of  the  target  files  (e.g..  object  and  executable 
files)  with  those  of  the  dependent  files  (e.g..  source  files).  If  any  of  the  dependent  files 
has  been  changed  recently,  i.e..  the  modification  dates  of  the  dependent  files  are  more 
recent  than  the  target  files,  nmake  will  execute  the  commands  in  mlayer.niak  to  update 
the  target  files.  Nmake  updates  only  those  outdated  target  files. 

To  execute  M LAYER  program,  type 

mlado  <  infitle 

to  execute  the  ml  a  20  version.  Infile  is  the  name  of  the  input  file  which  can  be  either  filein 
or  filein. eig.  If  filein  is  used,  MLAYER  will  search  for  all  modes  of  propagation.  If  fi¬ 
lein. eig  is  used.  MLAYER  will  compute  the  rnodsum  and  path  loss  without  a  search  for 
the  modes. 

2.  Wvgclean.mak 

This  is  the  makefile  for  the  creation  ofwvgclcan.exe 

3.  Wvgstrip.niak 

This  is  the  makefile  for  the  creation  of\wgstrip.exe. 

4.  Wvgclean.c 

Wvgclean.c  adds  comment  or  descriptors  to  each  datum  line  in  filein  and  fi¬ 
lein. eig.  This  will  help  in  editing. 

To  execute  the  program,  type 

wvgelean  <  filein. eig  >  temp 
ren  temp  filein. eig 

5.  Wvgstrip.for 

Wvgstrip.for  removes  comments  or  descriptors  that  arc  added  by  wvgclean.c. 

To  execute  the  program,  type 
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temp 


6. 


wvgstrip  <  filein.eig  > 
ion  temp  filein.eig 

Wvgstdin.for 

Wvgstdin.for  reads  input  data  from  lilein  or  filein.cig  on  execution  of  the 


M LAYER  program. 


112 


VI.  DISCUSSIONS  AND  RECOMMENDATIONS 

A.  ASSESSMENT  OF  MLAYER 

MLAYliR  was  developed  with  the  intention  of  providing  a  program  which  is  ca¬ 
pable  oflocating  all  propagating  modes  with  attenuation  rates  below  a  specified  value. 
However,  the  price  of  this  capability  is  the  extremely  long  execution  time.  For  instance, 
a  2-meter  evaporation  duct  takes  approximately  3  hours  and  2<)  minutes  (on  an  IBM 
PS  2  model  SO  with,  an  Intel  S03S6  processor  at  16  Mil/)  to  locate  9  modes  with  atten¬ 
uation  rates  less  than  5  decibels  per  kilometer.  A  14-meter  evaporation  duct  takes  about 
6  hours  to  locate  94  modes  with  attenuation  rates  below  2.1  decibels  per  kilometer. 
Samples  of  output  data  arc  found  in  Appendix  B. 

Calculations  of  field  strength  and  path  loss  with  MLAYER  were  found  to  agree 
favourably  with  experimental  measurements  |Rcf.  17  and  IS].  Typically,  the  model  un¬ 
derestimates  the  measured  values  by  about  10  dB.  The  sources  of  discrepancies  were 
probably  due  to  the  validity  of  the  surface  roughness  model  and  the  assumption  of  lat¬ 
eral  homogeneity  of  the  medium. 

In  the  following  sections,  proposed  areas  for  enhancement  of  the  capabilities  of 
MLAYER  tire  discussed. 

B.  SURFACE  ROUGHNESS  MODEL 

The  surface  roughness  model  used  in  MLAYER  is  based  on  the  KirchhoiT  approx¬ 
imation  which  is  a  single  scattering  theory.  This  model  is  only  valid  for  gently  undulating 
surfaces  where  shadowing  of  incident  and  scattered  fields  are  negligible.  When  shadow¬ 
ing  cannot  be  neglected,  the  KirchhoiT  approximation  tends  to  ovcrprcdict  the  scattered 
energy  which  leads  to  erroneous  results. 

Shadowing  occurs  when  incoming  or  outgoing  rays  are  blocked  by  the  rough  sur¬ 
face.  This  leaves  some  area  of  the  surface  in  shadow.  This  effect  is  associated  with  mul¬ 
tiple  scattering.  Therefore,  to  account  for  shadowing,  multiple  scattering  interactions 
must  be  included  in  the  model. 

A  more  accurate  model  called  the  phase  perturbation  technique  was  introduced  by 
Winebrcnner  and  lshimaru  [Ref.  19].  The  phase  perturbation  theory  incorporates  all  the 
multiple  scatterings  required  for  shadowing  of  the  incident  field.  The  theory  is  based  on 
the  extinction  theorem  which  is  a  consequence  of  Green's  theorem  [Ref.  20].  In  Refer¬ 
ence  19.  numerical  results  of  the  phase  perturbation  bistatic  scattering  cross  section  were 
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compared  with  exact  numerical  results.  In  the  region  examined,  it  was  found  that  the 
results  of  the  phase  perturbation  technique  agree  with  the  exact  results  except  at  low 
grazing  angles. 

The  phase  perturbation  technique  is  valid  for  a  wider  range  of  applications.  It  was 
shown  in  Reference  19  that  the  rcllection  and  backscattering  coefficients  reduce  to  that 
of  the  KircholT  approximation  when  Ij).  >  1,  i.c..  for  gently  undulating  surfaces,  where 
/  is  the  length  and  /.  is  the  wavelength  of  the  incident  wave. 

C.  LATERAL  HOMOGENEITY  OF  REFRACTIVITY  PROFILE 

In  M LAY  ER,  the  refractive  index  is  assumed  to  vary  only  with  height.  Although 
this  assumption  appears  to  be  adequate  most  of  the  time  [Ref.  17].  there  are  incidents 
where  the  validity  of  lateral  homogeneity  was  cited  as  a  possible  reason  l  or  discrepancies 
between  observations  and  predicted  results  [Ref.  IS  and  21]. 

For  refractivity  profile  that  varies  along  the  path,  the  laterally  inhomogeneous 
structures  can  be  approximated  by  several  homogeneous  sections.  The  exchange  of  en¬ 
ergy  between  the  propagating  modes  at  the  boundary  of  two  sections  are  analyzed  by 
mode-conversion  techniques  [Ref.  22]. 

D.  FIELD  CONTRIBUTION  FROM  THE  SOURCE 

In  the  formulation  of  the  propagation  model,  M LAY  ER  neglects  the  contribution 
from  the  direct  wave  of  the  source.  This  approximation  is  good  for  ranges  near  or  be¬ 
yond  the  horizon  where  the  major  field  contributions  arc  from  the  rcllected  waves.  For 
ranges  well  within  the  horizon,  the  contribution  from  the  direct  wave  is  significant  and 
should  not  be  ignored.  Hence,  for  wave  propagation  in  the  interference  region,  the  co¬ 
efficients,  A,  and  B,,  of  Equations  31  and  32  must  be  obtained  from  the  solutions  of 
Equation  46.  At  the  same  time,  the  surface  wave  field  contribution  from  the  integral  of 
the  contour.  C2  ,  in  Equation  73  may  become  significant. 


114 


MI.  CONCLUSION 


MLAYLR  is  a  useful  research  tool  for  conducting  case  studies  over  a  large  dynamic 
range  of  frequencies.  Calculation  of  field  strength  and  path  loss  arc  in  agreement  with 
observations  for  frequencies  from  0?  Mhz  to  94  GHz  (Ref.  17  and  IS).  It  can  serve  as  a 
yardstick  against  which  results  of  other  quicker  but  less  accurate  models  such  as  IRITS 
[Ref.  27]  and  LRU’S  [Ref.  24]  are  compared.  In  addition,  the  multiple  (as  opposed  to 
trilincan  piecewise  linear  refractivity  profile  makes  it  possible  to  model  the  simultaneous 
occurrence  of  an  evaporation  duct  and  elevated  ducts. 

finally,  modifications  of  MLAYLR  to  incorporate  conditions  of  horizontal  heter¬ 
ogeneity  of  refractivity,  effects  of  source  contribution  and  an  improved  me  idling  of 
surface  roughness  as  discussed  in  Chapter  VI  are  recommended  to  fully  exploit  its  ca¬ 
pabilities.  These  are  opportunities  opened  for  further  studies. 
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APPENDIX  A.  REPRESENTATION  OF  A,  (q)  AND  A,  (q)  IN  TERMS  OF 


The  Airy  function  is  related  to  the  modified  Bessel  function  of  the  third  kind  of  or¬ 
der  one-third  by  the  following  equation  [Ref.  4): 


Mh)  -  i  (f  db  . 


(356) 


where 


AV's(</)  =jf  ^//,;3(,,(‘/^/J)  (357) 

or 

A',/3  (<])  *  -./-T  e_/' *  H];?a)  {qe~Jr2)  .  (35$) 

In  Equations  357  and  358.  //,  ,(l>(.)  and  //,  ?  c>  (.)  are  llankcl  functions  of  the  first  and 
second  kind  of  order  one-third  respectively.  Substituting  Equation  357  into  Equation  356 
leads  to 


Thus. 


II 


1/3 


(1) 


Ai(-qe^)  =  i(-|-),^),/V2*/3) 


1/2  c  Jr! 6 


I  hi 


(t> 


3/2  —73-/2 

<7  *  ■ 


/ 

■> 


// 


1/3 


(1) 


(359) 


(360) 
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From  Equation  142. 


A ; ( a )  =  —j  2 ( 12]  Ai(  j 


I'M) 


Substituting  Equation  36U  inio  Equation  3M  gives 


=  h\  I'/) 


:  36: 


where  /;,  (</)  is  tiic  modified  Hankel  Function  of  the  first  kind  of  order  one-third  and  is 
defined  as  (Ref.  5] 


/i,  (</  )  = 


t 


To  derive  an  expression  for  the  .Airy  function.  .If  -  </).  can  be  expressed  as 

[Ref.  4) 


// 


1/3 


(!) 


.(364) 


From  Equation  143. 


k-.U])  =  2(12 )v*ejr!*  AH  -</) 


(365) 


Substituting  Equation  364  into  365  leads  to 
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/<;(,/)  =  ( 1 2 ) 1  '  eJ~: 


</ 


+  c 


II 


C) 


1/3 


(366) 


<T 


/J- r/3 

c 


=  /(,  ((/) 


/»,  (</•) 


t 


where  /!,(</)  is  the  modified  Ilankcl  function  of  the  second  kind  of  order  one-third  and 
is  defined  as  |Ref.  5j 


/<:(</)  = 


(36' i 


Tor  large  1 1/ 1  .  the  Ilankcl  function  can  be  approximated  by  the  leading  term  of  its  as¬ 
ymptotic  expansion.  That  is  |Rcf  4 j. 


II 


1/3 


(1) 


l‘/> 


/  :  \ 

1  .'2 

t  ~  * 

\  -</  j 

exp 

V  6  4  )_ 

—  r  <  arg(i/l  <  2~ 


(36S) 


and 


II 


1/3 


/  :  \ 

i  '2  r  .  /  r 

(  -‘/  ) 

cxp[-./  (</-- 

4  )J 

—  2r:  <  arg  (</)  <  ff 


(369) 


1  IS 


I 


1  herefote. 


k  j  ( </ )  -  1: 


V  •  /  1 


(intPt^r- 


^  \i 


cxr  ■'  T  li 


,-»! -n  -i  :  -!  '4  7  J 

*  —  /?  pvn  il  - 


<1  c\p  ./  -rq 


2-  ,  4- 

- —  <  areb/)  <  — — 

J>  J 


To  obtain  the  asymptotic  expansion  for  k:(q)  .  it  is  more  convenient  to  express  in 
terms  of//,  ,  ,:’F)  •  1  rom  liquations  356  and  35S. 

.•i/n/i  =  4-  ( -j  )1;‘  (  -,/y )  <’ ‘ '  u\ '.I  (  y  </"•*  ) 


.in  -</)  =  -y  (y  ),':c"y"  ‘c  '*'v 


3  2  -;'-e 

V  ‘  "  < 


—HI s'- 


Substitute  liquation  372  into  365  gives 


m«>«  -(Ti,n-)'!2"T!'(t<T'-r-) 


From  liquation  369  and  373,  the  asymptotic  expansion  for  /,:  (</)  is 
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APPENDIX  B.  SAMPLE  PRINTOUT  OF  INPUT  AND  OUTPUT  FILES 


**fc**#**s ample  printout  for  a  2  m  evaporation  duct  input  file*' 


9ghz02m 

0 


9600. 00 

FQMZIN 

0. 000000 

DELFQ 

1. 000000 

NFREQ 

0. 000000 

MPOL 

5. 000000 

AL0SS 

15. 000000 

SEATMP 

35. 000000 

SEASLT 

1. 000000 

I FLAG E 

0. 000000 

AIRTMP 

0. 000000 

RH 

0. 000000 

WGPM3 

0. 25uOCC 

RMSBHT 

48. 000000 

2.TIN1T 

0. 000000 

DELZT 

1. 000000 

NZT 

46. 0000000 

ZRINIT 

7. 0000000 

DELZR 

2. 0000000 

NZR 

18. 5000 

XIN'IT 

9. 250000 

DELX 

3. 000000 

NX 

0. 000000 

ZREF 

16 

N'ZLAYR 

0. 000000 

zi[0] 

0. 000000 

zin[ 0] 

0. 000000 

zigab[  0] 

0. 020000 

zi[  1] 

-0. 400000 

zim[  1] 

0. 000000 

zigab[  1] 

0. 040000 

z-:[2] 

-0. 550000 

zim[  2] 

0. 000000 

zigab[  2] 

0. 079000 

zi[  3] 

-0. 700000 

zim[  3] 

0. 000000 

zigab[  3] 

0. 158000 

zi[4] 

-0. 850000 

zim[  4] 

0. 000000 

z.igab[4] 

0. 251000 

zi[5] 

-0. 940000 

zim[ 5] 

0. 000000 

zigabf  5] 

0. 631000 

zi[  6] 

-1. 100000 

zim[  6] 

0. 000000 

zigab[ 6] 

1. 259000 

zi[  7] 
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-1. 170000 

zim[  7] 

0. 000000 

zigabf  7] 

2. 000000 

zi  l  8] 

-1. 180000 

zim[  S] 

0. 000000 

zigab] 8] 

3. 981000 

zi[  9] 

-1. 090000 

zim[  9] 

0. 000000 

zigab[ 9] 

5. 012000 

zi[  10] 

-1. 010000 

zim] 10] 

0. 000000 

zigab] 10] 

7. 943000 

zi[  11] 

-0. 750000 

zim] 11] 

0. 000000 

zigab[ 11] 

12. 589000 

zi[  12] 

-0. 270000 

zim[ 12] 

0. 000000 

zigab[ 12] 

19. 953000 

zi[  13] 

0. 550000 

zim]  13] 

0. 000000 

zigab] 13] 

25. 119000 

zi]  14] 

1. 140000 

zim]  14] 

0. 000000 

zigab] 14] 

31. 623000 

zi] 15] 

1. 900000 

zim] 15] 

0. 000000 

zigab] 15] 

50. 000000 

z i] 16] 

4.  069000 

zim] 16] 

0. 000000 

zigab] 16] 

********  sample  printout  for  9ghz02.  eig******,v* 


9ghz02m 


1 

9600. 000 

FQMZIN 

0. 000000 

DELFQ 

1 

NFREQ 

0 

MPOL 

5. 000000 

ALOSS 

15. 000000 

SEATMP 

35. 000000 

SEASLT 

1 

IFLAGB 

0. 000000 

AIRTMP 

0. 000000 

RH 

0. 000000 

WGPM3 

0. 250000 

RMSBHT 

48. 000000 

ZTINIT 

0. 000000 

DELZT 

1 

NZT 

46. 000000 

ZRINIT 

7. 000000 

DELZR 

2 

NZR 

18. 500000 

XINIT 

9. 250000 

DELX 

3 

NX 
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0. 000000 

ZREF 

16 

NZLAYR 

0. 000000 

zl[  0] 

0.  000000 

zim[ 0| 

0.  000000 

zigabf  0] 

0.  020000 

zif  1] 

-0.  400000 

zimf 1] 

0.  000000 

zigabf  1] 

0.  040000 

zif  2] 

-0.  550000 

zimf  2] 

0.  000000 

zigab[  2] 

0.  079000 

zi[  3] 

-0.  700000 

zim[ 3] 

0.  000000 

zigabf 3] 

0.  156000 

zi(  4] 

-0.  850000 

zitn[  4] 

0.  000000 

zigab[ 4] 

0.  251000 

zi[5] 

-0. 940000 

zim[ 5] 

0. 000000 

zigabf 5] 

0. 631000 

zi[  6] 

-1. 100000 

zim[ 6] 

0. 000000 

zigabf  6] 

1. 259000 

zif  7J 

-1.  170000 

zimf  7] 

0. 000000 

zigab[ 7] 

2. 000000 

zi[  8J 

-1. 180000 

zimf  8] 

0. 000000 

zigabf  8] 

3.  981000 

zi(  9) 

-1. 090000 

zimf 9] 

0. 000000 

zigabf  9] 

5. 012000 

zi[  10] 

-1.  010000 

zimf 10] 

0. 000000 

zigab[  10] 

7. 943000 

zi[  11] 

-0. 750000 

zim[  11] 

0. 000000 

zigab(  11] 

12. 589000 

zi[  12] 

-0. 270000 

zini[  12] 

0. 000000 

zigab[ 12] 

19. 953000 

zi[  13] 

0.  550000 

zim[ 13] 

0. 000000 

zigabf 13] 

25. 119000 

zi[  14] 

1. 140000 

zimf 14] 

0. 000000 

zigabf  14] 

31. 623000 

zif  15] 

1. 900000 

zimf  15] 

0. 000000 

zigabf  15] 

50. 000000 

zif  16] 

4. 069000 

zimf 16] 

0. 000000 

zigabf 16] 

9 

nrmode 

( -1. 2695 569 83588969E-001 , 1.  637613838717031E-001)  [  1] 

( -9. 29 119430069269 IE -002,1.  251243717867209E-001)  [ 2] 
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(-4.  143216670274859E-002,8.  947793427709196E-002)  [ 3] 
(8. 2407061 15401963E -003, 6.  9 17378956 118456E-002)  [ 4] 
(5. 55 1 1 60 395 245 24 4E -002, 5. 3248181 13693566E-002)  [ 5] 
(7. 6826 113485 04449 E- 002 ,6.  708823745503940E-002)  [ 6] 
(1. 394420773218969E-001,9. 725689299213881E-002)  [ 7] 
(2.  102775 54249550 3E -001,1.  357828223763695E-001)  [ 8] 
(2. S645891S3106626E -001,1. 658090146540282E-001)  [ 9] 


VfiHnWc: 


Sample  printout  for  9ghz02m. out******** 


1 


frequency  =  9600. 0000  mhz 

horizontal  polarization 

maximum  mode  attenuation  rate  =  5. 0000  db/km 

seawater  temperature  =  15.0000  degrees  Celsius 

seawater  salinity  =  35. 0000  grams  salt/kg  seawater 

dielectric  constant  of  seawater  =  80. 8869 

conductivity  of  seawater  =  4.6400  si/m 

rms  surface  bump  height  =  . 2500  meters 

iflagb  =  1 

air  temperature  =  . 0000  degrees  Celsius 

relative  himidity  =  .0000  percent 

liquid  water  concentration  in  air  =  . 0000  grams/meter**3 

m(  . 0000  )  =  . 0000  m-units 

gab(  . 0000  )  =  . 0000  db/km 


tropospheric  modified  refractivity  profile 


z 

m(z) 

gab(z) 

dm/dz 

d(gab)/'dz 

meters) 

(m-units) 

(db/km) 

( m-units /meter ) 

( ( db/km)/meter) 

.  0000 

.  0000 

. 00000D+00 

-20. 0000 

. 00000D+00 

.  0200 

-.  4000 

. OOOOOD+OO 

-7.  5000 

. 00000D+00 

.  0400 

-.  5500 

. OOOOOD+OO 

-3.  8462 

. 00000D+00 

.  0790 

-. 7000 

. 00000D+00 

-1.  8987 

. 00000D+00 

.  1580 

-. 8500 

. OOOOOD+OO 

-. 9677 

. 00000D+00 

.  2510 

-.  9400 

. OOOOOD+OO 

-.4211 

.  00000D+00 

.  6310 

-1.  1000 

. OOOOOD+OO 

-. 1115 

. OOOOOD+OO 

1.  2590 

-1.  1700 

. 00000D+00 

-.  0135 

. 00000D+00 

2. 0000 

-1. 1800 

. 00000D+00 

.  0454 

. 00000D+00 

3. 9810 

-1.  0900 

. 00000D+00 

.  0776 

, OOOOOD+OO 

5. 0120 

-1. 0100 

. 00000D+00 

.  0887 

. 00000D+00 

7. 9430 

-. 7500 

. OOOOOD+OO 

.  1033 

. 00000D+00 

12. 5890 

-.  2700 

. OOOOOD+OO 

.  1114 

.  00000D+00 

19. 9530 

.  5500 

. 00000D+00 

.  1142 

. 00000D+00 

25. 1190 

1.  1400 

. OOOOOD+OO 

.  1169 

. 00000D+00 

31. 6230 

1. 9000 

. 00000D+00 

.  1180 

. 00000D+00 

50. 0000 

4. 0690 

. OOOOOD+OO 
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start  search  for  modal  eigenvalues 


tmesh=  1. 67934D-03  tol=  1. 00000D-04 


zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  0.00D+00  tright=  2. 69D-01 

ttop  =  1. 68D-01  tbot  =  -1. 00D-04 


qeigen(  1)  =  8.24071D-03  6. 91736D-02 

qeigen(  2)=  5.55116D-02  5. 32482D-02 

qeigen(  3)  =  2.  10278D-01  1. 35783D-01 

qeigen(  4)=  1. 39442D-01  9. 72569D-0 

qeigen(  5)  =  7.6S261D-02  6. 70882D-0 


zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  -2.69D-01  tright=  0. 00D+00 

ttop  =  1. 68B-01  tbot  =  -1.00D-04 


qeigen( 

1)  = 

-1. 26956D-01 

1. 63761D-01 

qeigen( 

2)  = 

-9. 29 119D-02 

1. 251240-01 

qeigen( 

3)  = 

-4. 14322D-02 

8. 94779D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tlef t=  -5.38D-01  tright=  -2. 69D-01 

ttop  =  1.68D-01  tbot  =  -1.00D-04 


zeros  found  in  expanded  search  rectangle  defined  by 
t lef t=  -8.06D-01  tright=  -5. 38D-01 

ttop  =  1. 6SD-01  tbot  =  -1.00D-04 


zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  2.69D-01  tright=  5. 38D-01 

ttop  =  1. 68D-01  tbot  =  -1. 00D-04 


qeigen(  1)  =  2. 86459D-01  1. 65809D-01 


zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  5. 38D-01  tright=  8. 06D-01 

ttop  =  1. 68D-01  tbot  =  -1.00D-04 
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ro  to 


zeros  rouriu  in  expanded  search  rectangle  defined  by 
tleft=  8. 06D-01  tright=  1. 08D+00 

ttop  =  1. 68D-01  tbot  =  -1.00D-04 


*****  modal  eigenvalues  nrmode=  9  ***** 


mode 

1 

eigenvalue 
-1. 26956D-01 

1. 63761D-01 

theta 

1. 16912D-03 

2. 38566D-03 

atnu  db/km 
4. 8743D+00 

2 

-9. 29119D-02 

1. 25124D-01 

1. 03533D-03 

2. 05835D-03 

3. 7243D+00 

3 

-4. 14322D-02 

8. 94779D-02 

9. 86784D-04 

1. 54436D-03 

2. 6633D+00 

4 

8. 24071D-03 

6. 91738D-02 

1.  15188D-03 

1. 02280D-03 

2. 0589D+00 

5 

5. 55116D-02 

5. 32482D-02 

1. 50185D-03 

6. 03859D-04 

1. 5849D+00 

6 

7. 68261D-02 

6. 70882D-02 

1. 74517D-03 

6. 54734D-04 

1. 9969D+09 

7 

1. 39442D-0 1 

9. 72569D-02 

2. 29575D-03 

7. 21529D-04 

2. 8948D+00 

8 

2. 10278D-01 

1. 35783D-01 

2. 80081D-03 

8. 25696D-04 

4.  0415IJ+00 

9 

2. 86459D-01 

1. 65809D-01 

3. 24286D-03 

8. 70843D-04 

4:  9353D+00 

******** S amp le  printout 


for  9ghz02tn.  pit******** 


3  1  2 

frequency  =  9600. 0000  mhz 


18. 

5 

48.  0 

46.  0 

52.  90 

61.  21 

84.  54 

76.  23 

18. 

5 

48.  0 

53.  0 

69.  58 

73.  08 

67.  86 

64.  36 

27. 

8 

48.  0 

46.  0 

2.  17 

38.  48 

138. 80 

102.48 

27. 

8 

48.  0 

53.  0 

19.  57 

43.  39 

121.  39 

97.  57 

37. 

0 

48.  0 

46.  0 

3.  41 

22.  15 

140. 05 

121.  31 

37. 

0 

48.  0 

53.  0 

2.  69 

25.  91 

140. 77 

117. 55 

"iWfjVVfV 

r*Vi 

'Sample  printout  for 

a  14  m 

evaporation 

duct  input 

file*****' 

9ghzl4m 


0 

9600. 00 

FQMZIN 

0. 000000 

DELFQ 

1 

NFREQ 

0 

MPOL 

2. 100000 

ALOSS 

15. 000000 

SEATMP 

35. 000000 

SEASLT 

1 

IFLAGB 

0. 000000 

AIRTMP 

0. 000000 

RH 

0. 000000 

WGPM3 

0. 250000 

RMSBHT 

25. 000000 

ZTIMT 

0. 000000 

DELZT 

1 

NZT 

56.  5 
58.  6 
56.  5 
58.  6 
56.  5 
58.  6 


126 


3. 0000000 

ZRIN’IT 

7. 0000000 

DELZR 

2 

NZR 

18. 5  COG 

XINIT 

9. 250000 

DELX 

3 

NX 

0. 000000 

ZREF 

20 

NZLAYR 

0.  oooooo 

ail  0] 

0. oooooo 

zim[  Oj 

0. oooooo 

zigab[  0] 

0. 040000 

zif  1] 

-3. 900000 

zim[  1] 

0. oooooo 

zigabf 1] 

0.  100000 

zi[2] 

-5.  340000 

zim[  2] 

0.  OOOOOO 

zigabf  2] 

0.  200000 

zif  3] 

-6.  400000 

zimf  3] 

0.  OOOOOO 

zigabf  3] 

0.  398000 

zif  4] 

-7.  460000 

zimf  4] 

0.  oooooo 

zigabf  4]* 

0.  794000 

zif  5] 

-8. 490000 

zimf  5] 

0. OOOOOO 

zigabf  5] 

1.  585000 

zif  6J 

-9. 470000 

zimf  6] 

0. OOOOOO 

zigabf 6] 

3. 162000 

zif  7] 

-10. 350000 

zimf  7] 

0. OOOOOO 

zigabf  7] 

6. 310000 

zif  8] 

-11. 040000 

zimf  8] 

0. oooooo 

zigabf  8] 

12. 589000 

zif  9] 

-11. 320000 

zimf  9] 

0. oooooo 

zigabf  9] 

14. OOOOOO 

zif  10) 

-11. 330000 

zimf  1 0] 

0. OOOOOO 

zigabf  10] 

25. 119000 

zif  11] 

-10. 870000 

zimf 11] 

0. oooooo 

zigabf  11] 

39. 811000 

zif  12] 

-9. 750000 

zimf 12] 

0. OOOOOO 

zigabf 12] 

50. 119000 

zif  13] 

-8. 820000 

zimf  13] 

0. OOOOOO 

zigabf 13] 

63. 096000 

zif  14] 

-7. 560000 

zimf 14] 

0. oooooo 

zigabf 14] 

79. 433000 

zif  15] 

-5. 880000 

zimf 15] 

0. oooooo 

zigabf 15] 

100. 000000  zi[ 16] 

-3. 670000  zim] 16] 
0.000000  zigab] 16] 
125.  633000  zi]  17] 

-0. 800000  zim[ 17] 
0.000000  zigab]  17] 
158. 4S9000  zi[ 18] 
2.920000  zim[  18] 
0.000000  zigab]  18] 

199. 526000  zi( 19] 
7.690000  zim]  19] 
0.000000  zigab] 19] 

209. 526000  zi[ 20] 
8.870000  zim]  20] 
0.000000  zigab]  20] 


-Wovg ample  printout  for  9ghzl4m.  eig,v,v****** 


9ghzl4m 


1 

9600. 000 

FQMZIN 

0. 000000 

DELFQ 

1 

NFREQ 

0 

MPOL 

2. 100000 

ALOSS 

15. 000000 

SEATMP 

35. 000000 

SEASLT 

1 

IFLAGB 

0. 000000 

A1RTMP 

0. 000000 

RH 

0. 000000 

WGPM3 

0. 250000 

RMSBHT 

25. 000000 

ZTINIT 

0. 000000 

DELZT 

1 

NZT 

3. 000000 

ZRINIT 

7. 000000 

DELZR 

2 

NZR 

18. 500000 

XINIT 

9. 250000 

DELX 

3 

NX 

0. 000000 

ZREF 

20 

NZLAYR 

0. 000000 

zi[  0] 

0. 000000 

zim] 0] 

0. 000000 

zigab] 0] 

0. 040000 

zi[l] 

-3. 900000 

zim]  1] 

0. 000000 

zigab] 1] 

0. 100000 

zi[2] 

-5. 340000 

zim] 2] 

0. 000000 

zigab] 2] 

0. 200000 

zi[3J 

-6. 400000 

zim] 3] 

0. 000000 

zigab] 3] 
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0. 398000 

zi]  4] 

-7. 460000 

zim]  4] 

0. 000000 
n  r-  '  fs  r\  ^ 

zigab] 4] 

\J  .  t  J  \J  -7  \j 

-8. 490000 

zl[  jj 
zim[  5] 

0. 000000 

zigab] 5] 

1. 585000 

zi[  b] 

-9. 470000 

zi’n[  6] 

0. 000000 

zigab] 6] 

3. 162000 

zi[  7] 

-10. 350000 

zim[ 7] 

0. 000000 

zigab]  7] 

6. 310000 

zi[S] 

-11. 040000 

zim[  8] 

0. 000000 

zigab] 8J 

12. 589000 

zi[  9] 

-11. 320000 

zim[ 9] 

0. 000000 

zigab]  9] 

14. 000000 

z  i  [  10] 

-11. 330000 

zim[ 10] 

0. 000000 

zigab] 10] 

25. 119000 

zi[  11] 

-10. 870000 

zim[  11] 

0. 000000 

zigab]  11] 

39. 811000 

zi[  12] 

-9. 750000 

zim[ 12] 

0. 000000 

zigab]  12] 

50. 119000 

zi]  13] 

-8. 820000 

zim[ 13] 

0. 000000 

zigab]  13] 

63. 096000 

zi]  14] 

-7. 560000 

zim] 14] 

0. 000000 

zigab] 14] 

79. 433000 

zi]  15] 

-5. 880000 

zim] 15] 

0. 000000 

zigab] 15] 

100. 000000 

zi[  16] 

-3. 670000 

zim] 16] 

0. 000000 

zigab] 16] 

125. 893000 

zi]  17] 

-0. 800000 

zim] 17] 

0. 000000 

zigab] 17] 

158. 489000 

zi]  18] 

2. 920000 

zim] 18] 

0. 000000 

zigab] 18] 

199. 526000 

zi]  19] 

7. 690000 

zim] 19] 

0. 000000 

zigab] 19] 

209. 526000 

zi]  20] 

8. 870000 

zim] 20] 

0. 000000 

zigab] 20] 

94 

nrmode 

(-1. 14139217438230 3E -001, 2.  381981532549529E-002)  [ 1] 
(-9. 79 147 5 19610 1274E -002, 2.  146753762359050E-002)  [  2] 
( -8. 353759867 17646 3E -002,1.  997 15324143537 IE-002)  (  3] 

(-7. 039783 16 709 507 8E- 00 2,1.  889606906408722E-002)  [4] 
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-S.  814228S37496516E-002 
(-4. 65662S790733621E-002 
(-3. 5853306454522 12E -002 
( -2.  52774 200023 139 '  E- 00 2 
(-1.  5469013631 5-5 011E-002 
(-5. 54195607 1147925E-003 
(3. 551321843705835E-003, 
(1.  27539720058021 0E -002, 
(2. 182707485088627E-002, 
(3.  01 13060101 24 74 SE -002, 
(3.  842803 3420 65 884E -002, 
(4. 694372655547260E-002, 
(5. 485441606530146E-002 , 
(6.  227012336327741E-002, 
(7.  02977 85 01629384E -002, 
(7.  803056175576653E-002, 
(8. 555171 30 117415 5E -002, 
(9. 226312 04896058 6E -002, 
(9.  94029268030099 3E -002, 
(1. 06848500279037  IE-001, 
(1.  138252208S58628E-001, 
(1. 207736300828756E-001, 
( 1. 270094738 8 12559E-001, 
( 1.  3 34 153 34307738 OE -001 , 
(1.403978164901863E-001, 
(1. 467611 45 2238766E -001, 
(1.5340136148S5507E-001, 
(1. 596149004805981E-001, 
(1. 652780520579286E-001, 
(  1.  709514646617846E-001, 
(1. 772646284812436E-001, 
(  1.  836482895835790E-001 , 
(  1.  89172S805  710848E-001, 
(  1.  9515017291 45 819E-001, 
(2.  006781595657978E-001, 
(2. 05956690374031 OE -001, 
(2.  106263141 043572E-001, 
(2. 1575 1875 15 1627 4E -001, 
(2. 213117492755677E-001, 
(2. 27329423807333 6E -001, 
(2.  303987147631 242E -001, 
(2. 3189601987  70044E-001 , 
(2.  34888 16 11845628E-001 , 
(2.  402844396949530E-001, 
(2.  455147471247609E-001 , 
(2.  499562362348160E-001 , 
(2. 545987673232816E-001, 
(2.  600542439477720E-001 , 
(2.  662691013974199E-001, 
(2.  7 264 1 86 1 7 644 147 E -001 , 
(2.  787900254603262E-001 , 
(2. 849476 1645444 7 7E -001 , 
(2.  9 1374792272681 6E -001, 
(2.  980 15 1522337 324E- 001 , 
(3. 042675049575047E-001 , 
(3.  113767357254904E-001 , 


,1 

. 8067654 5 06251 47 E -002) 

[5] 

,1 

..  728118397954079E-002) 

[6] 

,1 

. 6633631812 05 5 10E-002) 

[7] 

,1 

..  628 34 8275801875 E -002) 

[8] 

,1 

. 5716551S6158130E-002) 

[9] 

,1 

. 540239335815212E-002) 

[  io] 

1. 

51885120797 24 15E -002) 

[  11] 

1. 

4605 5 6 29 9546 134E- 00 2) 

[  12] 

1. 

45585433572026 OE -002) 

[  13] 

1. 

4276 131407 12 134E-002) 

[  1^] 

1. 

36589230 666012 6E -002) 

[  15] 

1. 

365320803509313E-002) 

I  16] 

1. 

365 0909 239244 9 3E- 002) 

[  17] 

1. 

336276100865 187E-002) 

[18] 

1. 

2820 89 22 7 6466 39E- 002) 

[  19] 

1. 

285 5 7 22206050 75E- 002) 

[20] 

1. 

287860792288040E-002) 

[21] 

1. 

278978888881959E-002) 

[  22] 

1. 

210154187763285E-002) 

[23| 

1. 

20 84861 02 13981 6E -002) 

[24] 

1. 

19544412 1483925E-002) 

[25] 

1. 

220952877564303E-002) 

[26] 

1. 

22238449 7 7 15 7 26E- 002) 

[27] 

1. 

139133734459426E-002) 

[26] 

1. 

1 1 65 7 1468 6332 9 6E- 00 2) 

129] 

1. 

129961336802001E-002) 

[30] 

1. 

11 65438664 15548E-002) 

[31] 

1. 

136779885 13801 2E -002) 

[32] 

1. 

15315 64 7882978 OE -002) 

[33] 

1.  066116492781759E-002) 

[34] 

1.  02367 186 05 5196 2E -002) 

[35] 

1. 

049763 117 145 163E- 002) 

[36] 

1. 

02577386996145  IE -002) 

[37] 

1. 

007563477610222E-002) 

[38] 

1. 

041321459422345E-002) 

[39] 

1. 

0431721S9566190E-002) 

[40] 

9. 

926428640075490E-003) 

[41] 

8.  7982157 19693539E-003) 

[42] 

9. 

0 1502467 16 11194E -003) 

[43] 

8. 

32107 1087819404E-003) 

[44] 

1. 

18555414 2465963E-003) 

[45] 

8.  46479997 06 3716 6E -003) 

[46] 

7. 

882321847941567E-003) 

[47] 

8.  719487457459410E-003) 

[48] 

9. 

69 2237 12037 6 304E -003) 

[49] 

1.  0 108 88 14700002 9E -002) 

[50] 

1. 

010149239614689E-002) 

[51] 

1. 

049907S56151861E-002) 

[52] 

1. 

044438363485 770E-002) 

[53] 

1. 

092098859083637E-002) 

[54] 

1. 

175 19609 926 180 3E -002) 

[55] 

1. 

227969492073 18SE -002 ) 

[56] 

1. 

24614909459 1 169E-002) 

[57] 

1. 

329 7 305 9431 82 35E -002) 

[58] 

1. 

31965757751941  IF, -002) 

[59] 

1. 

3 1437 7 7026 78 184E -002) 

[60] 

(3. 1878197803 3 51 6 7E-001, 1.  333680648010814E-002)  [61] 
(3. 264578S39391167E-001 , 1.  334739430048439E-002)  [ 62] 
(3. 338 11249 72 164 9 0E -001 , 1 . 464166747417008E-002)  [ 63] 
(3. 41ub4l^ob084o28E-001,l.  49015263S527356E-002)  [ 64] 
( 3. 4S 7 67 7 27 008609  IE -001 , 1. 5244S2539938137E-002)  [ 65] 
(3. 56381 123895335 3E -001,1.  55 135302 1082856F -002)  [ 66] 
(3. 643646265942290E-001 , 1.  538942161 737 147E-002)  [ 67] 
(3. 7297  75181063331E-001, 1.  559989625580370E-002)  [ 68] 
(3. 8 1493108562107 4E -001,1. 61S605966683290E-002)  [ 69] 
(3. 89889 1425944480E-00 1 , 1.  689462812245438E-002)  [ 70' 
(3. 980319033943261E-001, 1.  714047724320524E-002)  [71] 
(4. 0672S3853697590E-001,1.  729643131530S17E-002)  [ 72] 
(4. 15494943360938 2E -00 1,1.  75 1352926449 17 IE-002)  [ 73] 
(4. 246400421526860E-001 , 1.  750551050950182E-002)  [ 74] 
(4. 341S31113214738E-001 ,1.  7890424O1852184E-002)  [ 75] 
(4. 435147531419498E-001,1.  866 150403101799E-002)  [ 76] 
(4. 525289981 946153E-001 , 1.  909787541066863E-002)  [ 77] 
(4. 61863009 84 07780E-001,!. 9 1926073S0S2309E-002 )  [ 78] 
(4. 716181770 152394E-001 , 1.  927318684881558E-002)  [ 79] 
(4. 815020432944205E-001 , 1.  95239 1397029549E -002)  [ SO] 
(4. 917304699196995E-001,1. 972766534431C11E-002)  [ 81] 
(5. 020225407 04978 6E -001, 2.  044077417123271E-002)  [ 82] 
(5. 118252276793764E-001 , 2.  102681780869346E-002)  [ 83] 
(5. 21894868568015 3E- 001, 2.  102179462275660E-002)  [ 84] 
(5. 32481519971042 OE -001, 2. 111602176832850E-002)  [ 85] 
(5. 432604 19 6 7 40 195E -001 , 2. 13684712907547SE-002)  [  86] 
(5. 54181116512579 6E -001, 2. 1 7227635446 1229E-002)  [ 87] 
(5. 651832495410506E-001 ,2.  227384315438282E-002)  [ 88] 
(  5. 75668824 777 7 175E-001, 2.  293313938401881E-002)  [ 89] 
(5. 66653340S029618E-001 ,2.  286336319718634E-002)  [  90] 
(5. 981 376451 733623E-001, 2.  292879  795570507E-002)  [91] 
(6. 097267273995461E-001,2.  326416356305734E-002)  [ 92] 
(6. 2137 12082 142 159E-001 , 2.  370559660462316E-002)  [ 93] 
(6. 330414592077247E-001,2.  426761624643026E-002)  [ 94] 


*****vr-.v*s ample  printout  for 
1 

frequency  =  9600. 0000 

horizontal  polarization 
maximum  mode  attenuation 


9ghz 14m. out******** 

mhz 

rate  =  2. 1000  db/km 


seawater  temperature  =  15. 0000  degrees  Celsius 

seawater  salinity  =  35.0000  grams  salt/kg  seawater 

dielectric  constant  of  seawater  =  80.8869 

conductivity  of  seawater  =  4.6400  si/m 

rms  surface  bump  height  =  . 2500  meters 

iflagb  =  1 

air  temperature  =  . 0000  degrees  Celsius 

relative  himidity  =  . 0000  percent 

liquid  water  concentration  in  air  =  . 0000  grams/meter**3 

m(  . 0000  )  =  . 0000  m-units 

gab(  .0000  )  =  .0000  db/km 
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tropospheric  modified  refractivity  profile 


z 

m(  z) 

gab(z) 

dm/dz 

d( gab) /dz 

meters ) 

( m -units ) 

( db/km) 

(m-units/meter) 

( ( db/km)/meter) 

.  0000 

.  0000 

. 00000D+00 

-97. 5000 

.  OOOOOD+OO 

.  0400 

-3. 9000 

. OOOOOD+OO 

-24. 0000 

.  OOOOOD+OO 

.  1000 

-5.  3400 

. OOOOOD+OO 

-10. 6000 

. OOOOOD+OO 

.  2000 

-6.  4000 

.  OOOOOD+OO 

-5. 3535 

. OOOOOD+OO 

.  3980 

-7. 4600 

. 00000D+00 

-2. 6010 

.  OOOOOD+OO 

.  7940 

-8.  4900 

. 00000D+00 

-1. 2389 

.  OOOOOD+OO 

1.  5850 

-9.  4700 

. 00000D+00 

-.  5580 

. OOOOOD+OO 

3. 1620 

-10.  3500 

. OOOOOD+OO 

-. 2192 

. OOOOOD+OO 

6.  3100 

-11.  0400 

.  OOOOOD+OO 

-. 0446 

. OOOOOD+OO 

12. 5890 

-11. 3200 

. OOOOOD+OO 

-. 0071 

. OOOOOD+OO 

14. 0000 

-11.  3300 

. OOOOOD+OO 

.  0414 

. OOOOOD+OO 

25.  1190 

-10.  8700 

. OOOOOD+OO 

.  0762 

. OOOOOD+OO 

39. 8110 

-9.  7500 

.  OOOOOD+OO 

.  0902 

.  OOOOOD+OO 

50.  1190 

-8.  8200 

. OOOOOD+OO 

.  0971 

. OOOOOD+OO 

63. 0960 

-7.  5600 

. OOOOOD+OO 

.  1028 

. OOOOOD+OO 

79. 4330 

-5.  8800 

. OOOOOD+OO 

.  1075 

. OOOOOD+OO 

100. 0000 

-3. 6700 

. OOOOOD+OO 

.  1108 

. OOOOOD+OO 

125. 8930 

-.  8000 

. OOOOOD+OO 

.  1141 

. OOOOOD+OO 

158. 4890 

2. 9200 

. OOOOOD+OO 

.  1162 

. OOOOOD+OO 

199. 5260 

7.  6900 

. OOOOOD+OO 

.  1180 

. OOOOOD+OO 

209. 5260 

8.  8700 

. OOOOOD+OO 

*****  start 

search  for  modal 

eigenvalues 

tmesh^ 

5. 84278D-04 

to 1—  5. 84278D- 

zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  1. 87D-01  tright=  2. 80D-01 


2. 45D- 

02 

tbot  =  -5. ) 

54D-05 

qeigen( 

1) 

= 

1. 89173D-01 

1. 

02577D-02 

qeigen( 

2) 

= 

2. 10626D-01 

9. 

92643D-03 

qeigenf 

3) 

= 

2. 30399D-01 

1. 

18555D-03 

qeigen( 

= 

2. 34888D-01 

7. 

88232D-03 

qeigenf 

5) 

2. 78790D-01 

1. 

17520D-02 

qeigenf 

6) 

2. 72642D-01 

1. 

09210D-02 

qeigenf 

7) 

= 

2. 66269D-01 

1. 

04444D-02 

qeigenf 

8) 

= 

2. 60054D-01 

1. 

04991D-02 

qeigen( 

9) 

2. 54599D-01 

1. 

01015D-02 

qeigenf 

10) 

= 

2. 49956D-01 

1. 

01089D-02 

qeigenf 

11) 

= 

2. 45515D-01 

9. 

69224D-03 

qeigen( 

12) 

= 

2. 40284D-01 

8. 7 1949D-03 

qeigenf 

13) 

= 

2. 31896D-01 

8. 46480D-03 
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t  lef t= 
ttop  = 


qeigenC 

14) 

= 

2. 27329D-01 

8. 321070-03 

qeigen( 

15) 

= 

2. 213i2D-01 

9. 015020-03 

qeigenC 

16) 

= 

2.  157520-01 

8.  798220-03 

qeiger.( 

17) 

= 

2.  t’59.3 7D-01 

1. 04317D-02 

qeigenC 

lb) 

= 

2. 00678D-01 

1. 04 132D-02 

qeigenC 

19) 

1. 95150D-01 

1. 00756D-02 

id  in  expanded 

search  rectangle  defined  by 

9. 25D- 

02 

tright=  1.  i 

S7D-01 

2. 45D- 

02 

tbot  =  -5. i 

B4D-05 

qe igen( 

1) 

- 

9. 94029D-02 

1. 21015D-02 

qeigen( 

2) 

= 

1. 13825D-01 

1. 19544D-02 

qeigen( 

3) 

= 

1. 70951D-01 

1. 06612D-02 

qeiger.C 

4) 

= 

1. 83648D-0 1 

1. 04976D-02 

qeigenC 

5) 

= 

1. 7  7  265D-0 1 

1. 02367D-02 

qeigenC 

6) 

=: 

1. 65278D-01 

1. 15316D-02 

qeigenC 

7) 

= 

1. 59615D-01 

1.  1367SD-C2 

qeigenC 

£) 

= 

1. 534010-01 

1.  11654D-02 

qeigenC 

9) 

1. 467610-01 

1.  12996D-02 

qeigenC 

10) 

= 

1. 40398D-01 

1.  11657D-02 

qeigent. 

11) 

= 

1. 33415D-01 

1. 1 39 14D-02 

qeigenC 

12) 

= 

1. 27009D-01 

1. 2223SD-02 

qeigenC 

13) 

= 

1. 20774D-01 

1. 220950-02 

qeigenC. 

14) 

= 

1. 06849D-01 

1. 20S49D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tlef t=  0. 00D+00  tr ight=  9. 35D-C2 

ttop  =  2. 450-02  tbot  =  -5. 840-05 


qeigenC 

1) 

= 

3.  55132D-03 

1. 5 1885D-02 

qeigenC 

2) 

= 

3. 8428GD-02 

1. 365890-02 

qeigenC 

3) 

= 

9.  22631D-02 

1. 2789SD-02 

qeigenC 

4) 

= 

8.  55517D-02 

1. 287860-02 

qeigenC 

5) 

= 

7.  80306D-02 

1. 28557D-02 

qeigenC 

6) 

= 

7.  02978D-02 

1. 28209D-02 

qeigenC 

7) 

= 

6.  22701D-02 

1. 33628D-02 

qeigenC 

8) 

=r 

5.  46544D-02 

1. 36509D-02 

qeigenC 

9) 

= 

4.  69437D-02 

1. 36532D-02 

qeigenC 

10) 

= 

3.  01131D-02 

1. 42761D-02 

qeigenC 

ID 

= 

2.  1827  ID-02 

1. 45585D-02 

qeigenC 

12) 

= 

1.  27540D-02 

1. 46056D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tlef t=  -9.35D-02  tright=  0. 00D+00 

ttop  =  2.45D-02  tbot  =  -5. 84D-05 


qeigenC 

1)  = 

-8. 35376D-02 

1. 997 15D-02 

qeigenC 

2)  = 

-4. 656630-02 

1. 728120-02 

qeigenC 

3)  = 

-5. 54196D-03 

1. 540240-02 

13.1 


tleft= 
ttop  = 


qei gen( 

4) 

-1. 54690D-02 

1. 57 166D-02 

qeigenC 

5) 

= 

-2. 52774D-02 

1. 62835D-02 

qeigen( 

6) 

= 

-3. 53533D-02 

1. 66336D-02 

qeigen( 

71 

= 

-5. S1423D-02 

1. 80677D-02 

qeigen( 

8) 

-7. 03978D-02 

1. 88961D-02 

id  in  expanded 
-1. 87D-01 
2.45D-02 

search  rectangle 
tr ight=  -9. 35E 

tbot  =  -5. 84E 

defined  by 
(-02 

1-05 

qeigenC 

1) 

= 

-9. 79 148D-02 

2. 14675D-02 

qe igenC 

2) 

ZZ 

-1. 14139D-01 

2. 3819SU-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  -2. 60D-01  tright=  -1. 87D-01 

ttop  =  2.45D-02  tbot  =  -5. 84D-05 


zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  -3. 74D-01  tright=  -2. 80D-01 

ttop  =  2. 45D-02  tbot  =  -5. 84D-05 


zeros  found  in  expanded  search  rectangle  defined  by 
t left=  2. S0D-01  tright=  3. 74D-01 

ttop  =  2. 45D-02  tbot  =  -5. 84D-05 


qeigenC 

1) 

3.  72978D-01 

1. 

55999D-02 

qeigenC 

2) 

= 

3.  64365D-01 

1. 53894D-02 

qeigenC 

3) 

= 

3.  56381D-01 

1. 

55 135D-02 

qeigenC 

4) 

3.  48768D-01 

‘  1. 

52448D-02 

qeigenC 

5) 

= 

3.  41094D-01 

1. 

49015D-02 

qeigenC 

6) 

= 

3.  3381  ID-01 

1. 46419D-02 

qeigenC 

7) 

= 

3. 26458D-01 

1. 

38474D-02 

qeigenC 

8) 

= 

3. 18782D-01 

1. 33368D-02 

qeigenC 

9) 

= 

3.  11377D-01 

1. 

31438D-02 

qeigenC 

10) 

= 

3.  04288D-01 

1. 

31966D-02 

qeigenC 

ID 

= 

2.  98015D-01 

1. 

32973D-02 

qeigenC 

12) 

= 

2. 9 1375D-01 

1. 

24615D-02 

qeigenC 

13) 

= 

2.  84948D-01 

1. 

22797D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  3. 74D-01  tright=  4. 67D-01 

ttop  =  2.45D-02  tbot  =  -5. 84D-05 


qeigenC  1)=  4.  61863D-01  1. 91926D-02 
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qeigen( 

2) 

~ 

4. 52529D-01 

1. 90979D-02 

qeigen( 

3) 

= 

4.  43515D-01 

1. 86615D-02 

qeigen( 

4) 

= 

4.  34183D-01 

1. 78904D-02 

qeigen( 

5) 

= 

4.  24640D-01 

1. 75055D-02 

qeigen( 

6) 

= 

4.  15495D-01 

1. 75135D-02 

qeigen( 

7) 

= 

4. 06728D-01 

1. 72964D-02 

qeigen( 

6) 

= 

3. 98032D-01 

1. 7 1405D-02 

qeigen( 

9) 

= 

3.  89889D-01 

1. 68946D-02 

qeigen( 

10) 

= 

3.  8 1493D-01 

1. 61861D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
t lef t=  4. 67D-01  tright=  5. 61D-01 

ttop  =  2. 45D-02  tbot  =  -5. 84D-05 


qei gen( 

1) 

= 

5. 54181D-01 

2. 

17228D-02 

qeigen( 

2) 

= 

5.  43260D-01 

2. 

13685D-02 

qeigen( 

3) 

= 

5. 32482D-01 

2. 

1 1 160D-02 

qeigen( 

4) 

= 

5. 21893D-01 

2. 

1C218D-02 

qeigen( 

5) 

= 

5. 1 1825D-0 1 

2. 

1026SD-02 

qeigen( 

6) 

= 

5. 02023D-01 

2, 

0440SD-02 

qeigen( 

7) 

4. 91730D-01 

1. 

97277D-02 

qeigen( 

8) 

= 

4. 81502P-01 

1. 

95239D-02 

qeigen( 

9) 

= 

4.  7 1618D-01 

1. 

92732D-02 

zeros  found  in  expanded  search  rectangle  defined  by 
tleft=  5.61D-01  tright=  6. 54D-01 

ttop  =  2. 45D-02  tbot  =  -5. 84D-05 


qeigen( 

1) 

= 

6.  44509D-01 

2. 47339D-02 

qeigen( 

2) 

= 

6.  33041D-01 

2. 42676D-02 

qeigen( 

3) 

= 

6.  21371D-01 

2. 37056D-02 

qeigen( 

4) 

= 

6. 09727D-01 

2. 32642D-02 

qeigen( 

3) 

5. 98138D-01 

2. 29288D-02 

qeigen( 

b) 

= 

5. 86653D-01 

2. 28634D-02 

qeigen( 

7) 

= 

5.  73869D-01 

2. 29331D-02 

qeigen( 

8) 

= 

5. 65183D-01 

2. 2273SU-02 

zeros  found  in  expanded  search  rectangle  defined  bv 
tleft=  6.54D-01  tright=  7.48D-01 

ttop  =  2. 45D-02  tbot  =  -5. 84D-05 


2  inodes  found  on  same  phase  line 

qeigen(  1)=  6. 68467D-01  2.47500D-02 

qeigen(  2)  =  6. 56153D-01  2.47347D-02 


zeros  found  in  expanded  search  rectangle  defined  by 
tlef t=  7.48D-01  tright=  8.41D-01 


ttop  = 


2. 45D-02 


tbot 


-5. 84D-05 


2  modes  found  on  seme  phase  line 
*****  modal  eigenvalues  nrmode=  94  ***** 

mode  eigenvalue  theta  atnu  db/km 


1 

-1. 

14139D-01 

2. 

38198D-02 

3. 

47000D-04 

3. 

36133D-03 

2. 

0384D+00 

2 

-9. 

79 148D-02 

2. 

14675D-02 

3. 

37466D-04 

3. 

11498D-03 

1. 

8371D+00 

3 

-8. 

35376D-02 

1. 

99715D-02 

3. 

39522D-04 

2. 

88036D-03 

1. 

7091D+00 

4 

*  /  « 

03978D-02 

1. 

8896  ID-02 

3. 

49316D-04 

2. 

64885D-03 

1. 

6170D+00 

5 

-5. 

81423D-02 

1. 

80677D-02 

3. 

66463D-04 

2. 

41421D-03 

1. 

5461D+00 

6 

-4. 

65663D-02 

1. 

72812D-02 

3. 

89814D-04 

2. 

17081D-03 

1. 

478SD+C0 

7 

-3. 

58533D-02 

1. 

66336D-02 

4. 

23955D-04 

1. 

92120D-03 

1 

4234D+00 

8 

_2. 

52774D-02 

1. 62835D-02 

4. 

84349D-04 

1. 

64625D-03 

1. 

3935D+00 

9 

-I. 

54690D-02 

1. 

57 166D-02 

5. 

67769D-04 

1. 

35548D-03 

1. 

3450D+00 

10 

-5. 

54196D-03 

1. 

54024D-02 

7. 

28131D-04 

1. 

03582D-03 

1. 

3181D+00 

11 

3. 

55132D-03 

1. 

51885D-02 

9. 

68348D-04 

7. 

68051D-04 

1. 

2998D+00 

12 

1. 

27540D-02 

1. 

46056D-02 

1. 

25460D-03 

5. 

70059D-04 

1. 

2499D+00 

13 

1827  ID-02 

1. 

45585D-02 

1. 

53413D-03 

4. 

64688D-04 

1. 

2459D+00 

14 

3. 

01131D-02 

1. 42761D-02 

1.  76251D-03 

3. 

96631D-04 

1. 

2217D+00 

15 

3. 

S4280D-02 

1. 

36589D-02 

1. 

96946D-03 

3. 

39607D-04 

1. 

1689D+00 

16 

4. 

69437D-02 

1. 

36532D-02 

2. 

16626D-03 

3. 

08626D-04 

1. 

16S4D+00 

17 

5. 

48544D-02 

1. 

36509D-02 

2. 

33540D-03 

2. 

86226D-04 

1. 

1682D+00 

18 

6. 

22701D-02 

1. 

33628D-02 

2. 

48351D-03 

2. 

63474D-04 

1. 

1435D+00 

19 

7. 

02978D-02 

1. 28209D-02 

2. 

63465D-03 

9 

A.  • 

38289D-04 

1. 0972D+00 

20 

7. 

80306D-02 

1. 

28557D-02 

2. 

77370D-03 

2. 

26958D-04 

1. 

1001D+00 

21 

8. 55517D-02 

1. 

28786D-02 

2. 

9027  ID-03 

2. 

17257D-04 

1.  1021D+00 

22 

9. 

2263  ID-02 

1. 27898D-02 

3. 

01313D-03 

2. 

07852D-04 

1. 

0945D+00 

23 

9. 

94029D-02 

1. 

21015D-02 

3. 

12585D-03 

1. 

89575D-04 

1. 

0356D+00 

24 

1. 

06849D-01 

1. 20849D-02 

3. 

23999D-03 

1. 

82645D-04 

1. 0342D+00 

25 

1. 

13825D-01 

1. 

19544D-02 

3. 

34337D-03 

1. 

75087D-04 

1. 0230D+00 

26 

1. 

20774D-01 

1. 22095D-02 

3. 

44356D-03 

1. 73620D-04 

1. 

0448D+00 

27 

1. 

27009D-01 

1. 

22238D-02 

3. 

53092D-03 

1. 

69523D-04 

1. 

0461D+00 

28 

1. 

33415D-01 

1. 13914D-02 

3. 

61798D-03 

1. 54177D-04 

9. 

7483D-01 

29 

1. 

40398D-01 

1. 

11657D-02 

3. 

7 1 101D-03 

1. 

47335D-04 

9. 

5552D-01 

30 

1. 

46761D-01 

1. 

12996D-02 

3. 

79399D-03 

1. 

45S41D-04 

9. 

6698D-01 

31 

1. 

53401D-01 

1. 

11654D-02 

3. 

87856D-03 

1. 

40966D-04 

9. 

5550D-01 

32 

1. 

59615D-01 

1. 

13678D-02 

3. 

95622D-03 

1. 

40704D-04 

9. 

7282D-01 

33 

1. 

65278D-01 

i 

1. 

15316D-02 

4. 

02569D-03 

1. 40268D-04 

9. 

8683D-01 

34 

1. 

7095  ID-01 

1. 

06612D-02 

4. 

09370D-03 

1. 27526D-04 

9. 

1235D-01 

35 

1. 

77265D-01 

1. 02367D-02 

4. 

16832D-03 

1. 

20257D-04 

8. 7602D-01 

36 

1. 

83648D-01 

1. 

04976D-02 

4. 

24268D-03 

1. 

21161D-04 

8.  9835D-01 

37 

1. 89173D-01 

1. 02577D-02 

4. 

30584D-03 

1. 

16655D-04 

8.  7782D-01 

38 

1. 

95 150D-01 

1. 

00756D-02 

4. 

37319D-03 

1. 

12820D-04 

8. 6224D-01 

39 

2. 

00678D-01 

1. 

04132D-02 

4. 

4347  ID-03 

1. 

14982D-04 

8. 9113D-01 

40 

2. 

05957D-01 

1. 

04317D-02 

4. 

49259D-03 

1. 

13703D-04 

8. 9271D-01 

41 

2. 

10626D-01 

9. 

92643D-03 

4. 

54304D-03 

1. 

06994D-04 

8.  4947D-01 

42 

2. 

15752D-01 

8. 79822D-03 

4. 

59766D-03 

9. 

37065D-05 

7. 

5292D-01 

43 

2. 

21312D-01 

9. 

01502D-03 

4. 

65652D-03 

9. 

48020D-05 

7. 

7148D-01 

44 

2. 

27329D-01 

8.  32107D-03 

4. 

7 1922D-03 

8. 63419D-05 

7. 

1209D-01 
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45 

2. 

303990-01 

1. 

18555D-03 

4.  75019D-03 

1. 22215D-05 

1. 

0 146D-01 

46 

2. 

31896D-01 

8. 46480D-03 

4.  76638D-03 

8. 69643D-05 

7. 

2439D-01 

47 

2. 

348S8D-01 

7. 

8S232D-03 

4.  7969 ID-03 

8. 04647D-05 

6. 

7454D-01 

48 

2. 

40284D-0 1 

8. 

71949D-03 

4.  85181D-03 

8. 80035D-05 

7. 

4619D-01 

49 

2. 

45515D-01 

9. 

69224D-03 

4.  90448D-03 

9. 67707D-05 

8. 

2945D-01 

50 

2. 

499560-01 

1. 

01089D-02 

4.  94869D-03 

1. 00029D-04 

8. 

6509D-01 

51 

2. 

54599D-01 

1. 

01015D-02 

4.  99440D-03 

9. 9041  ID-05 

8. 6445D-01 

52 

2. 

60054D-0 1 

1. 

0499  ID-02 

5.  04766D-03 

1. 01853D-04 

8. 

9848D-01 

53 

2, 

66269D-01 

1. 

04444D-02 

5.  10756D-03 

1. 00 134D-04 

8.  9380D-01 

54 

2. 

72642D-01 

1. 

09210D-02 

5.  16837D-03 

1. 03472D-04 

9. 

345813-01 

55 

o 

78790D-01 

1. 

17520D-02 

5.  22643D-03 

1.  10108D-04 

1. 

0057D+00 

56 

2. 

84948D-01 

1. 

22797D-02 

5.  28389D-03 

1. 13802D-04 

1. 

0509D+00 

57 

2. 

91375D-01 

1. 

24615D-02 

5.  34313D-03 

1. 14206D-04 

1. 

0664D+00 

58 

2. 

98015D-01 

1. 

32973D-02 

5.  40378D-03 

1. 20498D-04 

1. 

1379D+00 

59 

3. 

04288D-01 

1. 

31966D-02 

5.  46027D-03 

1. 18348D-04 

1. 

1293D+00 

60 

3. 

1 1377D-01 

1. 

31438D-02 

5.  52345D-03 

1.  16527D-04 

1. 

12“8D+00 

61 

3. 

18782D-01 

1. 

33368D-02 

5.  58872D-03 

1. 16857D-04 

1. 

1 41313+00 

62 

3. 

26458D-01 

1. 

38474D-02 

5.  65564D-03 

1.  19895D-04 

1. 

1850D+00 

63 

3. 

3381  ID-01 

1. 

46419D-02 

5.  7 1907D-03 

1.  2536813-04 

1. 

2530D+00 

64 

3. 

41094D-01 

1. 

49015D-02 

5.  78111D-03 

1. 26222D-04 

1. 

2752D+00 

65 

3. 

48768D-01 

1. 

524480-02 

5.  84578D-03 

1. 27702D-04 

1. 

3046D+00 

66 

3. 

56381D-01 

1. 

55 135D-02 

5. 90923D-03 

1. 28557D-04 

1. 

3276D+00 

67 

3. 

64365D-01 

1. 

53894D-02 

5. 97497D-03 

1. 26126D-04 

1. 3170D+00 

68 

3. 

72978D-01 

1. 55999D-02 

6.  04515D-03 

1. 26366D-04 

1. 3350D+00 

69 

3. 

81493D-01 

1. 

61861D-02 

6.  11381D-03 

1. 29642D-04 

1. 

3852D+00 

70 

3. 

89889D-01 

1. 

68946D-02 

6. 18078D-03 

1. 33851D-04 

1. 4458D+00 

71 

3. 

98032D-01 

1. 7 1405D-02 

6. 24497D-03 

1. 34403D-04 

1. 

4668D+00 

72 

4. 

0672SD-01 

1. 

72964D-02 

6. 31279D-03 

1. 34169D-04 

1. 

4802D+00 

73 

4. 

1549513-01 

1. 

75 135D-02 

6. 38044D-03 

1.  3441313-04 

1. 

4988D+00 

74 

4. 

24640D-01 

1. 

75055D-02 

6.  45021D-03 

1. 32898D-04 

1. 

4981D+00 

75 

4. 

34183D-01 

1. 

78904D-02 

6.  52228D-03 

1. 34320D-04 

1. 

5310D+00 

76 

4. 

43515D-01 

1. 

86615D-02 

6.  59206D-03 

1. 38626D-04 

1. 

5970D+00 

77 

4. 

52529D-01 

1. 

90979D-02 

6.  65873D-03 

1. 40447D-04 

1. 

6344D+00 

78 

4. 

61863D-01 

1. 

9 1926D-02 

6.  72700D-03 

1. 3971 ID-04 

1. 

6425D+00 

79 

4. 

7 1618D-01 

1. 

92732D-02 

6.  79763D-03 

1. 38840D-04 

1. 

6494D+00 

80 

4. 

81502D-01 

1. 95239D-02 

6.  86847D-03 

1. 39196D-04 

1. 

6708D+00 

81 

4. 

91730D-01 

1. 

97277D-02 

6.  9410113-03 

1. 39 179D-04 

1. 

6683D+00 

82 

5. 

02023D-01 

2. 

04408D-02 

7.  01331D-03 

1. 42723D-04 

1. 

7493D+00 

83 

5. 

1 1S25D-01 

2. 

10268D-02 

7.  08148D-03 

1. 45402D-04 

1. 

799413+00 

84 

5. 

21895D-01 

2. 

10218D-02 

7.  15075D-03 

1. 43959D-04 

1. 7990D+00 

85 

5. 

32482D-01 

0 

t—  • 

1 1 160D-02 

7.  22286D-03 

1. 43161D-04 

1. 

8071D+00 

86 

5. 

43260D-01 

2. 

1 3685D-02 

7.  29558D-03 

1. 43428D-04 

1. 

8287D+00 

87 

5. 

54181D-01 

2. 

17228D-02 

7.  36854D-03 

1. 4436313-04 

1. 

8590D+00 

88 

5. 

65183D-01 

2. 

22738D-02 

7.  44134D-03 

1. 46577D-04 

1. 

9062D+00 

89 

5. 

75869D-01 

2. 

29331D-02 

7.  51138D-03 

1. 49509D-04 

1. 

9626D+00 

90 

5. 

86653D-01 

2. 

28634D-02 

7.  58133D-03 

1. 47679D-04 

1. 

9566D+00 

91 

5. 

98138D-01 

2. 

29288D-02 

7.  65513D-03 

1. 46674D-04 

1. 

9622D+00 

92 

6. 

09727D-01 

2. 

32642D-02 

7.  72892D-03 

1. 47398D-04 

1. 

9909D+00 

93 

6. 

21371D-01 

2. 

37056D-02 

7.  80238D-03 

1. 48781D-04 

2. 

0287D+00 

94 

6. 

3304 ID-01 

2. 

42676D-02 

7.  87532D-03 

1. 50898D-04 

2. 

0768D+00 

******** Sample  printout  for  9ghzl4m. pit******** 
3  12 

frequency  =  9600. 0000  mhz 
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ro  K) 


1 


18.5  25.0 

18.5  25.0 

7.8  25.0 

7.8  25.0 

37.0  25.0 

37.0  25.0 


3.  0 

3.  71 

10.  0 

-8.  31 

3.  0 

1.  66 

10.  0 

-1.  61 

3.  0 

.  08 

10.  0 

.  85 

-.  11 

133. 73 

137. 55 

27.  7 

.  94 

145. 74 

136. 49 

33.  6 

-.  54 

139. 30 

] 41.  50 

27.  7 

.  84 

142.  57 

140. 12 

33.  6 

-.  58 

143. 38 

144. 04 

27.  7 

.  90 

142. 61 

142.56 

33.  6 

13S 
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